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Appendix A 

NUCLEAR WEAPONS EFFECTS 
By J. H. Iverson 

The purpose of this appendix is to describe the derivation of Figs. 

2-1 through 2-8. In all but one case, these figures were based on Ref¬ 

erence 1,* hereinafter referred to as ENW. 

Figures 3.67a and 3.67b of ENW were used in the preparation of Table 
A-l. Heights-of-burst (HOB) were optimized to obtain the maximum range 
for each overpressure shown in the air burst portion of the table. HOB 

equal to zero was used to obtain the contact surface burst ranges. 


Table A-l 

DATA FOR A 1 kt WEAPON YIELD 


Peak 

Overpressure (psi' 


_ Air Burst _ 

Range (ft) HOB ( ft" 


Contact Surface Burst 
Range (ft) 


6 

2000 

917 

1367 

8 

1707 

830 

1133 

10 

1456 

727 

1021 

15 

1194 

650 

818 

20 

945 

614. 

709 

30 

672 

565 

584 

40 

565 

500 

510 

50 

525 

466 

459 


* References are listed at the end of this appendix. 
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Figure 2-1 


The air burst curves were derived as follows: 

1. The data in Table A-l were used in ENW Figure 3.69 to obtain 
positive phase duration on the ground versus overpressure and 
versus dynamic pressure at eight pressures each. 

2. The best least squares curve fit of six curve types was chosen 
and the positive phase duration was cube root scaled to 1 Mt. 

3. Curves for overpressure and dynamic pressure were plotted, each 
versus positive phase duration. The equation of each curve is: 

a. For overpressure: t = 1/(0.282831 + 0.009655p) 

b. For dynamic pressure: t = 3.71018 - 0.013855p 

where t and p are positive phase duration and pressure, 
respectively. 


The surface burst curves were derived as follows: 

1. Ground distances were recorded from ENW Figure 3.69 at the point 
where each positive phase duration curve intersected HOB = 0. 

2. These ground distances versus the positive phase duration values 
for overpressure and dynamic pressure were each subjected to a 
least squares curve fit of six different curves; the best fit 

in each case was: 

a. For overpressure: t = d/(3110.49 + 1.07867d) 

b. For dynamic pressure: t = 0.419881 - (57.8975/d) 

where d is the distance from ground zero for 1 kt. 

3. The data shown in Table A-l were then substituted into these 
equations, which resulted in eight relationships between over¬ 
pressure and positive phase duration and eight relationships 
between dynamic pressure and positive phase duration. 

4. The best least squares curve fit of six curve types was then 
obtained for both overpressure and dynamic pressure versus posi¬ 
tive phase duration. The equations (scaled to 1 Mt) for the 
curves plotted are: 

—0•398083 

a. For overpressure: t = 6.03979p 

b. For dynamic pressure: t = 1/(0.258683 + 0.0016909p) 

where t and p are positive phase duration and pressure, 
respectively. 
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Data obtained from Brode, 2 * Figs. 24 and 25, were also plotted on 
Fig. 2-1. Data points for positive phase duration versus peak overpres¬ 
sure were extracted from Brode for overpressure and dynamic pressure. 

Each set of data was subjected to a least squares curve fit of six dif¬ 
ferent curves; the best fit in each case was: 

-0.430266 

a. For overpressure: t = 5.03615 p 

b. For dynamic pressure: t = 1/(0.262811 + 0.0121185p) 

where t and p are positive phase duration and peak overpressure, respec¬ 
tively. The curves of these two equations were plotted on both parts of 
Fig. 2-1 since Brode indicates that positive phase duration for a given 
peak overpressure is independent of HOB. 


Figures 2-2 Through 2-8 


Figures 2-2 through 2-8 required that air burst slant range be cal¬ 
culated from the data in Table A-l (for a contact surface burst, the 
slant range equals the range). Cube-root scaling provided the scaled, 
slant ranges shown in Table A-2 needed for the derivation of these figures. 

Table A-2 

SLANT RANGES FOR 200 kt AND 1 Mt 


Air Burst Contact Surface Burst 


Overpressure 

200 

kt 

1 

Mt 

200 

kt 

1 

Mt 

(psi) 

Yards 

Miles 

Yards 

Miles 

Yards 

Miles 

Yards 

Miles 

6 

4290 

2.44 



2666 

1.51 

4560 

2.59 

8 

3701 

2.10 



2211 

1.26 

3780 

2.15 

10 

3173 

1.80 

5424 

3.08 

1989 

1.13 

3400 

1.93 

15 

2651 

1.51 

4531 

2. 57 

1597 

0.91 

2730 

1.55 

20 

2197 

1.25 

3756 

2.13 

1381 

0.78 

2360 

1.34 

30 

1712 

0.97 

2926 

1.66 

1141 

0.65 

1950 

1.11 

40 

1471 

0.84 

2515 

1.43 

994 

0.56 

1700 

0.96 

50 

1368 

0.78 

2340 

1.33 

895 

0.51 

1530 

0.87 


* Superscript numerals are related to the reference list at the end of 
this appendix. 
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Figure 2-2 depicts thermal energy as a function of peak overpressure 
To prepare the air burst curves, the Table A-2 values were used in ENW 
Fig. 7.105 to obtain values of radiant exposure (thermal energy). These 
data were then scaled to 200 kt and 1 Mt by multiplying by 200 and 1000, 
respectively (ENW page 364). To prepare the surface burst curves, an 
existing computer program 3 was used to obtain data for plotting. 

Figures 2-4 and 2-7, the free-field gamma radiation intensity curves 
were derived by applying the data in Table A-2 to ENW Figs. 8.27a and 
8.27b. The units of roentgens were changed to rads in accordance with 
note 9, ENW, page 579. 

Figures 2-5 and 2-8 values were obtained as follows: 

1. The integrated neutron flux, obtained by entering ENW Fig. 8.61 
with the data in Table A-2, was multiplied by 1.8 X 10~ 9 rad to 
obtain dose in rads (ENW paragraph 11.90). 

2. These data were scaled to 200 kt and 1 Mt by multiplying by 200 
and 1000, respectively, as described in ENW paragraph 8.61. 

Figure 2-3 is merely a summation of the data found in Figs. 2-4 and 
2-5, and Fig. 2-6 is a similar summation of the data found in Figs. 2-7 
and 2-8. 
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Appendix B 


FIRE HAZARD REDUCTION 
By R. G. Carroll 


Section 1 - Siting for Fire Hazard Reduction 


Introduction 


The incendiary characteristics of a nuclear explosion are associated 
with both thermal and blast effects of a weapon. 

Air burst megaton weapons, under conditions of 6 to 12 miles visibil¬ 
ity, are capable of causing thermal ignitions of newsprint and similar 
tinder materials at a ground radius out to where the blast overpressure 
is approximately 2 psi. The types of fuel susceptible to ignition by 

thermal radiation are thin combustible materials such as newspapers, dry 
grass and leaves, and window-exposed fuels such as curtains and drapes. 

In general, the variables affecting the critical ignition energy, or the 
amount of thermal energy necessary to cause ignition, of tinder fuels are 
thickness, composition, color, and moisture content. Ignition of tinder 
fuels does not necessarily mean that a fire will develop. Other, heavier 
combustible materials must be in proximity to sustain combustion, and this 
condition is prevalent in residential and commercial buildings. There¬ 
fore, it is probable that ignition of window fuels, i.e., drapes, shades, 
and curtains, would determine the range of sustained ignitions and multi¬ 
ple fires. 

Variables that tend to mitigate the effects of the thermal pulse are 
atmospheric obscuration and visibility, number of density of ignition 
points, degree of fuel shielding from thermal radiation, window glass, and 
window screens. Ordinary, clean, single-pane window glass, for example, 
has a thermal transmittance in the range of 0.75. Since thermal radiation 
is emitted within a few seconds (80% of total thermal energy within 8 sec¬ 
onds for a 1 MT air burst 1 ), window glass and screens would survive long 
enough before the arrival of the blast wave to be effective in attenuating 
the thermal radiation. 
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When all of these factors and variables are considered, the net ef¬ 
fect is a decreased probable range of sustained ignition. As an example, 
for a condition of 6 to 12 miles visibility and a 1 MT air burst, the 
ignition ground radius for beige cotton curtains shielded by window glass 
and screens is 2 to 3 miles 6 (8-12 psi range). As another approach, 
under the same conditions, Crowley et al. 7 estimate that for a completely 
exposed dwelling room, the probability of sustained room ignition is 50%. 
Therefore, it can be expected that sustained ignitions sufficient to cause 
rapid build-up of extensive fires would occur within the limits of the 
5 psi radius in areas where buildings and structures are not otherwise 
destroyed or suffering extremely severe damage. 

In the regions beyond the 5 psi perimeter, the density of ignitions 
would diminish with distance. Although large fires, even mass fires or 
conflagrations, may occur in scattered locations, the rate of build-up 
and fire-spread would be slow enough that sufficient time would be avail¬ 
able for occupants of threatened shelters to take remedial action and, if 
necessary, to evacuate, hopefully to other shelters. 

At all overpressure levels out to the 2 psi range, the fire raising 
potential of the thermal pulse would be augmented by blast-induced, sec¬ 
ondary fires. According to Moll, 8 secondary ignitions are estimated to 
have a frequency of 0.006 per 1,000 sf of total floor area that has sus¬ 
tained damage by overpressures of 2 psi or more. This corresponds to about 
one ignition per 100 houses. 

The worst fire conditions for a shelter would be within the 5 psi 
perimeter. Even so, survival in this region is considered probable. Of 
some 280,000 people in the Hamburg firestorm of July 27-28, 1943, 85% 
survived and the majority of casualties were suffered by those in rela¬ 
tively primitive, basement-type shelters. 9 Of those people who occupied 
specially designed bunkers and splinter-proof shelters, 100% survived the 
firestorm. This is not meant .to suggest that firestorms are probable in 
the United States as a result of a nuclear attack. Rodden et al. 10 state 
that the optimum combination of the many conditions necessary to initiate 
a firestorm would be encountered infrequently. 

No matter how fires are started, subsequent development and spread 
would be determined by the amount and distribution of combustible materi¬ 
als in the vicinity. Since measures would be taken to protect the slanted 
shelter building from thermal effects and blast induced fires, the primary 
concern is to minimize the probability of ignitions stemming from burning 
buildings adjacent to the shelter building. Studies have shown that the 
probability of fire-spread depends on building density, or "built-upness," 
of an area and therefore the distance between buildings. 11 
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Given that most buildings would not be equipped with protective win¬ 
dow screens or other reflective devices and the interior fire load of 
combustible furnishings and materials in all buildings would be sufficient 
to sustain combustion after ignition, 4 there is a possibility that struc¬ 
tures and buildings adjacent to a shelter building would be on fire an 
hour or two after an attack. However, even though a shelter building may 
be equipped with thermal protection devices over all window openings, these 
would be blown away by the blast wave immediately subsequent to the thermal 
pulse. The shelter building would therefore be exposed to ignition by 
thermal radiation and convection from adjacent building fires. 


Objective 

The objective of this section is to provide a method for determining 
separation distances to minimize the chances of fire-spread from adjacent 
buildings to the shelter building. The method is basically that developed 
by Margaret Law. 12 The Law method provides separation distances only for 
thermal radiation emanating from the vertical wall openings of adjacent 
burning buildings; therefore, additional methods have been introduced to 
give consideration to the effects of flame and to bending of flames by 
wind. Exposure to the effects of convection currents has not been con¬ 
sidered. At the end of this appendix, there are (l) a reference list, 
including all references used in the preparation of this appendix; (2) a 
bibliography, listing the most pertinent other literature reviewed; and 
(3) Annex A, describing alternative methods for developing building separa¬ 
tion distances. 


Background 


Two major factors enter into consideration when planning building 
separation distances: 

• The level of radiation received from a fire in an adjacent build¬ 
ing. 

• The level of radiation that will ignite materials both on the out¬ 
side of buildings and within rooms exposed by window openings. 

The basis for determining these factors is oven-dry, unpainted wood 
whose behavior has been shown by tests to be representative of most com¬ 
bustible materials. Spontaneous ignition of such wood occurs only at in- 

O 

tensities of radiation above 0.8 cal/cm -sec. In the open, spontaneous 
ignition takes place rather quickly after exposure to radiation and either 
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occurs within about 2 minutes or not at all. 12 Pilot ignition occurs when 
a secondary source of ignition is brought sufficiently near to ignite the 
gases emitted by a heated surface, but for pilot ignition to occur in the 
open, heating times of about 10 minutes are needed. Above an intensity 
of 0.3 cal/cm^-sec, pilot ignition of wood from sparks or flying brands 
can take place. 12 Accordingly, 0.3 cal/cm^-sec is taken as the upper limit 
of radiation exposure for determining separation distances. 

Significant attenuation of incident radiation heat flux through window 
glass has been noted and measured. However, for purposes of this guide, 
no credit for such attenuation can be taken for radiation heat flux ema¬ 
nating from adjacent burning buildings since all windows must be assumed 
to have been destroyed by the blast wave. 


Intensity of Thermal Radiation 

The intensity of thermal radiation, I, emitted by a radiator can be 
found from: 

I = e a T 4 

where 


e = emissivity =1.0 for this guide 
a = Stefan-Boltzman constant 

= 1.36 X 1(T 12 cal/cm 2 -sec 1 -deg.C -4 
T = absolute temp. = (°C + 273°) 

The temperature of a fire depends on the rate of burning, and fires 
may be divided into two types:. 

• Those in which ventilation is restricted and the rate of burning 
depends on window size. 

• Those in which the window area is approximately the same as the 
floor area, where the rate of burning depends on the fire load; 
these are considered to be fully ventilated fires. 

Fires with restricted ventilation are not relevant to this guide, and 
therefore only fully ventilated conditions are considered further. 
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Peak intensities for test fires are shown in Fig. B-l. From this fig¬ 
ure, it is seen that the intensity peaks at about 4 cal/cm 2 -sec for all fire 
loads greater than 12 psf in cubical rooms. In practice, rooms are seldom 
cubical: therefore, Law recommends using 4 cal/cm 2 -sec for all fire levels 
in excess of 5 psf of floor area and 2 cal/cm -sec where fire loads are 
less than 5 psf. 


Configuration Factor 

The intensity of radiation, I, on a receiving surface can be related 
to the intensity, I Q , emitted at the surface of a radiator by the 
expression. 

I = 0I O 

where 0 is the configuration factor. This factor takes into account the 
geometry of the radiator as "seen" by the receiving surface and the sepa¬ 
ration distance between the radiator and the receiving surface. In the 
usual case, many openings in a burning building will be radiating; there¬ 
fore, the above expression for n openings becomes 

I = (0 + 0 _ +. 0)1 

12 no 

o 

Using 0.3 cal/cm -sec as the maximum tolerable I on a receiving surface, 
we have 


2X' 


Thus the separation required becomes the distance where the maximum value 
of \ is less than the ratio 2^.5. ^ Since 0 is proportional to S/C , 

where S = area of radiating surface and C = the separation distance to 
the receiving surface, it is possible, given the dimensions and configura¬ 
tions of windows and other openings, to compute the separation distance 

corresponding to the required value of ^ ^ 0^ . Except for simple radia¬ 
tor geometries, calculation of the magnitude of configuration factors is 
complicated. Nevertheless, circumstances may arise where the criteria 
for employing the simplified Law method to determine separation distance 
cannot be met, or perhaps greater precision in computing distances may be 
required. In this event, the only recourse is to determine separation 
distance as a function of 0 of the radiation surface. Annex A presents 
two ways in which this may be accomplished: (l) a computational method 
and (2) an optical analog method. 
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Law's Method 


On the premise that most cases will be accommodated if it is assumed 
that the radiator and receiving surface are parallel to each other. Law 
has developed a simplified method for determining separation distance. 

This method embodies the concept of an "equivalent radiator," which 
is merely the area on the facade of a building that encloses all facade 
openings. As used herein, facade or wall openings refer to the area of 
openings, such as windows and combustible siding, that are capable of 
transmitting or emitting radiation from the wall of a burning building. 
Separation distances have been calculated for various percentages of wall 
openings and the dimension ratio of the equivalent radiator. These are 
shown in Figs. B-2 and B-3 for fire loads less than, and equal to or 
greater than, 5 psf. The dimension ratio, N, is a function of the length 
and height of the equivalent radiator. If the smaller of the two dimen¬ 
sions is designated "a," then N is the ratio of the larger dimension to 
A- Figure B-4 illustrates the derivation of A and N. 

Entering Fig. B-2 or B-3 with the appropriate N and the percentage 
of wall openings in the equivalent radiator (as shown in Fig. B-4) a value 
of C/A can be read directly, and, C, the separation distance, is then eas¬ 
ily calculated. 

As an example using data from Fig. B-4 and assuming a fire load 
greater than 5 psf: 

• N = 2.0 

• Percentage of openings = 56% 

• Enter Fig. B-2 with N 

• Read up to the 56% curve (interpolate) 

• Then C/A = 1.95 

• C = (l.95)A = (1.95)32 = 62.4 ft. 

Building Facade Openings 

The separation distance, C, is very sensitive to the percentage of 
openings in the adjacent burning buildings. Obviously, the most conserva¬ 
tive approach would be to assume that a given building has upwards of 100% 
openings, and this condition will probably prevail at those ground ranges 
within the 8 psi overpressure region. However there are several reasons 
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FIG. B-2 SEPARATION DISTANCE "C" FOR FIRE LOAD OF 5 PSF OR MORE 12 
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FIG. B-3 SEPARATION DISTANCE "C" FOR FIRE LOAD LESS THAN 5 PSF 12 
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L = 64 ft 
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ELEVATION (Scale: 1" « 20ft) 


• H < L, therefore A = H = 32 ft 

• L = 64 ft 

• N = 64/32 = 2 

• Area ABCD = 2048 sf 

• Area ABCD encloses all openings and is termed the "equivalent radiator." 

• Area of openings (16 x 18) x 9 = 1152 sf 

• Area of openings 1152 

Area of ABCD “ 2048 ~ 56% 



ELEVATION OF EQUIVALENT RADIATOR 

FIG. B-4 DERIVATION OF DIMENSIONS FOR EQUIVALENT RADIATOR 
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40 ft 








for not taking this approach, which would lead to impractically large 
separation distances. 

In this region, blast overpressures would be accompanied by a peak 
wind velocity exceeding 200 mph. Moll 8 cites studies indicating that many 
incipient fires would be extinguished in fine fuels by blast winds with 
a peak velocity greater than 120 feet per second (associated with about 
2.2 psi overpressure) and in thicker fuels at 150 feet per second veloc¬ 
ities (3.0 psi overpressures). In addition, if curtain walls and windows 
are demolished in the region inside the 8 psi perimeter, the extremely 
high blast wind would carry away and scatter the lighter ignition fuels. 
Admittedly, this is intuitive and a generalization, since in some cases 
building fuels may be concentrated within a building by blast and thereby 
aggravate the fire hazard. 

These circumstances suggest that, within the 8 psi region, the error 
caused by selection of building openings as solely the percentage of win¬ 
dow and door openings in a facade, ignoring openings caused by curtainwall 
damage, would be offset in most cases by blast effects on potential igni¬ 
tion points. 

Outside the 8 psi region, out to about 2 psi overpressure, there 
would be regions where framed buildings would also be stripped of curtain- 
walls and completely exposed by blast—i.e., have 100% openings—without 
any mitigating effects of blast on ignition points and building fuels. 

Here again, it is suggested that blast effects on adjacent buildings that 
offer a fire threat to the shelter building can be ignored. The reasons 
lie in the following assumptions made by Law: 

1. For fuel loads of 5 psi or greater, the thermal radiation is 

o 

4 cal/cm -sec. This corresponds to a limiting fire temperature 
of about 1100°C (2000°F). 

2. The emissivity of the radiating building is constant at unity. 

3. Heat is being radiated simultaneously at the peak intensity from 
all openings in a building facade. 

These assumptions are perfectly suited to the conservative measures 
advocated for peacetime fires. However, for purposes of fire protection 
in a nuclear attack environment, they are more restrictive than necessary. 
For example, thermal radiation is directly proportional to the fourth power 
of the absolute temperature of the radiating fire. Variations from the 
limiting fire temperature can cause significant differences in radiation 
intensity. Also, the emissivity of a building fire would generally be some¬ 
what less than unity, and it would infrequently be the case when all wall 
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openings are radiating heat of the same intensity in a multistory building. 
This is so because the combustible contents of buildings and the rate of 
burning will vary from room to room. 

If a mass fire occurred, it is possible that buildings adjacent to the 

shelter building could ignite and present burning facades that are one 

expanse of flame. The Law Method, as cited herein, does not accommodate 

this eventuality. Radiation from such an exposure would exceed the 4 cal/ 
2 

cm -sec used by Law as a maximum. However, these conditions would evolve 
and advance relatively slowly and permit evacuation of an exposed shelter 
building as an ultimate alternative. 

Accordingly, it is recommended that, in the application of the Law 
Method, determination of openings in exposing, adjacent building facades 
be limited to the ratio of the area of windows and combustible siding to 
the area of the facade (this would represent a more hazardous fire situa¬ 
tion than the higher overpressure regions where walls have been destroyed 
but also where fuel loads have been blown away or reduced to debris or 
rubble). 


Recessed Elevations 


If a portion of the radiating building facade is set back, as in 
Fig. B-5, an adjustment in separation distance may be possible. 

Generally two types of recesses will be found: (l) a recess with 
openings on all three sides and (2) a recess with openings only on the 
rear wall. 

With openings on all three sides, the recess will appear as a radi¬ 
ating enclosure. In this case, the area of openings in the recess should 
be totaled and expressed as a percentage of the area of the aperture. If 
the total area of openings in the recess is equal to or greater than the 
area of the aperture, the aperture should be considered as a radiator 
with 100% openings. In either instance, the aperture is then treated 
simply as another element of the equivalent radiator representing the 
facade as a whole, and the corresponding separation distances is found 
as before. Figure B-5 includes an example problem. 

When openings occur only on the rear wall of the recess, a reduction 
in separation distance is possible. First a value of the separation dis¬ 
tance C 2 is found, assuming, as before, that all of the openings within 
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PARTIAL PLAN 

1. Area ABCD = 136 x 64 = 8704 sf 

2. Area of openings in AECG = 1024 sf 

3. Area of openings in FBHD = 1920 sf 

4. Area of openings in recess = 2304 sf 

5. Area of EFGH = 2048 sf which is less than 2304 sf of recess openings 

6. Therefore EFGH = 100% openings 


7. 

8. 
9. 

10 . 

11 . 


Total area of openings = items 2 + 3 + 5 = 4992 sf 
item 7 4002 

Percentage openings = = 57% 


From Fig. B-2, C/A = 2.1 


12. C = A x 2.1 = 64 x 2.1 = 134 ft 


FIG. B-5 SEPARATION DISTANCE FOR RECESSED ELEVATION WITH 
OPENINGS ON ALL RECESS SIDES 12 
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80 ft 







the recess are radiating from the aperture. Then the area of openings in 
the recess are reduced by the factor. 


v c l+r 




where r is the depth of recess. The reduction in the area of openings 
is possible because of the apparent reduction when viewed from a point 
at the separation distance. Figure B-6 illustrates the rationale. With 
the reduction in the area of openings for the recess, a second distance, 
C 2 , is found. An example problem is also given in Fig. B-6. 


Elevation with Setback 


When a part of a building is set back as shown in Fig. B-7, there 
can be a corresponding adjustment in separation distance. First the sep¬ 
aration distance is found assuming the elevation is in one plane. Then 
by constructing an equivalent radiator in the manner shown in Fig. B-7, 
a second separation line is drawn. The final line of separation is then 
traced as illustrated. 

The computation for Fig. B-7 is: 

• 21% openings 

• A = 30 

.* More than 5 psf fuel load 

• N for ABCD = 110/30 =3.67 

• From Fig. B-2, C/A = 1.18 and C = 35.4 ft 

• N for equivalent radiator A'B' = 115/30 = 3.83 

• From Fig. B-2, C/A = 1.20 and C' = 36 ft 

The space saved is hatched in Fig. B-7. 


Elevation with Widely Spaced Openings 

If openings, or a group of openings, in a radiating facade are spaced 
widely apart, a point opposite one may receive a negligible contribution 
of radiation from the other. When this is true, the openings may be con¬ 
sidered separately. 
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• Given a percentage of openings = 40% for areas AEDH, EFHG, and FBGC 

• For the equivalent radiator ABCD, N = 5, A = 20 

• From Fig. B-2, a 40% opening gives C 1 = 44 ft. 


• Reduction factor = 



0.54 


• Therefore EFGH has 0.54 x 40% - 22% opening 

• Then for ABCD overall openings = 29% 

• For N = 5, A = 20, from Fig. B-2, C 2 « 35 ft 


FIG. B-6 SEPARATION DISTANCE FOR RECESSED ELEVATION WITH 
OPENINGS IN REAR WALL ONLY 12 
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Partial 









The separation distance is calculated first for an equivalent radia¬ 
tor enclosing all openings in the facade. Then, if the horizontal dis¬ 
tance between two groups of openings is greater than twice the calculated 
separation distance, they should be treated as separate radiators. An 
example is shown in Fig. B-8. 


Factors Modifying Separation Distances 


Separation distances may be reduced for those facades of shelter 
buildings that have no openings and are constructed of noncombustible mate¬ 
rials. In the extreme, the separation can be zero if a 3-hour fire wall 
exists between the two buildings. 13 The amount of reduction for other 
than fire walls is a matter of judgment. Phung 14 states that an examina¬ 
tion of German cities suffering fire bombings indicates that a 10-foot 
separation between two brick buildings "had about a 50% chance of prevent¬ 
ing fire spread." Presumably the exposed wall had no openings and, also, 
the heights of both buildings were comparable. 

If exposing, adjacent buildings are one story or more taller than the 
shelter building, there would be a danger of exposing the shelter building 
roof to a dangerous level of heat radiation when the calculated separation 
distance is reduced. Obviously, judgment is required in this eventuality. 

If the exposed roof is Class A or B, as defined by the National Fire Code, 15 
this hazard is nominal. However, for Class C, D, E, or F, danger of igni¬ 
tion is present. 


Corrections 

The remainder of this section deals with corrections caused by the 
effects of flames that can be applied to separation distances calculated 
in accordance with the foregoing Law Method. Application of these correc¬ 
tions will serve to increase separation distances to the extent that an 
otherwise eligible building may have to be withdrawn as a slanting candi¬ 
date. Therefore, it is recommended that, as a general rule, flame correc¬ 
tions should be applied only to the extent that candidacy (at least for 
more detailed consideration) is not obviated. The rationale for this rec¬ 
ommendation lies in the fact that 100% protection from firespread, on the 
basis of building separation distances, is neither practical nor possible 
in urban areas (e.g., by flying embers). Also, as previously stated, the 
objective of this appendix is to minimize the chance of firespread on the 
premise that shelter evacuation is reserved as the ultimate alternative 
in those instances where the shelter building is dangerously threatened 
by fire. 
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• For equivalent radiator ABCD, Area = 176 x 24 = 4224 sf 


24, N = 17 = 7 3 


• Area of each opening = 7 x 24 = 168sf , 

• Area of all openings in ABCD = 1008 sf 

» „ 1008 

• Percentage openings = 7 —- = 24% 

4224 

• From Fig. B-2 C/A = 1.45, C = 35 ft 

• Since 80 ft is greater than twice 35 ft consider each area, AECG and FBHD, separately 

• Area AECG has 39% openings, N = 3, A = 24. Then C from Fig. B-2 = 44 ft 

• Area FBHD has 58% openings, N = 1, A = 24. Then C from Fig. B-2 = 35 ft 


FIG. B-8 SEPARATION DISTANCE FOR BUILDING WITH WIDELY SPACED 
OPENINGS 12 
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24 ft 




Adjustment for Window Flames 


The Law Method does not consider the contribution to the total thermal 
radiation from a burning building made by flames emitted from windows and 
other openings. In those instances where an exposing building may have a 
high interior fuel load and a large percentage of wall openings, a fully 
involved fire in the building may present a single expanse of flames over 
100% of the exposing facade. One way in which to account for this circum¬ 
stance is an appropriate increase in the percentage of wall openings in 
the exposing facade. 


Adjustment for Roof Flames 

Adjustment for roof flames must be approached on a case-by-case basis. 
The type of roof construction, its age and condition, the number and type 
of roof openings, and the presence or absence of combustible roof 
structures—such as cooling towers—are all factors that must be consid¬ 
ered. 


It is suggested that recommendations of experienced fire department 
personnel be obtained where the combustibility of roofs may be an issue. 
Their suggestions as to potential flame height may be accounted for by 
increasing the area of a burning facade by an amount equal to average flame 
height multiplied by the horizontal dimension of the facade. The flame 
area would of course be incorporated into the equivalent radiator. 

Cohn 11 recommends roof flame heights as follows: 

• flat or slightly peaked wood roofs 30 ft 

• medium peaked and bow-string roofs 45 ft 

• high peaked roofs 60 ft 


Bending of Flames Due to Wind 

A rational method to account for the increased separation distance 
due to bending of flames has not been developed. However, it has been 
suggested 11 that exposure separation distances should be increased as 
follows: 

• For "light" winds - no increase 

• For "average" winds - add 25 feet to the separation distance 
for light winds 
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For "higher" winds - add 40 feet to the separation distance for 
light winds 


Section 2 - Interior Fires and Associated Biological Hazards 


Background 


Fire hazards and behavior in a shelter building, following an enemy 
attack, must be regarded differently than peacetime fires. During peace¬ 
time, the rate of fire spread and smoke production assumes primary impor¬ 
tance in terms of establishing a safe exit time for evacuation of building 
occupants. Furthermore, building component fire ratings imply that evacua¬ 
tion of the building is possible, organized fire-fighting forces are avail¬ 
able, and an adequate water supply exists. 16 ? 17 

On the other hand, evacuation of a shelter building would be hazardous 
following nuclear attack because of fallout. Since a shelter building fire 
may compel evacuation if left to grow unchecked, shelter occupants would 
have to be prepared to suppress incipient fires throughout the shelter 
building. 


Fire Hazard Factors 


The most pertinent biological fire hazard factors affecting shelter 
occupants are: 

• Heat 

• 0 depletion; CO and CO build-up 

• Smoke 

In considering these factors, the shelter area is assumed to be lo¬ 
cated in the basement of a multistory building for purposes of this guide. 


Heat 


In circumstances where a shelter wall or ceiling is exposed to a 
building fire, the hazard to shelter occupants from heat would depend on 
the fire intensity and duration and the diffusivity and thickness of the 
exposed shelter component. The most hazardous condition would occur when 
the shelter ceiling, functioning as the floor slab of the first story, is 
exposed to a rubble fire. 
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Fuel, or fire, loads in commercial, industrial, and public buildings 
commonly exceed 10 psf in a significant proportion of spaces. 4 > 17 

Leutz 18 states that with a building fire load of approximately 7 psf 
German investigations determined that it is possible for the debris from 
a building fire to maintain a temperature of 300-400°C (572-752°F) on a 
floor slab for periods of up to 5 to 6 hours. Under the assumption that 
temperatures in the debris will obtain at the exposed surface of the floor 
slab, Leutz states that modeling techniques indicate a temperature of 50°C 
(122°F) is reached on the underside of an 8-inch, reinforced concrete floor 
slab within 2 to 3 hours after a temperature of 300-400°C (572-752°) is 
reached and sustained on the upper surface of the slab. 

It is not known if the effect of ventilation on the underside of the 
exposed slab was considered by Leutz. Further, his assumption that tem¬ 
peratures inside a burning debris pile, located on a concrete slab, are 

also developed at the exposed surface of the slab may be overly conserva¬ 
tive and needs to be tested by appropriate experiments. However, under 
these conditions survival of shelter occupants may depend heavily on ade¬ 
quate ventilation or other countermeasures such as removing or cooling off 
the burning debris on the exposed side of the slab. 

Tests have been conducted by the National Bureau of Standards 17 in 
fire resistive buildings using various fire loads to simulate occupancy 
types. The time-temperature curves of these tests are plotted against 
the NBS Standard Time-Temperature Curve and are shown in Fig. B-9. Even 

for the lightest fire load tested, 8.5 psf, the "cooling out” rate after 

the peak of the fire gives a residual room temperature of 700°F after 1-1/2 
hours, and a conservative extrapolation would still show a residual tem¬ 
perature in excess of 500°F after 2 hours.* However, if the burnout of 
a building were not accompanied by collapse of structural components, it 
is unlikely that sufficient heat would be developed in the floor slab to 
jeopardize shelter occupants.. On the other hand, if burnout were accom¬ 
panied by collapse of structural components, and an accumulation of hot 
rubble blanketed the floor slab, countermeasures might have to be taken. 
Rubble fires may smolder for days at a time. A full scale building fire 
test conducted by the National Bureau of Standards 16 measured temperatures 
in the rubble pile that persisted "in the vicinity of 1000°F for a period 
of 2 to 3 days. '‘f 

* These temperatures were not measured at the surface of floor slabs but 
rather were taken at some intermediate level within the rooms involved 
in fire. 

t Tests conducted by the U.S. Forest Service, Dept, of Agriculture, on a 
municipal building registered 1900°F in the rubble pile nearly 24 hours 
after the building was set on fire. 20 Slab surface temperatures were 
not reported. 
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Actual ti me-temperature curves recorded in tests as compared 
with the Standard Time-Temperature Curve. This indicates the relative fire 
severity depending upon different fire loading. The Standard Time-Temperature 
Curve, while used as a convenient measure for general fire testing, actually 
represents a condition of high fire severity met in the early stages of actual fires 
only where combustible materials are of such a character as to favor rapid 
development of high temperatures. 

Reprinted by permission from the Fire Protection Handbook, 12 th edition 
copyrighted by the National Fire Protection Association. 

FIG. B-9 TIME-TEMPERATURE TEST CURVES 
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Table B-l lists equivalent fire severities covering a wide range of 
fire loadings for "offices and light commercial occupancies." 


Table B-l 

ESTIMATED FIRE SEVERITY FOR OFFICES 
AND LIGHT COMMERCIAL OCCUPANCIES 17 
(Data for Fire Resistive Buildings 
with Combustible Furniture and Shelving) 


Combustible Content 
Total, Including Fin¬ 
ish of Floor and Trim 
(psf) 

Heat 

Potential, 

Assumed* 

(Btu/sf) 

Equivalent Fire Severity, Ap¬ 
proximately Equivalent to that 
of Test Under Standard Curve 

for the Following Periods: 


5 


40,000 

30 

min 


10 


80,000 

1 

hr 


15 


120,000 

1-1/2 

hrs 


20 


160,000 

2 

hrs 


30 


240,000 

3 

hrs 


40 


320,000 

4-1/2 

hrs 


50 


380,000 

7 

hrs 


60 


432,000 

8 

hrs 


70 


500,000 

9 

hrs 


* Heat of combustion of contents taken at 800 Btu/lb up to 40 psf, 
7600 Btu/lb for 50 psf, and 7200 Btu/lb for 60 psf and more to 
allow for relatively greater proportion of paper. The weights 
contemplated by the tables are those of ordinary combustible 
materials, such as wood, paper, or textiles. 


Safe exposure periods for humans at various combinations of tempera¬ 
ture and relative humidity are given in Fig. B-10. If a critical tempera¬ 
ture of 50°C (122 °f) is assumed, it is readily seen that even under 
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ENVIRONMENTAL TEMPERATURE 



FIG. B-10 SAFE PERIODS OF EXPOSURE FOR EXTREME ENVIRONMENTAL 
TEMPERATURES 
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moderate conditions of 30% relative humidity, those exposed could not sur¬ 
vive much beyond 3 hours. Of course, higher humidities can be expected 
in shelters, and, consequently, safe exposure periods would be considerably 
less than 3 hours unless special provisions for cooling were undertaken. 

As safe exposure periods are exceeded, the physiological effects are 
disastrous. 19 Body temperatures rapidly rise, metabolism increases, and 
heat stroke quickly follows. 19 


0 9 Depletion; CO 2 and CO Build-up 

Dangerous depletion of oxygen and CO 2 build-up in a shelter in an 
external fire environment would be unlikely. On the other hand, there 
could be a danger from CO poisoning. During the Hamburg fire storm of 
World War II, occupants of bunker type shelters were compelled to shut 
down ventilation plants until the fire storm had subsided to avoid draw¬ 
ing in heat, smoke, and CO. Yet, in spite of gross overcrowding in the 
bunkers and extremely uncomfortable conditions, survival was nearly 100%. 
The bunkers were massive concrete structures and were gasproof, since they 
were fitted with air locks and gasketed doors. Occupants of the relatively 
primitive, neighborhood basement shelters were not so fortunate, and deaths 
from heat, smoke, and CO were common although even here more than 80% of 
shelter occupants survived. 19 > 21 Earp 23 concludes that: 

When the ventilating systems (including filtration) of public 
shelters were worked during incendiary raids and especially during 
fire storms, there was a danger of carbon monoxide being drawn in 
because it is not absorbed by the filters; nor does the gas mask 
protect against it. Hence in the Hamburg raids it was very fortu¬ 
nate for the shelter occupants that so many of the ventilating 
systems were closed down after a short time. Even when the ven¬ 
tilating plant was stopped, carbon monoxide might get in when 
doors were opened to admit people during a raid. 

It can be concluded that there is a real danger of carbon monoxide 
poisoning both in shelters and in the streets. The calm appearance 
of the features and the natural position of many of the corpses 
found in the Hamburg shelters was a strong indication that death 
was due to carbon monoxide poisoning or to heat exhaustion. 

Nevertheless, it is believed that such conditions would be less severe 
in the United States, primarily because building densities in this country 
are substantially less than those in Hamburg during World War II. 
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Fires within a shelter building, however, would present a danger from 
C *2 depletion and CO 2 and CO build-up with the maximum hazard occurring when 
the fire was in the shelter. It is for this reason that shelter occupant 
countermeasures are so important. 

Table B-2 summarizes the biological effects on man of various levels 
of concentration of 0 ^ C0 2 , and CO. It should be noted that the effects 
of a varying proportion of each gas assumes an environment with normal 
concentrations of the other two gases. It is known that 0 depletion when 
accompanied by increased levels of C0 2 and CO has a synergistic effect; 
however, quantification of these aggravated effects is lacking. 23 

Sufficient experimentation to provide an adequate basis for analysis 
of the behavior of fire-generated gases within structures has not been 
performed. However, full scale experimental building burns conducted by 
Waterman 16 do give some specifics in this regard. Some of his conclusions 
are: 

1. Oxygen depletion in an active fire zone of a burning building 
will be reproduced throughout interconnecting spaces on the same 
or higher stories of a relatively tight structure or structural 
portion (unless very large volumes are being considered). Oxygen 

levels as low as 5% were measured (normal air has 20.9% 0 ). 

2 

2. In a burning building shelters remote from the fire zone within 
that building—i.e., not exposed to a direct threat from thermal 
effects and located on the same or higher floor level as the 
fire would be exposed to CO concentrations equal to 75% of those 
recorded in the fire zone. Since CO concentration in the fire 
zone attained a maximum of 3.3%, the shelters were exposed to 
levels up to 2.5% CO, well above the lethal limit of 1.28%. 

3. In fire resistive construction, downward movement of CO, emanating 
from fires on floors above a shelter, will probably not occur. 

4. However, wind pressures can drive fire gases into shelter spaces 
even if fire generated pressures are not sufficient to do so. 
Smoldering fires may produce significant infiltration of toxic 
gases by this means. 

5. Fire gases infiltrating a shelter having otherwise cool floor, 
walls, and ceiling do not significantly heat the shelter air 
space and would represent a hazard only because of toxicity or 
the ignition of combustibles within the shelter. 
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Table B-2 


PHYSIOLOGICAL EFFECTS OF OXYGEN, 
CARBON DIOXIDE, AND CARBON MONOXIDE 24 


_ Effects of Oxygen Deficiency 

Oxygen Content 
of Inhaled Air 

( percent) __ Effects 


20.9% No effects; normal air. 

15 No immediate effect. 

10 Dizziness; shortness of breath; deeper and more 

rapid respiration; quickened pulse, especially on 
exertion. 

7 Stupor sets in. 

5 Minimal concentration compatible with life. 

2-3 Death within one minute. 


_ Effects of Carbon Dioxide 

Carbon Dioxide 
Content of 
Inhaled Air 

(percent) __ Effects 


.04% 

2.0 

4.0 


4.5-5 


7-9 

10-11 


25-30 


No effects; normal air. 

Breathing deeper; tidal volume increased 30%. 

Breathing much deeper; rate slightly quickened; 
considerable discomfort. 

Breathing extremely labored, almost unbearable for 
many individuals. Nausea may occur. 

Limit of tolerance. 

Inability to coordinate; unconsciousness in about 
10 minutes. Symptoms increase, but probably not 
fatal in 1 hour. 

Diminished respiration; fall of blood pressure; 
coma; loss of reflexes; anesthesia. Death 
after some hours. 
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Table B-2 (concluded) 


Carbon Monoxide 
Content of 
Inhaled Air 
(percent) 

0 . 02 % 

0.04 

0.08 

0.16 

0.32 

0.64 

1.28 


Effects of Carbon Monoxide 


Effects 


Possible mild frontal headache after 2 to 3 hours. 

Frontal headache and nausea after 1 to 2 hours. 
Occipital (rear of head) headache after 2-1/2 to 
3-1/2 hours. 

Headache, dizziness and nausea in 3/4 hour. Col¬ 
lapse and possible unconsciousness in 2 hours. 

Headache, dizziness, and nausea in 20 minutes. 
Collapse, unconsciousness, and possible death in 
2 hours. 

Headache and dizziness in 5 to 10 minutes, uncon¬ 
sciousness and danger of death in 30 minutes. 

Headache and dizziness in 1 to 2 minutes, uncon¬ 
sciousness and danger of death in 10 to 15 minutes. 

Immediate effect. Unconsciousness and danger of 
death in 1 to 3 minutes. 
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6 . 


Although much more data are needed on real fire environments 
and movement of fire gases through complex structures, the 
above conclusions lead to a further statement that fires within 
a "tight" envelope, such as a basement containing a shelter, 
cannot be tolerated and must be extinguished unless the barrier 
between the fire and the shelter is perfectly gas tight. 

7. A ventilating system equipped to maintain a slight overpressure 
(1/4 to 1/2 inch of water) in the shelter would normally ensure 
against infiltration of smoke and toxic gases from fires exter¬ 
nal to the shelter. (However, if outside fires were present, 
the ventilating system might have to shut down for a period any¬ 
way - a situation that did occur in the Hamburg fire storm. As 
an alternative to maintaining an in-shelter overpressure, toxic 
gas infiltration must be avoided by taking appropriate measures 
to seal all openings in the shelter.) 


Smoke 


Smoke is the most significant of the lesser known hazard factors in 
a fire environment. Although each year many hundreds of casualties from 
smoke inhalation are reported, Lee 23 states that little is known about 
the actual biological mechanisms involved. This is primarily because 
smoke is such an ill-defined term. The generic term, "aerosol," defines 
smoke, fog, hazes, and fumes as a suspension of particles, liquid or solid, 
in a gaseous phase. However, the chemical constituents of these gaseous 
suspensions are difficult to evaluate. For example, smoke from the same 
sample of fuel will vary in the composition of the products or compounds 
in suspension depending on the rate at which it is heated and burned. 

Although opinions vary, some observers state that true cases of fa¬ 
talities from smoke inhalation are rare and that most casualties ascribed 
to smoke really are because of asphyxiation or CO poisoning. 23 Neverthe¬ 
less, a quantitative evaluation of the threat posed by smoke is not possi¬ 
ble, in contrast with 0^ depletion and CO^ and CO build-up, under the 
emergency conditions that would prevail in a shelter. Therefore, because 
of the high probability that smoke and fumes would be toxic, it must be 
assumed that all smoke and fumes are potentially lethal in the immediate 
environment of a shelter. 25 
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Causes of Fire 


The principal cause of initial building fires during nuclear attack 
would be ignition of interior fuels by radiation from the thermal pulse. 
Following the attack, shelter buildings could be exposed to ignitions stem¬ 
ming from the radiation effects of adjacent burning buildings, flying 
brands, and secondary ignitions caused by blast damage. 

Therefore, by limiting the radiant heat energy that could penetrate 
a building, a significant reduction in ignitions may be obtained. 

Flying brands and embers would pose a continual threat during the 
entire time in shelter. The thermal radiation from a weapon threatens 
only for a few seconds and that from adjacent burning buildings no longer 
than the total burnout time of the buildings, a matter of a few hours to 
a half day under the worst conditions. However, flying brands emanating 
from fires elsewhere in the city could be carried great distances. They 
would become particularly prevalent and dangerous in conjunction with mass 
fires and conflagrations, which invariably develop large quantities of 
brands and embers. Usually the most serious threat is to roofs, but in 
a mass fire environment, wind and convection currents frequently develop 
a horizontal component sufficiently strong to carry brands into window 
openings. Fortunately, the threat from flying brands is amenable to con¬ 
trol by fire watch teams, as discussed in Section 3 of this appendix. 

Building codes commonly require Class A or B roof coverings within 
the fire limits of a city or wherever fire-resistive building construction 
is involved. 17 Class C roof coverings are often specified in some cities 
as a minimum standard. All three classes "possess no flying brand hazard. 117 
Nevertheless, the latent protection offered by these classes of roof cover¬ 
ings could be severely taxed in the vicinity of a mass fire, and monitoring 
and corrective action would have to be taken and should be anticipated by 
shelter personnel. 


Build-up of Fires 


When a sustained ignition occurs in a room, there follows a period 
of heat build-up that eventually leads to total involvement of the room 
in fire. This is called a ' flashover,” and it is accompanied by a sudden 
and drastic surge in temperature. 

The period of time required to reach flashover depends on the con¬ 
tents of the room and its size. The higher the fuel loading and the 
smaller the room, the sooner flashover occurs. 25 In a typical residential 
room, the time elapsing to flashover is about 15 to 20 minutes following 
ignition. 2 6 
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Following room flashover, the fire is capable of penetrating other 
rooms in a building. Horizontal and vertical firespread largely depends 
on a building's geometry and construction. The presence and integrity of 
firedoors and firewalls are important in suppressing the spread of the 
fire. In multistory buildings, the construction of the floors and ceilings 
is also important. Where they are made of noncombustible materials, ver¬ 
tical fire spread will be limited to stairwells, elevator shafts, and 
other openings. Use of noncombustible materials can greatly contribute to 
the effectiveness of firefighting attempts. Other factors that affect 
the spread through the building are the fuel or fire loading, the position 
and compartmentalization of hallways, and the location of the initial 
ignition. 

Under building conditions where the buoyant, hot gases of a fire are 
not confined, fires will usually spread to higher levels. Downward spread 
from the floor of ignition will subsequently occur only if structural 
failures expose lower floors to falling embers and debris. Thus downward 
spread depends on fire loading and structural fire resistance. Provided 
that the fire resistance ratings of structural components conform to the 
criteria in the following section, it is unlikely that downward spread of 
fire will occur . 25 

However, even with good fire resistive construction, fires should 
not be allowed to develop. Possible blast damage to fire doors and the 
threat presented by smoke and toxic gases would demand that fires be sup¬ 
pressed as soon as possible after ignition. Fortunately, studies have 
shown that small self-help teams of untrained men, using portable fire 
extinguishers, are effective in fighting room fires up to a point just 
before flashover. 27 After flashover, the only recourse is to attempt to 
seal off the rooms involved. 


Section 3 - Design Countermeasures 


Objectives 


The objectives of design countermeasures against the threat of fires 
are to: 


• Minimize ignitions of the shelter building by radiation from the 
thermal pulse of a weapon and from adjacent burning buildings. 

• Provide portable fire extinguishers for the use of shelter fire 
teams enabling them to extinguish incipient fires. 
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• Minimize spread of fire by appropriate construction directed to 
containment of fire on the floor of origin. 

• Minimize downward spread of fires from the floor of origin. 

• Prevent spread of smoke, toxic gases, and heat. 

• Restrict shelter building candidates to those occupancies having 
moderate or lesser fire loadings. 

Wherever the following criteria are exceeded by the National Building 
Code 13 for the type of fire-resistive construction and occupancy selected 
for the shelter buildings's primary use, that code or any other applicable 
codes shall govern. 


Radiation 


Exterior wall openings, and roof openings such as skylights, should 
be provided with aluminized blinds or reflective drapes of noncombustible 
composition to provide protection from the thermal pulse of the weapon. 28 } 29 
Ordinary painted metal Venetian blinds, when closed, will also suffice, 
even though the paint will char, provided that they are hung with noncom¬ 
bustible tapes. As a minimum, shielding should be provided on the first 
floor of the shelter building. 

Protection from radiation from adjacent burning buildings should be 
provided by siting to obtain appropriate separation distances as described 
in Section 1 of this appendix. 


Portable Fire Extinguishers 

One 5-gallon stirrup pump type extinguisher with extra length hose 
(at least 6 feet long) should be provided for each 1,250 sf of shelter 
area with a minimum of 4 extinguishers per shelter. Also one 5-gallon 
stirrup pump type extinguisher should be provided for each 5,000 sf of 
nonshelter area. 29 


Structural Criteria 


The minimum fire resistance ratings for structural components of 
shelter buildings should be that specified for "Ordinary Construction 
as defined in Section 707 of the National Building Code, except that 


B-34 


separation distances calculated by the Law Method should govern where ex¬ 
posed facades (or walls) have window openings. Horizontal separation as 
defined by the National Building Code is not the same as the separation 
distance used by Law. Rather, under the code, horizontal separation is 
established for each type of building, then the separation distance be¬ 
tween two buildings becomes the sum of their respective horizontal separa¬ 
tions. 


Table B-3 summarizes the fire resistance ratings for "Ordinary Con¬ 
struction. " 


Smoke and Toxic Gases 


One alternative for protection of the shelter against infiltration 
of smoke and toxic gases would be to provide an independent, emergency 
shelter ventilation system designed to provide a positive pressure in 
shelter spaces of 1/4 to 1/2 inch of water. 25 

The other alternative is to seal the shelter against infiltration. 
Penetrations of the shelter walls and ceilings by cables, pipes, and ducts 
should be grouted or equipped with sealants or gland seals to prevent in¬ 
filtration. Ventilation ducts extending from the normal building ventila¬ 
tion system into the shelter should have vaportight fire dampers installed 
at the point of entrance to the shelter. 17 These should be automatic, 
heat-fused dampers capable of being operated manually by shelter personnel 
from inside the shelter. 

The alternative chosen will depend on cost-effectiveness considera¬ 
tions. 


Combustible Contents 


It is recommended that the maximum combustible fire load for any floor 
level posing a threat to the shelter should not exceed an average of 15 psf. 
Fire loads are composed of combustible furnishings, equipment, and struc¬ 
tural components. Building occupancy generally governs the fire load. 

To select candidate shelter buildings, the proposed occupancy of a building 
can be checked against National Bureau of Standards data on the combustible 
contents of representative buildings by occupancy. 

If candidate buildings are to be equipped with automatic sprinklers, 
to the extent that such sprinklers exceed the minimum requirements of the 
codes 13 , 30 a proportionate increase in the 15 psf maximum average fire 
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Table B-3 


MINIMUM FIRE RESISTANCE 

RATINGS FOR 

ORDINARY CONSTRUCTION 13 



Maximum 



Horizontal 

Percentage of 

Minimum 


Separation 

Windows to 

Fire Resistanc 

Component 

(feet) 

Wall Area 

Rating 

Exterior bearing walls* 

0-3 

0% 

3 hours 


3-20 

20 

2 hours 


20-30 

30 

2 hours 


Over 30 

40 

2 hours 

Exterior nonbearing walls* 

0-3 

0 

3 hours 


3-20 

30 

2 hours 


20-30 

40 

1 hour 


Over 30 

100 

None required 

Interior bearing walls* 



2 hours 

Interior partitions enclos- 



2 hours 


ing elevator shafts, stair¬ 
ways, and so forth, through 
floors of buildings four 
stories or more in height 


Interior partitions enclos¬ 
ing elevator shafts, stair¬ 
ways, and so forth, through 
floors of buildings fewer 
than four stories in height 

Roof structure 


1 hour 


1 hour 


Roof covering 


Class A or B 


Interior finish, maximum flame 
spread rating in nonshelter areas 


75 1 


Interior finish, maximum flame 25^ 

spread rating in shelter areas 


* All exterior walls and interior load-bearing walls will be of approved 
noncombustible materials, 
t Index based on ASTM E-84 Tunnel Test Method. 
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load may be allowed provided that the building sprinkler system is served 
by a gravity tank system that may be expected to survive the air blast 
loading. 


Section 4 - Shelter Occupant Countermeasures 


Fire countermeasures in a shelter building must be viewed as a purely 
local problem. Outside assistance could not be depended on, and shelter 
occupants would have to be prepared to take appropriate self-help measures 
against the threat of fire. Detailed guidance in this regard has been 
developed by IIT Research Institute in a report covering shelters in exist¬ 
ing buildings. 29 The report is equally applicable to shelters in slanted 
buildings, and pertinent paragraphs therefrom have been incorporated in 
this section. 


Organization 


The shelter fire chief should be selected well before an emergency 
period and receive special shelter fire fighting training. He should be 
chosen from the normal peacetime occupants of the building who are famil¬ 
iar with its interior layout, location of utility system controls, and 
characteristics of its occupancy. In addition, he should possess some 
firefighting or fire protection background. He should receive training 
covering all phases of fire prevention, detection, and extinguishment, 
including planning, equipping, and implementation techniques. He should 
also be responsible for the necessary equipment and be familiar with its 
use. In the attack period, he will be charged with staffing, instruction, 
and supervision of firefighting and fire watch teams, as well as other 
personnel engaged in fire prevention activities. 

The successful suppression of all incipient fires that could occur 
within fallout shelter areas and other portions of shelter buildings would 
depend on the availability of sufficient personnel within each shelter 
who possess the ability to handle such fires utilizing only portable ex¬ 
tinguishing appliances. 

Fires of this nature could be effectively controlled by the use of 
two-man self-help teams, with each man equipped with an extinguisher and 
smoke mask. Therefore, the number of shelter firefighting personnel 
should be a multiple of the number of portable extinguishing appliances 
required for any given shelter building. It would appear reasonable to 
use a basis of three firefighters per extinguisher since this would allow 


B-37 


for three 8-hour shifts and would also permit periodic rotation of fire¬ 
fighting personnel where required to limit total radiation exposure. 

All personnel should be given training that includes actual suppression 
of interior fires utilizing portable extinguishing appliances and smoke 
masks. In addition, basic instruction should be provided with reference 
to the types and mode of operation of various types of automatic sprinkler, 
standpipe and hose, and other extinguishing systems to ensure their max¬ 
imum utilization when available in shelter buildings. 


Shelter Building Vulnerability 

Initially, the shelter building should be surveyed to determine its 
vulnerability to fire and potential countermeasures. The survey should 
include: 


• Identification of fire limiting areas capable of unsuppressed 
burnout. 

• Determination of measures necessary to prevent spread of smoke 
and toxic gases into shelter areas. 

• Potential hazard of fire exposure from adjacent buildings. 

• Development of measures necessary to prevent primary ignitions 
from a weapon's thermal pulse. 

• Determination of utility and process systems to be taken out of 
service in time of emergency. 

• Inventory of fire control equipment and facilities. 


The results of these surveys should be used by the fire chiefs and shelter 
managers to develop a fire control operational plan. 

A quick, empirical method of assessing fire vulnerability of a shel¬ 
ter building, in terms of firespread potential, has been developed by 
Varley and Maatman. 29 It is limited to nonoperational, undamaged struc¬ 
tures and may be used for an entire building or to draw conclusions about 
individual floors. Application of this method on an annual basis, and in 
conjunction with recommendations by professional firemen or fire engineers, 
would serve three valuable purposes: (l) forewarning the shelter staff 
of changed building conditions, e.g., fire load and structural modifica¬ 
tions; (2) acquainting shelter staff with the building's firefighting 
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facilities and their condition; and (3) orienting the staff’s thinking 
within a framework of basic fire engineering principles. 

Because the shelter building might sustain damage in a nuclear at¬ 
tack, contingencies might arise that could not be foreseen in advance. 

For example, fire walls, barriers, fire doors, and smoke vents could be 
damaged and rendered ineffective. For this reason, generalizations about 
fire behavior that may be inferred from a fire vulnerability survey should 
be avoided or, at the very least, they should be treated with circumspec¬ 
tion. Some fire behavior generalizations to guard against are: 

• "Downward movement of CO will probably not occur." This does not 
mean that such a condition cannot occur. In tests conducted by 
Waterman 2 6 fires generated gas pressures that caused CO to pene¬ 
trate floor levels below the fire. Also Earp 22 states that: 

"Thus it is quite likely that the occupants of basement shelters 
in Germany were poisoned by carbon monoxide in the smoke penetrat¬ 
ing from the streets or burning buildings above." 

• "a fire in the upper stories of a building will burn upwards" im¬ 
plies relative safety for the occupants on lower stories. In fact, 
if fire barriers are lacking or have been damaged, it is possible 
for fire to spread downwards. 

• "Fire within a shelter space is only a problem of suppression 
since it will be discovered when incipient" is a dogmatic assump¬ 
tion and is highly dangerous. An "incipient" fire inside a shel¬ 
ter space may progress far enough to generate intolerable levels 
of CO and CO if its discovery depends on chance observations by 
occupants. Furthermore, before a small fire can be extinguished, 
the smoke generated may become intolerable. 

Consequently, the shelter manager should base planned remedial action 
on the somewhat subjective conclusion that (l) any fire on the floor lev¬ 
els above a shelter is potentially dangerous and (2) any fire on the same 
level as the shelter, either inside or outside, is potentially lethal. 

The nature of remedial action then becomes obvious - suppress fires as 
they occur. 


Fire Control Plan 


Both preattack and postattack operations must be covered in the fire 
control plan. As a minimum, the preattack plan should include a plan lay¬ 
out of each floor level of the building that indicates: 
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• Boundaries of fire limiting areas and location of fire control 
devices such as smoke vents, fire doors, fire walls, and barriers. 

• Location of firefighting equipment including sprinkler system 
control valves and water flow alarm indicators. 

• Location of all utility services controls. 

• Supplementary instructions as required for operation of the 
building's utility services. 


The postattack portion of the plan should provide for procedures to 
accomplish such tasks as the following: 

• Extinguishment of ignitions that may have been caused by nuclear 
burst. 

• Inspection of the entire shelter area for other possible ignitions. 

• Assignment of shelter firefighting teams to survey the shelter 
building for the purpose of locating and extinguishing fires in 
nonshelter areas. 

• Inspection of the entire shelter building for any other possible 
ignitions. 

• Specific assignment of certain shelter firefighting teams to 
survey adjacent buildings for the purpose of locating and extin¬ 
guishing fires. 

• Where necessary to the long term shelter operation, restoration 
of the various utilities to service after examination of the 
potential safe operability of each. 

• Where necessary, removal of any remaining combustibles to loca¬ 
tions well away from exterior wall openings to reduce the possi¬ 
bility of their subsequent ignition from potential fires burning 
in adjacent buildings. 


• Establishment of a fire watch routine to maintain the entire shel¬ 
ter building under constant surveillance. 


It would be necessary during the postattack operations to coordinate 
the efforts of the various shelter firefighting teams and fire watch per¬ 
sonnel with those of radiation monitoring personnel. This would be 
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especially important for evaluating the potential safe operating periods 
available for any fire control and surveillance efforts that were attempted 
outside the fallout shelter area in the shelter building and in any of the 
surrounding buildings. 


Equipment Required 


• One 5-gallon stirrup pump extinguisher (with a hose at least 

6 feet long) for each 1,250 sf of shelter area with a minimum of 
4 extinguishers per shelter area. 

• Smoke masks (U.S. Bureau of Mines approved type), together with 
spare canisters, provided in number equal to the allowance for 
stirrup pump extinguishers. 

• Gloves and protective clothing for each member of the firefighting 
teams. 

• Such tools as may be available in the shelter building, e.g., 
sledges, saws, chain falls, rope, and wrenches. 


Postattack Fire Surveillance 


With fire control activities limited to the self-help level, it is 
important to establish a continuing fire surveillance program for this 
period. 


Fire Watch in Nonshelter Areas . The problems of maintaining a fire 
watch tour under fallout radiation conditions could be met in two ways. 
With respect to personnel, cumulative dosages could be charted and appro¬ 
priate rotation of fire watches could be accomplished. The other factor 
that would enhance the feasibility of firewatch tours is proper routing 
through irradiated areas. Routes could be designed to follow the most 
shielded portions of each area. 


Fire Watch in Shelter Areas . Discovery of a fire within a shelter 
cannot be left to chance. Because such fires could be potentially lethal, 
their discovery should be followed by immediate extinguishment. Given 
the prevalence of tension and apprehension, it would not be wise to rele¬ 
gate this important task to shelter occupants since it would invite con¬ 
ditions of confusion if not outright panic. 
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In a buttoned-up shelter having no outside source of air, the effects 
of a fire could be critical. Posting of a continuous fire watch is sug¬ 
gested as the only measure that could adequately meet this situation. 

Fire Watch Tour Schedules . Fire watch tours must be scheduled so 
that the time interval between successive tours would not allow a fire to 
develop beyond the stage at which it is extinguishable by self-help teams. 
Therefore, the various functional areas or spaces in the shelter building 
should be classified by fire hazard in terms of the time interval during 
which extinguishment is possible. The time interval for each area, i.e., 
the allowable fire build-up time, then would become the maximum time allow¬ 
able between successive fire watch tours. 

Table B-4 provides allowable fire build-up times for representative 
building functional areas. Since specific buildings may have unique func¬ 
tional conditions, the table is useful only as a baseline for on—site 
evaluation and refinement. 
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Table B-4 


ALLOWABLE FIRE BUILD-UP FOR DORMANT AREAS 
LOCATED IN FIRE RESISTIVE BUILDINGS* 


Area 


Allowable 

Build-up Time 
(min) 

Dining rooms 


40 

Flammable liquids 


None 

Heating and mechanical rm. (no combustibles) 


120 

Janitors' closets 


60 

Kitchens 


12 

Laundry rooms 


30 

Locker rooms (metal lockers) 


60 

Lodging rooms 


30 

Metal working areas 


90 

Motor vehicle areas 


20 

Offices (metal furniture) 


120 

Offices and schools (wood furniture) 


75 

Packing - unpacking 


12 

Paper (rolls, prints, patterns) 


20 

Paper (cutting, gluing, folding) 


12 

Repair shops 


30 

Sales areas 


40 

Sewing areas 


20 

Solid fuel storage (soft coal) 


90 

Solid fuel storage (wood, scraps, dust) 


30 

Stationery supplies 


20 

Theaters and auditoriums (stages without scenery) 


60 

Trash rooms 


12 

Wholesale storage (wood crates, wood patterns) 


40 

Wholesale storage (tight 10' piles, cartons, It. 

crates) 

30 

Wholesale storage (loose 10' piles, cartons, It. 

crates) 

12 

Wood working 


30 

Upholstering 


20 


* Values presume uncongested contents and structural surfaces having 
flame spread not greater than finished wood, 27 
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Annex A - Alternative Methods for 
Developing Building Separation Distances 


Computational 


The basic expression of the configuration factor 0 can be obtained 
from the following diagram wherein dA^ is an elemental radiator, dA2 
is an elemental receiver, and 9 represents the angle between the line- 
distance, C, connecting the centroids of dA^ and dA£ and the normal to 
each area, respectively. 



The intensity of radiation received at any point from a radiator is: 

I = 0I o 

and for our purposes: 

2 

I = 0.3 cal/cm -sec as a maximum 
2 

1=4 cal/cm -sec for a fuel load equal to or greater 
than 5 psf 

-2 

I Q = 2 cal/cm -sec for a fuel load less than 5 psf. 
The energy received by A 2 per unit time from A 1 is: 
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and the intensity of radiation at dA2 from the entire radiator is then: 


I = 


*/ 


COS 0^ COS 0 2 


dA., 


and since 0 = 1/X 0 : 


0 = 


if 


cos 8 1 cos 0 2 


dA- 


0.3 


The configuration factor 0 will be either - “ = 0.075 (fuel load 
0 3 

£5psf) or —— = 0.15 (fuel load < 5 psf) and orientation angles can be 
2 

measured; therefore, the distance C can be computed. 


Optical Analog 31 


Figure B-ll indicates an optical analog that may be employed to 
determine the separation distance C. A model of the building facade 
(the radiator) is constructed of black photographic paper, covered with 
a diffusing material and then attached to the diffusing screen. The dif¬ 
fusing screen is uniformly illuminated. The interior of Section A is 
painted white to obtain high reflectivity. The interior of Section B is 
painted black to prevent reflected light from impinging the photocell. 

The cell has a sensitive area of about 1 cm^ and is used together with a 
galvanometer to measure light intensity. The intensity of light emitted 
can be held constant by employing a constant voltage transformer. A 
calibrated scale attached to the galvanometer can then be used to read 
0 directly. Thus, when the photocell is placed close to the radiating 
area, the configuration factor will be unity, and when placed in the 
position of the receiving element, the configuration factor can be read 
directly from the scale. Should there be reflecting surfaces between 
the radiating surface and the point at which the measurement is to be 
made, these, with suitable reflectivity coefficients, would have to be 
included in the illuminated model. When the model of the radiating area 
is set up, the photo cell is backed off to read 0.15 or .075, depending 
on fuel load, and the proportional distance C is measured directly be¬ 
tween the cell and the radiator. 
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Appendix C 

INITIAL RADIATION CALCULATIONS FOR ENTRANCEWAYS 
By C. K. Wiehle 


Reproduced, without change, from: Wiehle, C. K., 
Shelter Entranceways and Openings , prepared by 
Stanford Research Institute for OCD, final report, 
September 1967. (AD-662 749) 


NOTE: While fully outside entranceways are illus¬ 

trated, the analysis techniques are useful for all 
entranceway configurations. No one type should be 
inferred as the recommended type for all applica¬ 
tions . 
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Appendix C 


INITIAL RADIATION CALCULATIONS FOR ENTRANCEWAYS 
By C. K. Wiehle 


Introduction 


Although there are many unknowns concerning both the attenuation of 
radiation through entranceways and its distribution within the structure, 
the simplified procedures presented in Refs. 5 and 6 were felt to be ade¬ 
quate for the designs prepared here. Basically, there are three problems 
associated with the analysis of an entranceway for the initial nuclear 
radiations. These are the determination of the free-field initial radia¬ 
tion intensities, the maximum permissible dose, and the entranceway attenu¬ 
ation of the radiation. The methods used to calculate these quantities, 
together with an illustrative example, are given in the following subsec¬ 
tions. 


Initial Nuclear Radiation 

The initial nuclear radiation is defined in Ref. 1 as the radiation 
emitted within one minute after the detonation of a nuclear weapon, in¬ 
cluding the prompt radiation. For the purpose of this study, Ref. 1 was 
used to obtain the intensities of the free-field gamma and neutron radia¬ 
tion at the entranceway opening. Because of the yield-range-dose rela¬ 
tionships, it is necessary to specify a yield, range, and HOB (height of 
burst) before the radiation intensity in the free-field can be calculated. 
However, even with precise yield-range data, Ref. 1 cautions that any cal¬ 
culation of radiation exposure dose is subject to wide fluctuation in re¬ 
liability due to variation in weapon design and characteristics. Only 
the 20 psi overpressure level from a 200-kt weapon yield was considered 
for the radiation intensity calculations in this study. 


Maximum Permissible Radiation Dose 


The establishment of a total permissible dose for shelter occupants, 
although beyond the scope of this study, is a primary factor in the de¬ 
sign and cost of entranceways. Even at the modest overpressure level of 
20 psi examined in this study, the initial radiation requirement increased 
the cost of the blast shelter entranceways. The total dose in the shelter 
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is composed of the initial radiation and residual fallout radiation that 
penetrates both the entranceway and the shelter proper. For a specific 
shelter design problem, the relationship between the contributions of the 
radiation from these two sources can be studied to obtain an economical 
compromise for the entranceway and shelter design. However, for a study 
of only the entranceway portion of the shelter system, it is necessary 
first to establish a permissible dose and then to estimate the entrance¬ 
way contribution for both initial and residual radiation. For the pur¬ 
poses of this study, it was sufficient to adopt the criterion presented 
in Ref. 5. Essentially, this limited the total radiation dose in the 
shelter to 40 rads, with the assumption that the permissible dose was 
divided equally at 20 rads each for the entranceway and the shelter pro¬ 
per. Further, the entranceway portion was divided equally between initial 
and fallout radiation, which resulted in a permissible initial radiation 
dose of 10 rads through the entranceway. To adhere to the criterion men¬ 
tioned above for a specific shelter would require assuming an accumulated 
total dose for the free-field fallout radiation and determining its attenu¬ 
ation through the entranceway. However, as noted in the fallout radiation 
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analysis in Appendix 2, for a study of entranceways only, it was believed 
sufficient to determine the PF for the blast entranceways in the usual 
manner (Ref. 20) at a detector point 3 ft inside the shelter entry with¬ 
out calculating the actual radiation dose. 

Another factor that is important in describing the permissible ra¬ 
diation dose received by shelterees is the exposed and absorbed dose of 
nuclear radiation. To describe the effects of nuclear radiation on a 
biological system adequately, it is necessary to express the free-field 
radiation exposure dose as an absorbed dose. This has been accomplished 
by introduction of the rad, which is defined as the absorbed dose of any 
nuclear radiation that is accompanied by the liberation of 100 ergs of 
of energy/gram of absorbing material (Ref. 1). Because of the uncertainty 
of determining the biological effect of nuclear weapons and to simplify 
the radiation analysis, the absorbed dose in rads was used in this study 
to judge the adequacy of the entranceways for providing initial radia¬ 
tion protection. 


Entranceway Radiation Attenuation 


The attenuation of the free-field nuclear radiation intensities 
through an entranceway can conveniently be separated into three phases 
for purposes of analysis: the entrance reduction factor, entranceway 
bend and corridor attenuation, and barrier attenuation. 


* Appendix D herein. 
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Entrance Reduction Factor 

The important factors in determining the initial radiation dose with¬ 
in the entranceway are the orientation of the entrance opening relative 
to the burst point, the entranceway geometry, and the free-field radia¬ 
tion intensity. In general, Refs. 5 and 6 both conclude that the hazard 
from radiation streaming into the entrance opening would be maximum for 
a line-of-sight direction along the longitudinal axis of the entrance 
section between a detector located in the entranceway, usually at the 
first bend, and the burst point. This would result from the fact that the 
higher energy gamma rays and neutron particles are more likely to come 
from a line-of-sight direction than from a scattered direction. Although 
the effect on radiation dose of the angle between the line-of-sight and 
the burst point was not investigated, it was possible to examine limited 
burst positions on the vertical plane through the longitudinal axis of 
the entrance section. The locations were on the line-of-sight from the 
detector approximately along the longitudinal axis of the entrance sec¬ 
tion, on the grazing line-of-sight between the detector and the first step 
of the entranceway, and for a surface burst location. It was found that 
the most critical location was for the grazing line-of-sight. In general, 
as the HOB is increased, the range for a given yield and overpressure al¬ 
so increases to a maximum for some optimum HOB. Since the initial radia¬ 
tion dose is range dependent, it would decrease to a minimum at the opti¬ 
mum HOB for a given yield and overpressure. Although the magnitude of 
the entrance reduction factor is dependent on the angular deviation from 
the line-of-sight orientation, the reduction of this factor for the graz¬ 
ing line-of-sight is more than offset by the increase in free-field ra¬ 
diation intensity due to the decreased range; this was not the case for 
the surface burst location. It may well be that some intermediate angle 
between those selected may produce a greater entrance section radiation 
dose, but it was felt sufficient for the purposes of the entranceway de¬ 
signs and cost estimates in this study to use the grazing line-of-sight 
for the initial radiation analysis. 

To determine the fraction of the free-field radiation that would pen¬ 
etrate the entrance section, it is necessary to compute the solid angle 
fraction subtended by the entrance opening and the detector location. Fig¬ 
ure 1-1, obtained from Ref. 20, can be used to calculate the solid angle 
fraction for rectangular openings. To calculate the entrance reduction 
factor, R^ e for both gamma and neutron radiation, the solid angle fraction 
is used to enter Figure 1-2, which was obtained from Ref. 6. 
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VALUES OF SOLID ANGLE FRACTION 





















































































































































































ENTRANCE REDUCTION FACTOR 



Figure 1 -2 

ENTRANCE REDUCTION FACTOR FOR INITIAL NUCLEAR RADIATION 
VERSUS SOLID ANGLE FRACTION AND BURST POINT 

SOURCE: Ref. 6. 
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Entranceway Bend and Corridor Attenuation 


Gamma Radiation 

As noted in Refs. 15-19, even though considerable analytical and 
experimental effort has been expended in recent years to study the trans¬ 
mission of gamma radiation through tunnels and bends, detailed calculations 
for the initial nuclear gamma radiations are both tedious and of uncertain 
accuracy. For the purposes of this project, the following simplified meth¬ 
od presented in Ref. 5 was therefore felt to be adequate. The reduction 
factor, Rf, for the first 90° bend beyond the entrance section is given by 

R = 0. lu) 
f 1 

where m = solid angle fraction subtended by the 
corridor section at the next point of interest 

For the second and subsequent 90° bends, the reduction factor is given by 

R = 0.5(ju 
f n 

n 

for n = 2, 3... 


Neutron Radiation 


Because of the lack of theoretical and experimental information 
on neutron attenuation in entranceways, it was sufficient for this study 
to use the simplified procedures presented in Refs. 5 and 6 to determine 
the length of corridor required for neutron radiation attenuation. The 
method is based on the concept, of the length of corridor required to re¬ 
duce the neutron dose by one-half. Since neutron attenuation occurs by 
neutron collision with the corridor walls, it is assumed that the attenua¬ 
tion is a function of the average cross-sectional dimensions of the corri¬ 
dor and not of the bends in the corridor. The experimental evidence indi¬ 
cates that the corridor half-length increases with the neutron energy, 
although information is lacking for the higher energy levels associated 
with nuclear weapons. This is accounted for in the method by neglecting 
any neutron radiation attenuation by wall interaction in the first section 
of the entranceway and by assuming that most of the higher energy neutrons 
are not transmitted past the first corridor bend. 
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It is assumed that the neutron half-length for the corridor be¬ 
yond the first bend is given by 

L = 1/2 K (H + W) = 0.366 (H + W) 

1 /Z 

where L = half length of entranceway corridor, ft 
1 /2 

K = 0.732, experimentally determined ratio 
H = height of corridor, ft 


W = width of corridor, ft 


The number of neutron half-lengths for the corridor is given by 
L 


n = 


1/2 


n = number of corridor half-lengths 

L = total length of corridor beyond first bend to 
point of interest, ft 

The reduction factor, R fc , for neutron radiation for the corri¬ 
dor beyond the first bend is given by 

1 

R fc " ( 2 )n 

Because of a lack of adequate information, the method above 
does not specifically include a factor for secondary gamma rays result¬ 
ing from the absorption of thermal neutrons in the corridor walls. In 
Ref. 5, it is stated that the present degree of conservatism in design 
should help reduce this hazard. 


Barrier Attenuation 


Barrier at Entrance 

To determine the barrier reduction factor, B, at the outside en 
trance, Figures 1-3 through 1-5 have been reproduced from Ref. 5. Fig¬ 
ure 1-3 shows the barrier reduction factor in relationship to the mass 
thickness of the shielding material and the angle of incidence for ni¬ 
trogen capture gamma radiation. The use of nitrogen capture gamma 
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BARRIER REDUCTION FACTOR 



Figure 1 -3 

BARRIER REDUCTION FACTOR VERSUS MASS THICKNESS 
FOR NITROGEN-CAPTURE GAMMA RADIATION 

SOURCE: Ref. 5 
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Figure 1 -4 

BARRIER REDUCTION FACTOR VERSUS MASS THICKNESS 
FOR 14 Mev NEUTRONS 

SOURCE: Ref. 5. 
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MASS THICKNESS, X - psf 


Figure 1 -5 

BARRIER REDUCTION FACTOR VERSUS MASS THICKNESS 
FOR 2.5 Mev NEUTRONS 


SOURCE: Ref. 5. 
















radiation for shielding analyses is justified since it is the primary com¬ 
ponent of the initial gamma radiation at the ranges of interest in this 
study, and it is recommended as being on the conservative side in Ref. 1. 

A degree of conservatism for the barrier reduction factor is probably war¬ 
ranted for entranceways due to the previously mentioned neglect of the sec¬ 
ondary gamma rays. Figures 1-4 and 1-5 show the barrier reduction factors 
for fusion yield neutrons (~ 14 Mev) and fission yield neutrons (~ 2.5 Mev), 
respectively. 


Barriers Beyond First Corridor Bend 

Gamma Radiation . The energy level of gamma radiation, which 
has been scattered through an angle of 90°, cannot be greater than 0.51 
Mev, regardless of the initial energy. Therefore, the recommendation in 
Ref. 5 that Figure 1-6 be used for gamma ray barrier shielding beyond the 
first corridor bend, was adopted for this study. 

Neutron Radiation . As recommended in Ref. 5, the reduction 
factors for neutron attenuation through barriers located beyond the first 
bend were obtained from Figure 1-5 for 2.5 Mev neutrons. The use of the 
lower average neutron energy level for interior barriers seems justified 
on the basis of the degradation of the free-field energy level beyond the 
first 90° corridor bend (Refs. 5 and 19). 


Illustrative Example 

To obtain the cost data presented in this report, it was necessary 
to perform an analysis of the attenuation of the initial radiation for 
the six blast entranceways (Types J and K). To demonstrate the method 
outlined in the previous sections, a typical initial radiation analysis 
is presented for entranceway ,J-1 in Figures 1-7 and 1-8. 


Solid Angle Fractions 

The solid angle fractions subtended at the various points shown on 
Figures 1-7 and 1-8 are given in the following tabulation. The values 
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of W, L, and Z were determined graphically as noted on the figures and 
used to obtain the solid angle fractions, w, from Figure 1-1: 


Point 

W 

L 

Z 

e 

n 

U) 

1 

3.00 

3.96 

15.67 

0.76 

7.91 

0.008 

2 

3.46 

6.67 

5.54 

0.52 

1.66 

0.10 

3 

3.54 

6.83 

8.21 

0.52 

2.40 

0.055 

4 

3.42 

6.25 

2.83 

0.55 

0.91 

0.25 


where W, L = dimensions of the entrance opening 
projected on the plane normal to the se¬ 
lected line-of-sight, ft 

Z = distance from detector to plane of W 
and L, ft 

W 


e = — 

L 

2Z 



Assumed Weapon Parameters 


Peak overpressure, P so = 20 psi 
Weapon yield, W = 200 kt 
Wl/3 =5.85 


Free Field Initial Radiation 


As noted previously, the worst case orientation for radiation stream¬ 
ing in the entranceways in this study, was found to be for the grazing 
line-of-sight 0, as shown on Figure 1-7. 


Slant Range 


Knowing the peak overpressure and the angle 0, the horizontal 
range, R, to the shelter entranceway can be determined from Ref. 1 (Para. 
3.67): 


R = 780 W 1 / 3 = (780)(5.85) = 4560 ft 
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and the slant range, R , is 

s 

4560 

R = -- = 4900 ft = 1630 yd 

S COS 0 


Free-Field Initial Radiation 


Using the appropriate figures and scaling factors from Ref. 1 
(Para. 8.27, 8.61, and 11.90), the free-field initial radiation dose can 
be determined for a 200 kt weapon yield at the slant range of the entrance¬ 
way as follows: 

Gamma radiation = 4500 rads 
Neutron radiation = 580 rads 


Entranceway Reduction Factors 


Gamma Radiation 


At Point 1* 


From Figure 1-2, for uu^ = 0.008 and p = (0-0) 


= 7 


R = 0.15 
f e 


At Point 2 


R f " R fe ( °- ll V 


R f = (0.15)(0.1)(0,10) = 0.0015 


At Point 3 


R„ = R (0.1m )B 
f fe 3 


From Figure 1-6,. for X = —(150) = 100 psf 


B = 0.1 

R f = (0.15)(0.1)(0.055)(0.1) = 0.00008 
At Point 4 

R„ = R (0. lu) ) (0. 5cjo ) 
f fe 2 4 

= (0.15) (0.1) (0.10) (0.5) (0. 25) = 0.00019 


* The barrier reduction factor for the wood blast and entry doors is in¬ 
significant and is therefore not included. 
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Neutron Radiation 


At Point 1 

From Figure 1-2, for = 0.008 and g = 7° 

Rfe = 0.08 

At Point 2 

L l/2 = 0.366 (H+W) 

= 0.366 (7.00 + 3.67) = 3.90' 


n 


7.67’ 

_L_ _ 7.67 

^ 1/2 3 • 00 


= 1.97 


*fc - 


(2) n 


= 0.255 


R = R„ X R 
f fe fc 

R = (0.08)(0.255) = 0.020 
f 


At Point 3 


From Figure 1-5, for X = 100 psf and angle of incidence 
90° 

B = 0.35 

Rf = (0.08)(0.255X0.35) = 0.0071 


At Point 4 

L l/ 2 = 0.366 (7.65 + 3.67)* 
= 4.14' 


L 

n 


= 11.84' 

11.84 * 

4.14 


2.86 


* Average for corridor in entranceway and shelter 
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R fc = -(2 \z= °- 138 

R f = R fe X Rfc 

Rf = (0.08)(0.138) = 0.011 

Initial Gamma Radiation Dose 

Through Entranceway 

Dose at Point 2 
R f = 0.0015 

°2 = Rf X Do 

where D 0 = outside radiation dose 
D^ = inside radiation dose 

D 2 = (0.0015) (4500) = 6.8 rads 

Dose at Point 3 

R f = 0.00008 

D 3 = (0.00008)(4500) = 0.4 rad 

Dose at Point 4 

R f = 0.00019 

D 4 = (0.00019)(4500) =0.9 rad 

Through Entrance Roof Slab 
Dose at Point 2 

From Figures 1-7 and 1-8, W = 2.17', L = 3.42', Z = 11.67' 
From Figure 1-1, for e = 0.64, n = 6.82 
UUlA = 0.009 

From Figure 1-2, for u> 1A = 0.009 and p = = 0 = 20° 

R fe =0.07 
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From Figure 1-3, for X = 1/2 150 =75 psf, and angle of 
incidence = 0 = 21° 

B = 0.06 

Rf = (0. lu) 2 ) B 

Rf = (0.07) (0.1)(0.10) (0.06) = 0.000042 
Dg = R f X D Q 

D 2 = (0.00004) (4500) = negligible 

Dose at Point 4 

R f = R fe (0.1 u3 2 )(0.5o) 4 ) B 

R f = (0.07) (0.1)(0.10) (0.5) (.25) (0.06) = 0.000005 
Dj = negligible 


Through Roof Slab over Point 1 


Dose at Point 2 

’ X = <~) + (4.5 X 100) = 525 psf* ' 
z 

From Figure 1-3, for X = 525 psf and angle of incidence 
21 ° 

B«0.0001 
Rf = B (0. lu) 2 )• 

R f « (0.0001) (0.1) (0.10 )«0.000001 
^2 = negligible 

Therefore gamma radiation dose through roof slab over 
point 1 is negligible at all points of interest 


* Assumes an infinite plane barrier with no geometry reduction 
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Through Entranceway Walls 


Dose at Point 4 

X = (1.5 x 150) + (2.5 X 100) = 475 psf 

From Figure 1-3, for X = 475 psf and angle of incidence 
0 - 15° 

D 4 = negligible 

Initial Neutron Radiation Dose 

Through Entranceway 

Dose at Point 2 

R f = 0.020 

°2 = R f X D c 

P 2 = (0.020)(580) = 11.6 rads 

Dose at Point 3 
R f = 0.0071 

D 3 = (0.0071) (580) = 4.1 rads 
Dose at Point 4 

Rf = 0.011 

O 4 = (0.011)(580) = 6.4 rads 

Through Entrance Roof Slab 
Dose at Point 2 

From Figure 1-2, for u> = 0.009 and 0 = 20° 

Rfe = 0.07 

From Figure 1-4, for X = 75 psf and angle of incidence = 
0 = 21 ° 

B = 0.8 
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R f “ R fe X B x Rfc 

Rf = (0.07)(0.8X0.255) = 0.014 

°2 = R f X D 0 

D 2 = (0.014)(580) =8.1 rads 

Dose at Point 3 

R f = R fe X B X R fc x B 

Rf = (0.07)(0.8)(0.255) (0.35) = 0.0050 
D 3 = (0.0050)(580) =2.9 rads 
Dose at Point 4 

R f = R fe X B x Rf c 

Rf = (0.07) (0.8) (0.138) = 0.0077 

D 4 = (0.0077) (580) = 4.5 rads 

Through Roof Slab Over Point 1 
Dose at Point 2 

From Figure 1-4, for X = 525 psf and angle of incidence = 
21 ° 

B = 0.0007 

R = B X R 
f fc 

R f = (0.0007)(0.255) = 0.00018 

D = (0.00018)(580) = negligible 
2 

Therefore, neutron radiation dose through roof slab over 
Point 1 is negligible at all points of interest 
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Through Entranceway Walls 


Dose at Point 4 

From Figure 1-4, for X = 475 psf and angle of incidence = 
0 °* 


B = 0.001 


D 4 = (0.001)(580) =0.6 rads 


Summary 


The tabulation below summarizes the initial radiation dose. 


Type of _Dose in Rads 


Source 


Radiation 

Point 2 

Point 3 

Point 4 

Entranceway 


Gamma 

6.8 

0.4 

0.9 

Entranceway 


Neutron 

11.6 

4. 1 

6.4 

Entrance roof slab 


Gamma 

0 

0 

0 

Entrance roof slab 


Neutron 

8.1 

2.9 

4. 5 

Roof slab at Point 

1 

Gamma 

0 

0 

0 

Roof slab at Point 

1 

Neutron 

0 

0 

0 

Entranceway walls 


Gamma 

— 

— 

_ 

Entranceway walls 


Neutron 

— 

— 

0.6 

Total radiation 






dose 



27 

7 

12 


From the above tabulation, it is apparent that to reduce the 27 rad 
dose at Point 2 to a tolerable level in the shelter would require either 
lengthening the corridor or adding a barrier wall at the shelter entrance. 
The addition of a 6-ft long barrier wall reduces the dose to 7 rads at 
the interior face of the wall and to 12 rads at the entrance to the shel¬ 
ter proper. In view of the uncertainties in calculating the free-field 
nuclear radiation quantities and in the radiation analysis method, the 
entranceway design is considered adequate for the cost analysis performed 


* This assumes the walls are infinite in extent. If the geometry were 
considered, the reduction factor would be decreased. 
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in this study. If desired, it would be relatively inexpensive to reduce 
the radiation dose further by increasing the thickness of the entrance 
slab, by increasing the length of the barrier wall, or by increasing the 
length of the corridor between Points 1 and 2. For instance, increasing 
the thickness of the entrance slab from 6 in. to 9 in. would reduce the 
radiation dose for Point 4 at the shelter entrance from 12 to approxi¬ 
mately 8 rads, but would not significantly increase the total cost of the 
entranceway. 
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Appendix D 


FALLOUT SHIELDING ANALYSIS OF ENTRANCEWAYS 
By H. L. Murphy 


Reproduced, without change, from original publica¬ 
tion as Appendix 2 to: Wiehle, C. K., Shelter 
Entranceways and Openings , prepared by Stanford 
Research Institute for OCD, final report, September 
1967. (AD-662 749) 


NOTE: While fully outside entranceways are illus¬ 

trated, the analysis techniques are useful for all 
entranceway configurations. No one type should be 
inferred as the recommended type for all applica¬ 
tions. 
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FALLOUT SHIELDING ANALYSIS OF ENTRANCEWAYS 
By H. L. Murphy 


The entranceways developed in the project were examined for shield- 
ing capability against fallout contaminant from a nuclear detonation, as 
stated earlier. This somewhat repetitive FSA (fallout shielding analysis) 
work led to an approach that is presented in this appendix. Entrance- 
ways A-2 and C-l are used for demonstration. 

Entranceway pencil drawings to scale were used, since their use per¬ 
mits direct measurement of the many dimensions needed in FSA and thereby 
saves much working time over a method requiring any calculation of dimen¬ 
sions. 

The basic tool was that whose use is currently taught in Office of 
Civil Defense FSA courses for practicing engineers and architects (Ref. 20). 
Therein, a careful comparison of Figure 4.6, used for cleared (or contami¬ 
nated) finite circular areas, with the roof chart (Chart 4) values for a 
zero mass thickness roof showed identical results for the same solid angle 
fractions; however, Chart 4 with its entering argument of solid angle 
fraction (uu) is better suited to entranceway FSA purposes. 

Certain assumptions are necessary in entranceway FSA. Some used in 
this project were: 

1. Contamination factors: Fallout contaminant enters the entrance¬ 
way portal at the "standard" (outdoor) rate per unit (horizontally 
projected) area, and this is then spread uniformly over all hori¬ 
zontal areas of the entranceway but not beyond the first turn. 
Using entranceway A~2 (see Figure 2-1) as an example, the con¬ 
tamination factor would be the area of the portal (3.67 X 9.5 ft) 
divided by the total area of steps and particulate trap under the 
grating (3.67 X 19.42 ft), or 0.489. To assume that the contami¬ 
nant would only fall straight down into the entranceway would be 
nonconservative in PF (protection factor) results and would be 
assuming zero wind conditions. At the other end, it might be 
assumed that most if not all of the contaminant would accumulate 
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on the step at the level of the grating and on the trap under 
the grating. This assumption is conservative but was not used 
because shelter occupants could readily sweep the step contami¬ 
nant through the grating or flush down both step and trap to get 
rid of the contaminant entirely. 

This assumption was modified in entranceways such as B-5 as fol¬ 
lows: Each side of the entranceway was separately handled as 
described for entranceway A-2 down to and including the grating, 
resulting in a doubling of the assumed contamination under the 
grating. 

Another contamination factor was used in dealing with the steps 
located above the horizontal plane of the detector; it is dis¬ 
cussed in the Worked-Out Example 1. 

2. Steel gratings hold no falling contaminant and provide negligible 
shielding. 

3. Doors of ordinary wood or hollow metal construction provide neg¬ 
ligible shielding. 

4. Detector location is a point although the human body is not. 

5. For stairs located above the detector horizontal plane, neither 
the stairs nor stair nosings, if any, provide significant shield¬ 
ing. This assumption is demonstrated by the examples. 

6 . Where a contaminated area could both be "seen" (at least partially) 
by the final detector location and was in a tunnel-like location 
for example, the trap under the grating, entranceway A~2—contribu¬ 
tions were included for both (a) the contaminated area directly 

to the final detector, location and (b) for the same contaminant 
being "seen" at a detector location directly above the contami¬ 
nant, and then through a right-angle (tunnel) turn to the final 
detector location. This assumption was termed the "tunnel effect" 
for use in this study. It was sometimes applied as just described; 
sometimes only the larger of the two contributions (a) and (b) was 
used; or, for a relatively open entranceway (for example, entrance¬ 
way D“4), only contribution (b) was used. The rationale was based 
on the likelihood of gamma radiation scattering off walls in a 
somewhat long and narrow, relatively confined, tunnel-like situa¬ 
tion. 

Entranceway A-2 fallout shielding analysis was made using a pencil 
scale drawing illustrated by Figure 2-1. The computational steps shown 
in Table 2-1 are explained following the table. 
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Table 2-1 


COMPUTATIONS FOR WORKED-OUT EXAMPLE 1 



w 

w 

L 

Z 

U) 

Multipliers 


c 

o 

1. 

U) 

1 

3. 67 

18 

9.3 

.0860 

-0.5/0.489 


-.0036 

2. 


3. 67 

35. 2 

9.3 

. 1095 

0.5/0.489 


.0046 

3. 

(0 

3 

3.67 

10.4 

3.9 

. 2213 

-0.5 X 0.85 


-.0186 

4. 

U) 

4 

3. 67 

40.6 

3.9 

. 2746 

0.5 X 0.85 


.0238 

5. 

(JU 

5 

3.67 

8.3 

1. 1 

. 6223 

-0.5 


-.0791 

6. 

UU 

6 

3.67 

13 

1. 1 

. 6416 

0.5 


.0827 

7. 

U) 

7 

3. 67 

3. 67 

3 

. 1755 

-0.5 


-.0169 

8. 

CD 

8 

3.67 

8.3 

3 

. 2780 

0. 5 


.0284 

9. 

U) 

Q 

3. 67 

3.67 

4 

. 1112 

1 


.0207 

10. 







C = 
l 

.0472 

11. 

“>10 

3.67 

7 

7 

.0723 

°- 2C i“>io 


.0007 

12. 

“’ll 

3.67 

18. 7 

4 

. 2505 

-0.5 


-.0252 

13. 

“>12 

3. 67 

21.33 

4 

. 2553 

0.5 


.0257 








C = 
t 

.0012 X 0.489 








= 

.0006 








PF KS 

1700 


Lines 1 2: Core and peripheral area of a decontaminated, zero mass 
roof (Chart 10, Case 3). Dividing this one item by the contami¬ 
nation factor permits applying the factor to the overall total 
contribution, rather than to each step, because it applies to all 
other contributions except this one. 

Lines 3-on: Roof contributions all taken from Chart 4, zero mass 
thickness. 
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Lines 3-4: Steps changed to a plane, which is then treated as if 

horizontal; second contamination factor enters here, in that the 
plane used has a greater area than the horizontal area of the 
steps, thus the 0.85 multiplier. 

Lines 5-8: Steps below detector horizontal plane, simplified into 
two steps. 

Line 9: Particulate trap contribution to detector location 1. 

Line 11: Converts total contribution to detector location 1 through 
a turn and solid angle fraction into a contribution to detector 
location 2 . 

Lines 12-13: Contribution from portion of trap seen from detector 
location 2 . 

An entranceway C-l FSA was made, also using a pencil scale drawing 
illustrated by Figure 2-2. The computational steps shown in Table 2-2 are 
explained below. 


Table 2-2 

COMPUTATIONS FOR WORKED-OUT EXAMPLE 2 



U) 

w 

L 

Z 

(JU 

Multipliers 


V-/ 

O 

1 . 

‘“l 

3.7 

4 

9 

.0278 

+0.5 [X 0 =75] 


.0011 

2 . 

U>3 

3. 7 

23. 5 

9 

. 10 22 

+0.5 [X 0 = 75] 


.0040 

3. 







C ! = 

.0051 

4. 


3 

6 . 8 

3. 5 

. 1770 

0 . 2 ^ 11)3 

c t a = 

.0002 








pf a~ 

5000 

5. 

“>4 

3.67 

12 

6.4 

. 1207 

1 - 1 

O 

O 

rH 

II 

O 

X 

1 _ 1 

Ifl 

£ 

0 


.0006 

6 . 


2.5 

6 . 67 

3. 5 

. 1491 



— 

7. 

^6 

3.67 

17. 2 

9.3 

.0839 

0.5 [X Q = 100] 


.0021 


Cm = .0027 X 0.571 
= .0015 

PF b « 670 
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ENTR. 

STAIR 

DOOR 

CAPACITY 


NO. 

WIDTH 

W 

SIZE 

D 

PERS. MIN. 

STAIR 

C-1 

3 ' -8 

3'-0 . 6'-8 

90 

TREAD 


Figure 2-2 

WORKED-OUT EXAMPLE 2 







































Door A 


General: Radiation from contamination on stairs must come through 

the concrete steps, a 6-in. concrete wall acting like an interior 
partition, and the doorway restriction, or must be reflected 
through three turns, in either case making this contribution neg¬ 
ligible. Similar thinking may be applied to the trap under the 
grating. Thus, the only contribution left to be considered is 
that from the contaminated roof above the room where detector 
location 1 is shown, an ordinary off-center detector location 
problem. Contributions were all taken from Chart 4. 

Lines 1-2: Roof above off-center detector location 1 (ou and ud 2 not 
shown in Figure 2-2). 

Line 4: Turn and solid angle fraction from 1 to 2. 


Door B 


General: Stairs simplified as in Worked-Out Example 1. Because of 

tunnel effect, contributions taken directly to detector loca¬ 
tion 4, and via 3 to 4; that direct to 4 restricted to portion 
seen through doorway, because any other must come through con¬ 
crete stairs (used 8-in. or 100 psf, for mass) and through 6-in. 
more concrete in lintel, thereby becoming negligible. Any con¬ 
tribution from trap under grating is negligible because of both 
geometry and barriers. Contribution from roof over detector 
location 1 must come through several heavy barriers, so it is 
negligible. Contamination factor for steps and trap under grat¬ 
ing amounts to 9.5/13.17 and must be combined with a factor for 
the simplified stairs amounting to 9.5/12, giving a combined 
factor of 0.571. Contributions were all taken from Chart 4. 

Line 5: First right-angle turn is 0.2 and no turn is 1; for this par 
tial turn, 0.7 was used. 

Line 7: Ray from stairs to detector location 4 treated as normal to 
the simplified stairs, thus solution is for a half-roof. 

Some of the values in Tables 2-1 and 2-2 are shown to an unwarranted 
number of decimal places. They are shown that way, however, because 
(1) all solid angle fractions and all zero mass roof contributions were 
computer-calculated so that no extra work was involved and (2) such might 
be useful to anyone trying to check through the computations, for under¬ 
standing or other reasons. 
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Near the end of this work, two simplifications used (as in the first 
example problem solution herein) were subjected to more detailed examina¬ 
tion. Computer programs were used for computing solid angle fractions and 
for computing radiation contributions, both to save professional working 
time and, more specifically, to carry values out to several decimal places, 
take differences, and then round off the total to fewer decimal places, all 
aimed at reducing round-off values wherein "roof" cores were continually 
subtracted from larger roof" values. Simplifications examined were: 

1. Simplification of the stairs located below the detector horizontal 
plane to fewer stairs. In Figure 2-1, the total contribution from 
the actual steps below the detector horizontal plane, including 
the step at the level of the grating, amounted to 0.0164.* In 
Worked-Out Example 1 the contribution based on the simplification 
amounted to 0.0151,* or an error of -8 percent in this one item. 

In carrying this nonconservative error forward, however, the tun¬ 
nel turn reduces its effect so much as to leave the overall total 
contribution unchanged. 

2. Simplification of the stairs located above the detector horizontal 
plane to a sloping, plane, smooth, contaminated surface, then 
analyzed by rotating the plane of the detector. Each step in Fig¬ 
ure 2-1 was analyzed by considering it in five "roof" increments 
as indicated in the sketch, but the same Z distance was used for 

each step, as shown by the 
sketch arrow. The mass 
thickness of each increment 
was determined from its 
scaled thickness (on a scale 
drawing at l"=0.1'), using 
reinforced concrete at 
150 pcf. With a core and 
peripheral roof calculation 
for each of five increments 
and for each of 15 steps, 

150 roof contributions were 
calculated. In Worked-Out 
Example 1, the contribution 
from the simplification of 
the stairs above the detector 
horizontal plane was 0.0052.* 
The results from the 150 roof 
calculations was 0.00524,* 
indicating an error of 
< 1 percent. Table 2-3 gives 

* Value before applying general contamination factor, 0.489. 
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values obtained for stair widths of 22, 44, and 66 in* , and shows gross 
errors ranging from -9 percent to +10 percent in using the sloping 
plane approximation. Such errors would be indeed negligible after the 
radiation contribution values were multiplied by 0.489 (contamination 
factor), 0.2 (first tunnel turn), and 0.0723 (tunnel solid angle fraction), 
or a combined multiplier of 0.00707. Unfortunately, time was not avail¬ 
able for detailed calculations varying other parameters than stair width. 


Table 2-3 

CONTRIBUTIONS FOR VARIOUS STAIR WIDTHS 


Detailed 

Computations 


Stair 

Width 

C 

T 

9.5"' 

C 

T 

. 6" . 

22 " 

.00166 

.00263 

44" 

.00331 

.00524 

66 " 

.00487 

.00771 


Sloping Plane 
Approximation 


C 0 

o 

.85 C 

o 

Error 

.0028 

.0024 

-9% 

.0061 

.0052 

-1% 

.0100 

.0085 

+ 10% 


The matter of stair nosings was reviewed briefly during the work. 
Nosings are required by some building codes, for the particular tread- 
rise values used. One standard nosing detail, shown in Ref. 38, is a sim¬ 
ple one requiring only the sloping of the riser, the effect of which is to 
increase the tread by 1 in. Such a change would have only a trivial effect 
on the computations herein; e.g., use of this nosing in the examples in 
this appendix would mean only that each nosing would move the contaminant 
out 1 in. from the back of the tread below, leaving the same amount of 
tread contaminated as was assumed. 
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PREFACE 


This appendix is intended to serve both as a source of immediately 
applicable methodology and as a guide to the underlying gas dynamic 
theory. Those interested more in applying the methodology than in der¬ 
ivations and comparisons of calculated and observed results will find the 
following parts of this appendix of particular importance: 

Location _ Description _ 

Section IIA The simplest and fastest method of estimating 

average pressure in a room as a function of 
time; adequate for many purposes. 

Section IIB3 Two different methods for a step-by-step, hand 
calculation providing average room pressure as 
well as dynamic pressure in the opening, valid 
for all flows through a single opening into a 
single room when the outside pressure is known 
as a function of time. Either method may be 
used but Method F is simpler for inflow, 

Method D for outflow; they may be combined in 
use. 

Section IIC Formulas for geometrical extent of the jet 

created inside the room by inflowing air and 
for dynamic pressure distribution within it. 

Table E-3 A computer program to calculate average pres¬ 

sure within a single room and dynamic pressure 
in as many as eight openings as functions of 
time when a single room has openings into sev¬ 
eral different pressure fields (e.g., a room 
with front, rear, and side windows struck by 
a blast on the front wall). 

Section IVB A step-by-step hand calculation of average 

room pressure and dynamic pressure in each 
opening during filling of a single room 
through one or two openings into separate 
pressure fields. 
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When using the methodologies noted above, the user may find the mean¬ 
ings of symbols under Notation at the end of this appendix. Subscripts 
not explained there refer to physical spaces, i.e., subscript 1 indicates 
quantity is measured outside the room; other odd-numbered subscripts refer 
to interior of rooms and even-numbered subscripts refer to connecting ducts 
or openings. 
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Appendix E 


ROOM FILLING FROM AIR BLAST 


I Introduction 


When a blast wave strikes a building, even should the structure with¬ 
stand the initial impact, the resulting inflow of air through windows and 
other openings can be critical in determining the safety of any people 
sheltered by the structure and in determining the response of the struc¬ 
ture itself to the blast impact. Although the physical laws obeyed by 
moving gases are well known and the course of the inflow in filling the 
building can in principle be calculated completely, any such calculation 
is far too lengthy to be practical for most purposes; fortunately, sim¬ 
plifications can be introduced which greatly shorten the labor of estimat¬ 
ing effects of the blast inside the building and which give results in 
good to fair agreement with experiments done with small models. In gen¬ 
eral , the effect of the inflow is to provide a stream of fast moving air 
in the shelter space which may (1) endanger shelterees by hurling them 
against large relatively fixed objects or by hurling objects against them, 
and (2) provide a back pressure on the inner surfaces of structure walls 
countering the blast pressures on their outer surfaces. 

Several factors enter into the calculation: the pressure outside 
each wall with openings and the time each opening becomes available, the 
area occupied by each opening and the volume of each room, the number of 
connected rooms and the area of each connection, and the ambient pressure 
and temperature in the building before the blast strikes. Perhaps the first 
of these to consider is what proportion of the wall exposed to the blast 
is open. If this fraction is greater than one half, the shock front lead¬ 
ing the blast wave will pass into the building only slightly weakened and 
subsequent inside pressure should be estimated from a knowledge of shock 
position and the laws of shock reflection. Methods appropriate to this 
case are only touched upon here. On the other hand, should the fraction 
be less than one tenth, clearly the filling is not a shock process and 
the methods treated here are quite pertinent. Unfortunately, in many 
applications the fraction of open area will lie between these two ex¬ 
tremes and, in these cases after the room filling calculation has been 
completed according to the methods suggested here, some thought must be, 
given independently to the influence of the entering shock front. 






When the source of the blast wave is an explosion and the location of 
the building in relation to the point of explosion is known or postulated, 
"free field" pressure histories at the building site can be found in stan¬ 
dard references, 1,2 ' 1 ^ < and from these histories well known methods 1 * 3 are 
available to derive approximate histories on the outside of the walls of 
the building. Briefly, these methods account for a short-lived peak of 
pressure created.by the impact of the front upon the wall nearest the ex¬ 
plosion, the relatively fast erosion of this high pressure to a level which 
is the sum of the free field pressure plus a drag pressure on the wall 
due to the high winds behind the blast front. This quasi-steady pressure 
then decays slowly to zero as the blast wave moves onward past the struc¬ 
ture . 


Ordinary window glass breaks rather quickly, i.e., within 8 ms (milli-i 
seconds) or less when struck by blast overpressure of 1 psi (pound per / 
■-square inch) or more. 4 Doors may withstand outside pressure longer, or 
even altogether. The time an opening becomes available with respect to 
the first impact of the blast upon the building becomes, then, the break¬ 
ing time plus the time required by the wave to travel from the wall near¬ 
est the explosion to the opening. If the opening is in the wall nearest 
the explosion, travel time is of course zero. Strictly, the decay of the 
blast wave overpressure which occurs during this time must be taken into 
account, but when the blast arises from a nuclear explosion of yield 
greater than a few kt (kilotons) this decay is slight and negligible; that 
is, a single free field pressure history for all openings may be assumed. 

It should be emphasized that the methods given here are simplified 
and their use can lead only to estimates. They are intended to provide: 

(1) calculations applicable to hand computation by those untrained in gas 
dynamics and (2) approximate results useful until more careful calcula¬ 
tions are made. Only in the case of the simplest structural configura¬ 
tions and the simplest pressure history shapes can limits of error be sug¬ 
gested for these results. Such cases will be the subject of the discus¬ 
sion immediately below. 


II Classical Nuclear Blast Wave Incident upon a Single Roo m 
A. Estimation of Inside Pressure Histo ry 

The "classical" blast wave from nuclear explosions consists of a steep 
pressure front or rise followed by a long-lasting decay phase during which 
the pressure in the wave falls to zero. It is accompanied by high winds 
giving rise to dynamic pressure against objects in the stream. Striking 
a wall at normal or near normal incidence, it creates a high pressure 


* References are listed at the end of this Appendix. 
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zone at the surface, which however is rapidly eroded as relief waves move 
across the wall face from the edges. As a first approximation in the cal¬ 
culation of room filling, this reflected phase can usually be neglected. 
Following the decay of the high reflected pressure, the quasi-steady pres¬ 
sure (free field plus dynamic) remains against the wall for hundreds of 
milliseconds to several seconds, depending on explosive yield. Generally 
filling is complete before this quasi-steady pressure has fallen more than 
a few percent; hence, as a first approximation the outside pressure may 
often be considered constant and, if only a single wall with opening is 
exposed to the blast, the time AT (in milliseconds) to complete filling 
may be computed as the ratio 35 , where V is room volume in cubic feet and 
A is area of opening in square feet.+ The average room pressure at any 
time t during the filling process is then simply the fraction of the quasi¬ 
steady outside pressure given by the ratio . For the purposes of this 
calculation, areas of several openings in the same wall should be added 
together to form the quantity A.* 

In case there are two or more walls with openings exposed to the blast 
and each such wall sustains a different outside pressure history (as will 
happen, for example, when the drag coefficient is different for two walls), 
the calculation is more complicated but first estimates of filling time 
and average inside pressure during filling may be found by adding interior 
pressures calculated as if each wall alone were exposed. As an example, 
consider a room of volume 30'xl0'xl0' = 3000 ft 3 in which the front wall 
has total openings of 36 ft 2 and side walls have total openings of 
60 ft 2 . The ratios for the front and side walls are 41.7 and 25 ft, 

respectively. If the quasi-steady overpressure on the front wall is 10 psig 
(pounds per square inch - gauge) and on the side wall 8 psig and if further 
the side wall opening becomes available 10 ms after the first blast im¬ 
pact , then the average inside pressure will be approximately as shown in 
Figure E-l by the heavy line OAFG. In other words the room will fill in 
approximately 24 ms. Lines OAC and DE represent filling rates through 
front and side walls,respectively; and ordinates of OAC and DE are added 
to form the line OAFG. Of course after the average inside pressure ex¬ 
ceeds 8 psi there will be outflow through the side wall; to allow for 
this loss the line FG has been placed between the outside pressure at 
the side wall (8 psig) and the outside pressure at the front wall (10 psig). 
The ordinate at FG is closer to 8 psig than to 10 psig because the area 
of the opening in the side wall is greater than that in the front wall. 

The line FG is intended to represent the final quasi-equilibrium pressure 
in the room. 


* The experimental justification of most of the procedures described in 
this section is demonstrated later in Figures E-4, E-5, E -6 and E- 8 . 
t Empirical relationship, dimensionally inconsistent. Meanings of 

symbols as used in this Appendix are defined as introduced and under 
"Notation" at end of Appendix. 
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TIME-ms 

FIGURE E-1 APPROXIMATE FILLING RATES THROUGH TWO WALLS 

These simple calculations do not apply when the reflected pressure 
lasts an appreciable length of time, or when the wave is non-classical 
such as it would be were a precursor present. Under these conditions the 
more detailed methods set forth below must be followed. 


B. General Case 

Justification for use of the simple methods noted above rests upon 
experience with a step-by-step calculation and comparison of its results 
with experiments. This calculation applies the principles of steady isen- 
tropic flow in ducts in successive, small time intervals. Conditions com¬ 
puted for the end of one time step become initial conditions for the next 
step. Conservation of energy, momentum and mass along with the assumption 
that the air behaves as a perfect gas with constant specific heats deter¬ 
mine the thermodynamic variables, pressure, temperature and density, as 
well as the wind speed through the opening. Unique expressions lending 
themselves to simple calculation cannot be given for the laws of conserva¬ 
tion of energy and momentum, and alternate forms leading to somewhat dif¬ 
ferent results will be stated. All such expressions rely upon certain 
approximations to the conservation laws and these approximations usually 
introduce errors into the results in comparison with which the approxima¬ 
tions arising in the assumptions of isentropy, perfect gas behavior and 
constant specific heats are negligible. 
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1. Inflow. Figure E-2 shows the idealized room with a single open¬ 
ing struck head-on by a blast wave. Three regions are noted: the outside 
© ; the doorway (2) which serves as a duct connecting the outside with the 
room ( 3 ). In order to make the calculations tractable, uniformity of con¬ 
ditions in each of the two regions © and ( 3 ) and over the cross section of 
region ( 3 ) is assumed; furthermore, during each small time interval At, steady 
conditions are assumed in each region. During the aforementioned quasi¬ 
steady state outside the building these assumptions are probably valid for 
region© but they clearly introduce error if the reflection or diffraction 
phase lasts an appreciable time, for during that episode relief waves are 
moving into the region from the edges of the building as well as from the 
doorway itself causing rapid fluctuations in wind speed and pressure. 

(Some account is taken of changes in pressure during the diffraction epi¬ 
sode by the standard techniques of estimating outside pressure.) Similar 
remarks can be made concerning regions © and ( 3 ), but if we are content to 
deal with "average" pressure and speed in those two regions, we may apply 
the step-by-step isentropic analysis. However, our present methods do 
not provide for any apportionment of gaseous energy in region © between 
streaming kinetic energy and internal energy; for simplicity of calcula¬ 
tion it will be treated as entirely internal at all times, which will 
cause overestimation of pressure and neglect of winds within the chamber. 
In evaluating the wind threat, the speed and dynamic pressure in the 
duct (2) must be regarded as the upper bounds on wind speed and dynamic 
pressure in region©. Later, methods will be given for estimating change 
in dynamic pressure as the wind moves into the room. 



x - AXIS 


FIGURE E-2 THE FILLING CHAMBER 
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In writing the conservation equations, two views can be taken of con¬ 
ditions in region ( 5 ). On the one hand pressure, density and wind speed may 
be those of the free field behind the blast front or, on the other hand, 
the air -upstream of the opening may be treated as stagnate at a pressure 
above free field either by the amount of the reflected pressure or by the 
amount of the product of the drag coefficient and dynamic pressure. Pro¬ 
vided drag coefficients are known, the second view is more simply applied, 
especially when the blast front does not meet the wall head-on. In what 
follows,the pressure, P x , and density, Pl , in region© will be those of 

stagnate air outside the wall. The work done in moving a mass element 

Am in region © through a small distance Ax toward region ©is: 

P 

P A dx = P AV = — • Am 
11 1 1 Pi 

where is the cross sectional area, and AV is the volume occupied by 
Am in region©. The mass element carries with it the internal energy it 
had in region©, i.e., 

, P 

1 . _1 • A®. 

y - 3. Pi 

where y is the ratio of specific heat at constant pressure to the specific 
heat at constant volume. (The perfect gas equation of state is assumed; 
see Ref. 5.) If the flow into the room is steady, energy conservation 
requires that the same total energy, specifically, the sum 



be given up within region © during the same time interval. Furthermore 
mass conservation asserts that the mass element moving through region© 
toward region© equal Am. The work done in region © in pressing the mass 
element toward region© is 

P 

2 Am 


and the internal energy in the element is 
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1 


• Am 



where the subscript Q) denotes conditions in region (§). Since the air is 
flowing into the room, however, the element in (2)also carries (streaming) 
kinetic energy of amount 


where u designates particle or material speed. Thus, if conditions are 
not changing too fast, we can write (cancelling out the factor Am): 



2 12 

— + — u 

Po 2 2 


( 1 ) 


To apply conservation of linear momentum we consider a control sur¬ 
face , shown dashed in Figure E-3. Neglecting gravitational forces, the 
x-component of the force integrated over this surface must equal the rate 
of flow of the x-component of momentum out of the volume plus the rate of 
increase of x-momentum within the volume. 6 Neglecting frictional and 
viscous effects, the only force on the surface is the thermodynamic pres- 



FIGURE E-3 FIRST CONTROL SURFACE 
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sure, exerted normal to the surface in an inward direction. Across the 
throat of the duct this pressure is taken to be P 2 ; everywhere else on 
the surface it is approximately P x . (Around the entrance to the duct 
pressure on the wall will be less than Pi, but the mass flow rate turns 
out not to be highly sensitive to corrections made for this effect and, 
in order to simplify calculation, uniform pressure P x will be assumed. 

The mass speed u 2 however would be reduced by this correction.) During 
the quasi-steady episode momentum within the surface will change only 
slowly and we assume that it is in fact constant during each successive 
time increment At. (This assumption is clearly false during the shock 
diffraction episode due to the presence of fast moving reflected shock 
fronts within the surface. During the quasi-steady episode the assump¬ 
tion neglects the relief wave spreading into the stagnate air outside the 
opening.) Thus, carrying out the integration over the surface, we find 


f P cos 0 da = (P - P ) A 
J 12 2 

In this integral 0 is the angle between the inward normal to the surface 
and the positive x-direction and da is an element of area on the control 
surface. The total momentum in the positive x-direction passing through 
the surface is simply 

2 


Hence, momentum conservation reduces to 


P 

1 



p„ u 
2 2 


( 2 ) 


Finally we note that even in the presence of moderately strong or weak 
shocks the isentropic equatiori of state of a perfect gas is accurate enough 
for this approximate calculation; hence 


P 2 


P 1 



1/Y 


(3) 


(See Ref. 7) . 

Given P x and p 1; Eqs. (1), (2), and (3) may be solved for P 2 , p 2 , 
and u 2 . The result can be written as: 


* For diatomic gases like air, y=1.4; see Ref. 8. 
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Y + 1 


/y p 


Y + 1 


which is independent of p , 


2 . Y - 1 

11 * - r - A - 


and B = —,Eq. (4) can be put in the form 
Y + 1 

1 

By ^ = y + A (5) 

When A and B have the values stated above, Eq. (5) has two solutions, 
one of which is y = 1 and the other is y = 0.1912. The second solution 
is the only one of interest here and will be designated y . 


To continue the calculation pj_ must be known. This value can be 
found from the Rankine-Hugoniot relations and knowledge of the strength 
and angle of incidence of the original shock front (Ref. 1,2), or it can, 
with enough accuracy for incident shock strengths less than 15 psi, be 
computed from ambient conditions using the isentropic equation of state, 


i .e. , 


P x “I !/Y 


where P and n are ambient pressure and density, respectively. With p 
O h o 

known, air density and pressure in the opening can be calculated from: 


p 2 - p x y Q Vi 


P = y P 
2 o 1 


from which wind speed becomes: 


P - P 
1 2 


P (1 - y ) 
1 y o 


P-. y 

r l o 


The mass flow into the room (3) can be written 
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^ = P 2 U 2 V* 


1 1 

= p. d - y ) y Y p 2 a At do) 

j_ 1 o o J 2 

If V 3 is the volume of the room and the prime is used to denote conditions 
in the room at the beginning of time step At, then average density in the 
room Pg at the end of At can be written as: 



(ID 


To find pressure in the room at the end of the time step we assume that 
all the energy lost in region (l) appears as an increase of internal en¬ 
ergy of the gas in the room, i.e., 


-1— • — * Am = — 

y - i p 1 y 




P - P 
3 3 


which can be solved to give P in terms of known quantities 


v P 

P Id . *2 + P < 

3 Pi v 3 3 


At any time air temperature T g within the room can be calculated from 
the perfect gas law: 



where R is the gas constant for air in the appropriate units (e.g„, in 
metric units R = 0.3028 joule/g°C). The quantity T 3 can reach high values 
as a result of the compression existing behind the shock and within the 
room; however, if airflow to the outside (outflow) is maintained, the 
relaxation of pressure following the passage of the front will return 
room gas temperatures to safe levels before injury to occupants is likely. 
Only the long lasting increase in room temperature resulting from fires 
will normally be a threat to the shelterees. 
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Conditions existing in the chamber at the beginning of a time step 
are used in the foregoing calculations only in Eqs. (11) and (12) and do 
not influence duct parameters because transients have been omitted from 
consideration. Transient phenomena, for example, determine the direction 
of flow; that is, if P-j^ > P 3> flow is inward as discussed above, but 
otherwise flow is outward. Repeated neglect of signals originating from 
the room leads to an accumulation of error in the calculation of average 
pressure as can be seen from comparisons between calculation as above and 
measurement shown in Figures E-4, E-5, and E-6. The experiments were 
carried out in a 24-inch shock tube; the configuration of each model 
chamber is shown in an inset in the figure; and the foregoing calcula¬ 
tion produces Curve F of the figures. The curves in Figures E-5 and 
E-6 labelled "external history (A)" are measurements in the free stream 
by means of a pitot tube oriented with respect to the stream to conform 
with the orientation of the opening in the model room;"external history 

9 

in Figure E-4 is side-on overpressure in the unobstructed shock tube. 

In each case the calculation initially yields pressure in agreement with 
observation but eventually shows room pressure during filling in excess 
of measurement, although the maximum difference is 20% or less. Magni¬ 
tude of At for these calculations was one-quarter the transit time of a 
sound signal across the room. 
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OVERPRESSURE (PSI) 


14 



FIGURE E-4 FILL HISTORIES FOR SIDE-ON INCIDENCE 
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OVERPRESSURE ( PSI ) 



TIME ( MILLISECONDS) 


FIGURE E-5 FILL HISTORIES FOR A FACE-ON MODEL 
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OVERPRESSURE (PSI) 



TIME (MILLISECONDS) 


FIGURE E-6 FILL HISTORIES FOR A REAR FILL MODEL 
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Choice of the size of At is somewhat arbitrary except that (1) value 
much greater than sound transit time will give a false idea of the degree 
of irregularity in the fill process and (2) enough steps should be taken 
to make it possible for the influence of variations in P to be shown in 
the results. The length of the bar labelled "t" in Figures E-4, E-5, and 
E-6 represents the sound transit time across the longest room dimension. 

The existence of significant theoretical errors in our treatment of 
flow into a room by quasi-steady analysis is clearly revealed by consider 
ing the single control surface formed by superposition of that shown in 
Figure E-3 and that indicated with dashed lines in Figure E-7. Such a 
surface coincides with the inner surfaces of the room and passage and 
extends into quiescent air outside. Under our hypotheses there is no 
flow through this surface anywhere and no change of momentum within it, 
yet the surface integral of the x-component of pressure over the bound¬ 
aries does not vanish. 



FIGURE E-7 SECOND CONTROL SURFACE 
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This absurdity can be avoided in one or both of two ways. A term may 
be added to the right side of Eq. (2) to account for the changing flow 
pattern within the control surface shown in Figure E-3 or a term may be 
added to the left side of Eq. (2) to account for the possible nonuniform 
ity pressure over the boundaries of the surface. As the room begins 
to fill, a rarefaction wave moves back into the high pressure gas outsid 
the first doorway bringing more and more gas into motion toward the open 
ing. In other words, the neglect of the rate of change with time of the 
momentum within the control surface outside the room may be at least one 
cause of the contradiction noted above. Any attempt to calculate a cor¬ 
rection for this effect would certainly add to the complexity of these 
simplified procedures; furthermore, the degree of agreement between ob¬ 
servations and theory of Method F (shown in Figures E-4, E-5, and E-6) 
suggests that the added effort to account for the rarefaction wave may 
not be needed to achieve the desired degree of accuracy. 

Some writers 10 ' 11 > 12 » 17 use the equation 


instead of Eq. (2). Justification for Eq. (13) is based on analogy with 
the treatment of flow into a large chamber steadily being evacuated. 13 
Equation (13) of course provides continual coupling between flow condi¬ 
tions and conditions in the room. 


Using Eq. (13) we derive the pressure buildup inside the room in 
the following way. Substituting Eq. (3) into Eq. (1), then replacing 
P 2 with according to Eq. (13) and solving the resulting equation for 
u 2 , we find: 



from which we calculate the mass inflow in the time increment At: 



Whenever Eq. (13) is employed, empirical corrections must be made to re¬ 
duce the calculated inflow rate or the calculation will not yield realis 
tic values for room pressure in small experimental models; the simplest 
correction is the discharge coefficient, 13 represented by the factor K 
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in Eq. (14). Investigators at the IIT Research Institute 17 have found 
it necessary to use the value K=0.7 to reconcile computed pressure rises 
in small models with those measured; Curves D in Figures E-4 and E-5 have 
been produced by a calculation based on Eq. (13), with K=0.7.during 
inflow and K =1.00 during outflow. 


The value of the discharge coefficient is usually discussed in con- 

13 

nection with boundary layer thickness and the Reynolds number. It 
should be noted therefore that the relatively good agreement between 
the observed room pressures and Curves D was obtained in very small 
models and not in full-sized rooms. To provide an estimate of the in¬ 
fluence of the value of K on calculated pressure rise, Curve G, based 
on the value K = 1.00 during both outflow and inflow, has been entered 
in Figure E-4. In the flow into full-sized rooms, presumably, the Rey¬ 
nolds number will be larger and the discharge coefficient more nearly 
equal to 1.00 than in the flow into small models. 


Finally, the pressure increment during the interval At is found by 
substitution of Eq. (14) into Eq. (12): 


P - P' = K y P 

O «J -i- 


p 

I /p.vl 

r /p'x 1 -^! 

i -i — ) 

1 

2 A At 

2y _1 

2 / 3 \Y 

2 

Y" 1 " P x 

VV 

[ \ 1/ 

CO 

> 


(15) 


Melichar 11 ’ 12 omits both the factors K and y before the right-hand 
side of Eq. (15), which is equivalent numerically to making K = 0.7 and 
y =1.4. He attempts to justify this procedure on theoretical grounds 
unconnected with boundary layer theory. 12 Some of his numerical results 
are shown as Curves C in Figures E-4, E-5, E-6, and E-8. Melichar em¬ 
ploys a value of At equal to the transit time of sound across the room. 


Pursuing the analogy between room filling and steady flow into a 
chamber held at constant pressure, we expect to encounter the phenomenon 
of choking as the ratio of room pressure to outside pressure drops below 
the critical value. 13 For a given set of reservoir conditions, i.e., for 
each pair of values P-l and p-^, the critical pressure ratio is that for 
which isentropic flow into the room achieves the maximum mass rate and 
for which the flow speed equals local sound speed. Therefore, to find 
this critical ratio, we can differentiate Eq. (14) with respect to P^/P-^, 

set the result equal to zero and solve for (Pn/P, ) ,^ , which yields; 

J 1 crit. 



E-27 


For every value of P 3 below the critical value as given above, mass flow 
into the room will be that obtained by substituting the critical ratio 
into Eq. (14), i.e.: 



Because the mass flow rate is limited in this way independently of the 
value of P 3 it is called "choked” flow. Flow obeying Eq. (14) is called 
"unchoked" or "subsonic." Numerically, when y = 1.4 the critical ratio 
equals 0.5283; thus, assuming ambient pressure is 14.7 psia,* the critical 
outside pressure is 


14.7 + 

~ 5283 = 27.83 psia or 13.13 psig T 

Peak overpressure in none of the experiments reported in Figures E-4, 

E-5, and E -6 rose above 13.13 psig; hence, according to the foregoing 
theory, choked flow should not have occurred. 

A similar degree of comparison between calculation and measurement 
is found in the results shown in Figure E -8 stemming from a 27 ft 3 model 
exposed to a large chemical explosion, except that the transient fluctu¬ 
ations associated with the entering shock front are more easily dis¬ 
cerned in the larger model than in the small shock tube models. The 
parameter values used in the calculations summarized in Figure E -8 were: 

P = 13.58 psia ambient pressure 

3 

P 0 = 0.0672 lb/ft ambient density 

Y = 1*4 ratio of specific heats 

Ag = 0.821 ft area of opening 

3 . 

Vg = 27.0 ft volume of model room 

The fraction of the impacted wall area occupied by the opening is slightly 
less than one-tenth. 

From the results shown in Figures E-4, E-5, E- 6 , and E- 8 , some esti¬ 
mate can be made of the validity of the greatly simplified method of com¬ 
puting pressure rise in a filling room set forth in Section IIA. In each 
figure, the value of V/2A in feet has been entered. A constant pressure 


* Pounds per square inch, absolute, 
t Pounds per square inch, gauge. 
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FIGURE E-8 COMPARISON OF THE COMPUTED VERSUS THE EXPERIMENTAL FILL HISTORY 
OF THE MODELED CHAMBER TESTED IN PROJECT DISTANT PLANE 




rise from zero time and zero overpressure to a pressure equal to outside 
pressure at a time equal to V/2A msec overestimates the fill pressure and 
underestimates the fill time in the small models but seems to give results 
in good agreement with both calculation and measurement in the three-foot 
cube model reported in Figure E-8, although the presence of important os¬ 
cillations (caused by shock waves) in the pressure record in the large 
model makes clear assessment difficult. It should be noted also that the 
fraction of the wall struck by the blast that is occupied by the opening 
(opening fraction) is over 9 percent in the large model while only 1.2 
percent in the small four-in. cube models. The straight line estimate 
explained in Section IIA would appear to be good also in Figure E-4 for 
which model size is intermediate and opening fraction is almost 6 percent. 

The data from the models suggest that the simplified method of esti¬ 
mating pressure rise in a room set forth in Section IIA is adequate when 
outside pressure decay is slow as in a free-field nuclear blast wave. As 
will be noted later, when there exists an important diffraction phase in 
the blast wave interaction with the structure, more sophisticated methods 
may be justified. 

2 - Outflow . All the experimental records discussed so far show that 
room pressure eventually exceeds outside pressure by a small amount. If 
all openings are into the same outside pressure field and if the fall of 
outside pressure is steady and slow, this overshoot is not likely to be 
of practical importance. However, as will be seen later, if a room has 
openings that are affected by different outside pressure histories, a 
significant outward pressure may develop on one or more walls of the room. 
To calculate this pressure, the rate of outflow of air through each open¬ 
ing must be found. 


To compute inflow we treated the exterior atmosphere and the aperture 
leading into the room as parts of a system of ducts through which quasi- 
steady flow was maintained. Pressure increase within the room was com¬ 
puted as resulting from the transfer of mass and energy from the large 
outside reservoir into the room. Outflow may be treated in the same way, 
except that the direction of flow is opposite. For outflow, Eqs. (1) and 
(2) still apply, but the air in the duct now originates in the room and 
the adiabatic expansion law, Eq. (3), must be replaced by 


P 2 P 3 


P' 

L 3 J 


Y 


(16) 


Equations (1), (2), and (16) may be combined to form an equation identical 
to Eq. (5) except that now 

P 

2v ‘".SI 

B = 


2y 


y + 1 


p'. 
3 J 


Y 
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But, as before, 


y 


and 


A 


Y ~ 1 

Y + 1 


P 

P 


2 

1 


In this form, Equation (5) has two solutions in the range 1 > y > 0 
whenever y = 1.4 and 




> 1 


or, 



These two solutions merge into a single solution at y = 



(17) 


if 



(0.9094) 


P' 

3 



Furthermore, if y = 1.4, Equation (5) has two solutions but only one in 
the range of 0 ^ y <; 1 whenever 


Whenever 



P 

2: —^ > (0.9094) 
P 1 Y 



< (0.9094) 



(18) 


(19) 


no solution to Equation (5) exists. 

P 

Since the quantity ~y is constant along an isentrope and increases 
with entropy, when Inequality (17) is true during outflow, specific 
entropy outside is greater than or equal to specific entropy within the 
room. 


The specific entropy (i.e., entropy per unit mass) in the room at 
any time during filling can be formally calculated as follows. 

First we note that temperature of air in the room rises during the 
whole inflow period, as can be seen by treating Am as a true differential 
and writing Equation (12) as: 
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• dm 


dP = - 

3 P 1 v„ 


From the perfect gas law 


P V = m R T 
3 3 3 3 


(where m 3 is the mass of air in the room, and R the gas constant for air) 
we write 


R T m R 

3 3 

dP o = "T— • dm + - • dT 

3 V 3 V 3 3 


P 3 m 3 E Y P, 

dP = — * dm + - • dT = - 

3 m 3 v 3 3 P x V 5 


v P P 
1 2 _1 _3 

V 3 P 1 ~ p 3 


m R 

dm = • dT 

V 3 


[^1 “ T 3 ] dm = m 


dT 

3 3 


Clearly, at the start of filling T x > T 3 , hence at the start dT >0. 
Moreover Equation (20) shows that T 3 can increase above T-^ while dm > 0 
until T 3 approaches y T^. Thus as temperature outside (T^) falls due 
to the adiabatic relaxation behind the shock front, inflow and rising 
inside temperature continue. 

In an isentropic process involving a certain mass of ideal gas the 
quantity _P_ is constant; or, taking differentials, 

T V(y - 1) 

dp _ dT 
p ~Y t - T 

Rewriting Equation (20) as 


VT, - T 
3 Y 1 3 
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and dividing numerator and denominator by and noting that density 
p 3 = , we find 


v T - T 
Y 1 3 


In other words the process of room filling during the time period t to 
t + dt results in a temperature increase dT^ that bears the relation 


(dT ) . 

3 isentropic 


v T - T 

T 1 3 

v T - T 

Y 3 3 


to the temperature increase in an isentropic process providing the same 

■rl'T' 


density increase. When T-^ > T 3 , 


(dT3) igen t ro pi C 


Now the entropy change dS in a process at constant density is 


c dT 
v 


where c y is specific heat at constant volume. If we imagine a return to 
the isentrope (after an increment of filling) by a process at constant 
density or volume, the specific entropy change within the room resulting 
from the filling increment becomes 


'i [' 


dS = — dT - (dT ) . 

3 T 3 3 isentropic 


c y 
v • 

dS = —— 
3 T 3 


T - T 
1 3 

Y T - T 
_ Y 1 3 J 


which on substitution of the quantity dTg from Eq. (20) becomes 


C v Y 

ds - — (T - T )dm 

3 T m 1 3 

3 3 


Therefore, during inflow (i.e., when dm > 0) specific entropy within the 
room will increase until T 3 = or P 3 = whichever occurs first. 
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Because of the passage of the shock front, air outside has greater 
than ambient specific entropy. Hence, initially before inflow begins 


but Eq. (22) shows that the inflow process, according to the theory ap¬ 
plied here, creates entropy within the room. Should enough entropy be 
so created, Inequality (19) may eventually be satisfied and solution of 
the outflow equations presented here may be impossible. The rise in spe¬ 
cific entropy within the room as a result of inflow is: 


As = c v 
3 v Y 


P =P 

/ 3 1 T - T 

_1_3 dm 

T m dt 
-.0 3 3 


where t = time and , T 3 , and m 3 are functions of time. From Eq. (20) 
we know that T 3 rises asymptotically toward yT]^ until P 1= P 3 ; thus T 3 
may be larger than T-^ For a given function P-^t) [which determines 
T l( t )» T 3 (t), and m 3 (t)], the maximum specific entropy within the room 
will be reached when 

T 3 = T 1 

If we assume that at this time inflow is still under way, then 

P > P 
1 3 

which implies,when y > 1, that 

1- - 1- i 

P Y >P Y 

1 3 


I_1 I_! 

y 1 Y 1 

P < P 

1 3 


Now from the equality of the temperatures at this time, the perfect gas 
law implies 


P P 
_1_3 

P 1 = P 3 


Multiplying both sides by the last inequality above, we find 
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or 




P 


3 



P 


3 


P 


Y 

3 


Hence, in general, the possibility may not be ruled out that Inequality 
(19) will be satisfied. Whether the outflow equations as presented here 
have a solution will depend on the nature of the function P^(t). If out 1 
flow is to be calculated, then to provide an initial disparity between 
specific entropy outside and inside, the inflow period must be treated 
using the initial density behind the shock front and computed from the 
Hugoniot relation: 1 

(Y+ 1)P io + (Y-D p 0 

P lo = P o (v-l)P + (Y+l) p ( 23 ) 

lo o 


Subsequent 

law: 


outside air densities may be calculated from the adiabatic 


P 


1 


lo 



(24) 


In Eqs. (23) and (24), p lo and P lQ refer to the air density and ab¬ 
solute pressure immediately behind the shock front. When peak pressure 
at an opening is reached by reflection of an incident shock wave from a 
wall, Eq. (24) is not correct regardless of whether P lo is taken as peak 
incident absolute pressure or peak reflected absolute pressure, but the 
error in using Eq. (24) is small. We will arbitrarily consider p lo and 
P-^ o as representing conditions behind the free—field shock front. Use 
of Eqs. (23) and (24) may make it possible to satisfy the reverse of 
Inequality (19). 


Should such still not be possible, other calculational methods must 

11 12 

be used for the outflow phase, such as that proposed by Melichar ’ or 
that reported by IIT Research Institute . 17 These methods equate duct 
pressure P 2 with outside pressure during outflow; the IIT investigators 
then fit observed outflow pressure data by choosing values for the dis¬ 
charge coefficient and for the ratio of inside to outside pressure at 
the time of flow reversal. 
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3. Outline of Hand Calculation . In constructing from the foregoing 
equations a calculational scheme for estimating the parameters of flow 
into and out of a single room with a single opening we start with a series 
of values of P-^, one for each time step. These may be obtained by linear 
interpolation from a given table of outside pressure as a function of time 
and each value should pertain to the center of the time interval. The 
size of the time interval, At, itself is arbitrary, but it should be no 
greater than the quantity T; presumably up to a limit, greater accuracy 
results from smaller values of At. The size of At may be changed during 
the calculation when the rate of change of flow parameters changes. We 
also need values for ambient pressure P Q and density p Q . 

Two methods of calculation are shown below. The first is that used 
to produce Curves F in Figures E-4, E-5, E-6, and E-8, namely, that based 
on Eqs. (1), (2) and (3) for inflow or Eqs. (1), (2) and (24) for outflow, 
and in the outline below it is called Method F. This method has the ad¬ 
vantage of great simplicity and of not requiring knowledge of empirical 
constants; however, as will be explained later, values of wind speed and 
dynamic pressure computed by it are subject to doubt in some cases and 
for that reason a method given by IIT Research Institute is included also. 
The latter method is responsible for Curves D in the figures and, there¬ 
fore, in the outline below it is called Method D. As noted earlier, 

Method D is numerically equivalent to the calculation used by Melichar. 11 » 12 

Average pressure inside a room and dynamic pressure in the single 
opening to an outside reservoir whose pressure variation in time is known 
may be calculated as functions of time by the sequential application of 
the steps stated below. Each cycle through a series of steps completes 

the calculation for one time interval. The first three steps are executed 

only during the first cycle; subsequent passes begin with step (4), as in¬ 
dicated in the outline. Step (5) is a branch point to separate sequences 
for inflow and outflow, chosen according to a criterion given in step (5). 

There are further branches (a) to Method D or F, chosen at the discretion 

of the user at each time, and (b) under Method D to choked or unchoked 
flow, determined by stated criteria. Throughout the outline, the quantity 
Y has been set equal to 1.4. 

(1) Set P' = P and p' = p and t = 0. 

o o 3 o 

(2) Compute 6P + P 

lo o 

p = p - 

lo o P + 6P 
lo o 

where P is the absolute pressure immediately behind the 
shock front. 
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(3) Choose value of At (see opening paragraph of this Section IIB3). 

(4) Determine outside pressure at the current time; i.e., deter¬ 
mine P^, from the known reservoir pressure history (pressure 
variation with time). During first time interval p 1 = p lo ' 

(5) Determine direction of flow; i.e., if > P 3 , flow is inward; 
go to step ( 6 ). Otherwise flow is outward; go to step (28). 


Inflow 


(6) Compute 


P 1 = P lo 


lo 


0.7143 


Branch to selected Method below for step number (7D) (Method D) 
or step (7F) (Method F). 


Method D (Inflow) 

(7D) If P'/P ^ 0.5283 inflow is choked; go to step ( 8 D). Otherwise 
3 1 

inflow is unchoked; go to step (18D). 


Choked Inflow 


( 8 D) Am = K 

choked 


2,4 

2 \ 0.4 * 


/ 2 \u.4 

1 ‘ 4P l P l(i^) F 


1/2 


A At 
2 


, 1/2 


= 0.6847 K [p P F] A At 
11 2 


Using the recommended value 17 of K = 0.70 this becomes 

Am = 0.4793 [ pPF] 1 / 2 A 0 At 

choked 11 2 


(9D) P 2=P ; 

(1QD) p 2 = P x 


,0.7143 



* The factor F will often be necessary for consistency of units. For 

example, if u 2 is to be in ft/sec and P-^ is in lb/in 2 and p^ in lb/ft , 
then F = 32.X 144. = 4,608. 
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Am 
T 3 

If desired, dynamic pressure at the opening may be found 
from: 

2 

P u 

12 2 

< 13D > -2=i — 

< 14D > P 3 * >3 
< 15D > P' = P 3 

(16D) Advance time by amount At and return to step (4). 

Unchoked Inflow 


. . choked 

(11D) u = -- 

' 2 At A p K 

2 2 


(12D) P = P; + 1.4 - 
3 3 p. 


(18D) 

(19D) 


P = 


P' 

3 


Pn = P- 


0.7143 


(20D) u o = 7 


P 1 P 2 


(21D) Am = p„u ‘A At 

unchoked 2 2 2 


P 1 ^unchoked 

(22D) P = P' + 1.4 — - 

3 3 p. Vo 


(23D) 

(24D) 


Am 


P = p' + 
3 3 


unchoked 


If desired, dynamic pressure, q 2 , in the opening may be 
found from: 


1 u 2 
^2 _ 2 ^2 F 
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(25D) 

P 3 = P 3 



(26D) 

P 3 = P 3 



(27D) 

Advance 

time by amount At and return to step 

(4). 

Method 

F (Inflow) 


(7F) 

P = 0.1912 P 

2 1 


(8F) 

P 2 P 1 

-| 0 * 7143 

— = 0.3067 p 

L p J 1 


(9F) 

2 

U 2 = 7 

r p i p 2 1 r p ii 

-F = 2.637 — F 

_ P 1 p 2 J L P i- 


(10F) 

S =U 2 

1/2 

P 2 A 2 At = 0.498 (P 1 P 1 F) A g At 


(HF ) 

P 3 = P 3 

P Am 

, . 1 3 

+ 1.4 

P v„ 

1 3 




Am 

*3 


(12F) 

P 3 “ P 3 

J 

V 3 


(13F) 

If desired, the dynamic pressure, q^, in the 
be calculated: 

doorway may 

(14F) 

q 2 2 P 

2 

u 2 

1 - = 0.4044 P 

2 F 1 


(15F) 

3 M 3 



(16F) 

P' = P„ 

3 3 



(17F) 

Advance 

time by amount At and return to step 

(4). 


Outflow 

(28) Branch to selected Method below for step (29D) (Method D) 
or step (29F) (Method F). 
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Method D (Outflow) 


(29D) If P^/Pg ^ 0.5283, outflow is choked; go to step (30D), 
Otherwise outflow is unchoked; go to step (39D). 


Choked Outflow 


(30D) Am = -0.6847 K [p'P'F] 1/2 A At 

choked 33 2 


For K = 1.0 this becomes' 1 


1/2 

Am = -0.6847 [p' P' f] A At 

choked 33 2 


(31D) 


(3 2D) 


P = P 
2 1 


P = p — 
2 3 P 

| 1 


0.7143 


. choked 

(33D) u = -- 

2 At A K 
2 v 2 


(34D) 


P 3 = P 3 + X * 4 


Am , . P^ 

choked 3 


(35D) p = p' + — 
3 3 V 

3 


(36D) 


P' = P 
3 3 


(37D) 


p' = p 
3 M 3 


(38D) Advance time by amount At and return to step (4) 


Unchoked Outflow 


(39D) 


(40D) 


P 3 P' 


P 10.7143 
2 


(41D) 


r p pi 

_ 1 2 

7 —-F 

L P 1 P 2 
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The theoretical differences between the two methods outlined above 
are striking. At the present time there is not enough experimental evi¬ 
dence from full-sized rooms to permit a clear judgment between them. 

If the room fills slowly, that is, the filling time is long compared 
with the transit time of sound across the room, then the analogy under¬ 
lying Method D should be good; but even so there is doubt about the cor¬ 
rect value of the discharge coefficient, K. The possible influence of K 
on the calculated room pressure is large, as can be seen from the differ¬ 
ence between Curves D and G in Figure E-4. If the filling is rapid, per¬ 
haps Method F should be preferred. However, as will be seen in the fol¬ 
lowing paragraphs, this Method can lead to anomalous values of wind speed 
through the opening. 


C. Wind Speed and Dynamic Pressure 

Other flow parameters of interest are wind speed and dynamic pres¬ 
sure. The present simple methods yield values of these parameters only 
in the doorway, but, ordinarily, values at the doorway may serve as 
upper bounds for the whole room. 

Unfortunately, measurements of wind speed and dynamic pressures in 
shock-filling rooms are even fewer than observations of room pressure. 
The differences among the predictions of speed among the several calcu- 
lational methods are large, but the values of dynamic pressure are often 
in fair agreement. For estimates of the acceleration of objects in the 
stream, the dynamic pressure is the only pertinent parameter. 

During both inflow and outflow wind speed may be calculated from 
Eq. (9) and, if the isentropic assumption is made, Eq. (6). Dynamic 
pressure, q 2 , is customarily defined as 
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Using Eq. (9) this can be written: 



Because of the neglect of transient effects, Eqs, (9) and (25) do not re¬ 
flect the influence of changes in the room pressure, and values calculated 
from them may be unrealistic. Peak speeds found from Eq. (9) for the ex¬ 
periments reported in Figures E-4, E-5, E-6, and E-8 are supersonic, i.e., 
above local sound speed in the doorway. 

Because of its use of the choking analogy Method D never calculates 
supersonic speeds; but it does neglect transient shock fronts which could 
alter the flow substantially in fast filling rooms. 

Table E-l contains the several calculated values of peak speed and 
dynamic pressure corresponding to each of the experiments reported in 
Figures E-4, E-5, E-6, and E-8. The table suggests that use of Method F 
may result in overestimate of that important criterion for damage potential, 
dynamic pressure. As noted above, because signals from the filling room 
into the outside reservoir are neglected, Method F in addition to over¬ 
estimating peak dynamic pressure may overestimate the time-average dynamic 
pressure in the doorway. 

There are no direct observations of wind speed in the opening to com¬ 
pare with calculations. Coulter, 9 however, had reported (1) measurements 
of the acceleration of an ^ inch diameter nylon ball placed in the doorway 
of a small model struck head-on by a weak shock wave, and (2) measurements 
of the motion of smoke columns inside these models during shock filling. 

With the help of acceleration coefficients derived (theoretically and ex¬ 
perimentally) by Bowen and others 14 we can compute motion of the nylon 
sphere from Eq. (25). Coulter allowed a shock front of overpressure equal 

to 4.89 psig to strike a reflecting plate in which an entrance 1x4 inches 
was cut. Since the chamber behind the opening had a volume to opening area 
ratio (V/A) equal to 1.33, filling was essentially complete in 1.33/2 = 
0.667 ms, during which time there was little or no decay of the incident 
wave. If we assume ambient pressure equal to 14.7 psi, reflection at 
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Table E-l 


PEAK VALUES OF OTHER FLOW PARAMETERS 


, \ v Method 


C* 

F 


G 

Experiment 

q 2 (psi) 

“ 2 ( f Ps) 

S 

U 2 

q 2 U 2 

Fig. E-4 

7.98 

1045. 

9.93 

1650. 

8.02 985. 

Fig. E-5 

10.2 

1180. 

11.2 

1680. 

10.2 1110. 

Fig. E-q 

7.79 

972. 

9.74 

1640. 

7.80 972. 

Fig. E-8 

4.08 

747. 

7.37 

1634. 



* Correct adiabatic law has been substituted for the inverted form 
appearing in Melichar's FORAST code. 11 


normal incidence of a 4.89 psi shock produces a reflected overpressure 
equal to 11.1 psi. 1 Hence the value of P^^ in Eq, (25) becomes 11.1 + 14.7 
25.8 psi, and 

1 

q = “ x 25.8 (1 - 0.1912) 

^ Ct 

= 10.5 psi 

2 

or q = 48100 lb/sec ft 


In this expression "lb" is a unit of mass. 

A still object in a stream of moving air is accelerated according to 
the formula 


dv 

= a q ( 26 ) 

where & is an acceleration coefficient characteristic of the object and 
q is the dynamic pressure as defined here. Eq. (26) continues to hold 
for the moving object as long as its speed is small compared to air speed. 
Bowen and others 14 suggest a value of a = 0.138 ft 2 /lb for a 1/8 inch 
diameter steel sphere. Since a is inversely proportional to mass, the 
corresponding value for a nylon sphere of the same size is approximately 
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7.87 x 0.138 = 1.09 ft 2 /lb (specific gravity of iron equals 7.87). Hence 
acceleration of the nylon sphere should be, according to Eq. (26), 

d v 2 

— = 1,09 X 48100 = 52200 ft/sec 
dt 

and the final speed of the sphere when filling is complete 667 psec after 
impact becomes: 


™ 6 

v = 52200 X 667 x 10~ = 34.8 ft/sec 

1 2 -12 

and the displacement is s = X 52200 X (667) X 10 = 0.0116 ft = 

0.14 inch. Coulter measured terminal speed after 700 psec as 29.8 ft/sec. 
Observed displacement in this same time was approximately 0.10 inch, which 
means the ball did not move far from the center of the doorway during the 
period of observation and was subjected to only slightly attenuated peak 
dynamic pressure. 

Corresponding to Eq. (25), Method D produces the following equation 
for the calculation of dynamic pressure in the opening: 



The value of q in Eq. (27) is not constant in time even when P^ does not 
change appreciably. For simplicity, we may assume P 3 / varies linearly from 
14.7 psia to P-^ over the filling time 0.667 ms. Again, we take P-^ = 25.8 
psia and find the following variation of dynamic pressure with time: 


time (ms) 

0 

0.2 

0.4 

0.6 

0.667 

P 3 ' (psia) 

14.7 

18.0 

21.35 

24.7 

25.8 

q„ (psi) 

8.97 

6.80 

4.15 

1.08 

0.00 


The tabulation shows a nearly linear variation of qIntegrating Eq. (26) 
approximately, we find the final value of speed to be: 
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v = 1.09 


lb 


ft 



8.97 —g 0.667 X 10' 
in^ 


ft 


sec 


= 15.0 


ft 

sec 


which is about half the observed value. The relatively better agreement 
with experiment furnished by Method F than by Method D may reflect the 
fact that the shock filling in this case is "fast." 

However, it would seem that any relatively good agreement between 
calculated and observed acceleration in this example was fortuitous, since 
Coulter's photographs of events in the chamber make it clear that the 
diffraction episode very nearly coincided with the calculated filling 
time, that is, air flow during the period of observation was more a shock 
related process than a case of quasi-steady adiabatic flow. Some support 
for this view comes from the second set of wind observations reported by 
Coulter: motion of smoke in the same kind of chamber used for measure¬ 

ments of the motion of the nylon ball. Most of Coulter’s observed values 
of wind speed and dynamic pressure from smoke motion are one or two 
orders of magnitude lower than comparable values calculated by the methods 
presented here. However, Coulter could not by this smoke method measure 
speeds in the doorway and, in fact, as filling went on his observations 
were restricted to regions farther and farther from the core of the flow. 

Vertex Corporation 10 and Melichar 11 ’ 12 have applied the theory of 
steady subsonic jets to estimate the attenuation and spreading of the 
incoming stream of high speed air. Melichar describes the early period 
after shock impact as including jet evolution but his method of estimating 
duration of this preliminary stage requires consideration of shock front 
patterns. Presumably flow speeds and dynamic pressure during evolution 
are less than those quantities in the fully established jet, so use of 
quasi-steady state conditions' is conservative for our purposes. According 
to Melichar, the steady jet consists of four regions, illustrated in 
Figure E-9 (reproduced from Ref. 12). In the conical core whose base is 
the opening into the chamber, horizontal speed and density uniformly equal 
the speed and density in the opening. The tip of this core extends into 
the room between 3 and 6 times the opening diameter. Vertex Corporation 10 
assumes that the flux of horizontal momentum is constant from one vertical 
section through the jet to another and gives the equivalent of the fol¬ 
lowing formula for the radius of the jet: 
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t TRANSITION 

. REGION OF SIMILIAR 

THE CORE EXISTS 

REGION 

VELOCITY PROFILES 



Source: p.33, Ref.12 


FIGURE E-9 SCHEMATIC OF JET FLOW 
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R = 0.178 


(28) 


L& + 5.6 (A/tt) 1/2 ] 

where SL = distance along jet axis from (circular) opening of area A. 

Vertex also suggests that doors may be considered circular openings for 
the purposes of this estimate. Now, if the flux of horizontal momentum 
is constant at all cross sections through the jet, and if 


I 


r»R 2 

p u 2 tt rdr = 

o 



q 2 tt rdr 


(29) 


(where r is radial distance from jet axis, p density, u horizontal par¬ 
ticle speed) then I is independent of SL . Melichar provides speed pro¬ 
files, i.e., u as a function of r and SL , for one value of pressure dif¬ 
ferential across the orifice and suggests that p be taken constant through¬ 
out the jet. Vertex approximates the function u = u(r) as: 



(30) 


where u* is the speed on the jet axis. Vertex further assumes 


u* \SL + 5.6 (A/tt)' 


constant 


which is only very approximately true according to curves reproduced by 
Melichar, but which allows one to calculate u in terms of speed in the 
opening, since at the opening SL = 0. Vertex does not assume p independent 
of r and i but provides formulas for p as a function of T and for T as a 
function of r and SL . 


For a first estimate of dynamic pressure in the jet, an equation for 
average dynamic pressure q across the jet will be derived from Eq. (29). 

If q 2 is the dynamic pressure in the opening then by equating the 
value of I at the doorway to its value downstream we can write: 
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2 

q A = q TT R 
H 2 



when Eq. (31) has been applied and presumably "safe" areas in a shelter 
found, then more detailed consideration must be given to the possible en¬ 
hancement in those areas of the value of q along the axis of the jet. It 
is suggested that values of u* be found from Melichar's curves 3 - 2 and air 
density in the jet be taken equal to the average in the room at the time 
considered. 

The theory of steady jets, upon which the foregoing considerations 
are based, stems ultimately from observations of continuous flow through 
small orifices in large reservoirs kept at constant pressure into an es¬ 
sentially infinite reservoir, also kept at constant pressure. In the pres¬ 
ent case flow is into a relatively confined space. Indefinitely continued 
expansion and attenuation of the jet cannot be expected. In fact there 
is a strong possibility that in some situations even the expansion rela¬ 
tively near the doorway itself will be reduced in amount. Circulatory 
flow patterns with relatively little turbulent mixing could be set up in 
small rooms with large openings. This is a subject needing further study 
by civil defense investigators. 

Ill Multiple Rooms 


Although the algebra becomes increasingly cumbersome as rooms are 
added, the foregoing principles of calculation can be applied to a series 
of connected rooms. For example, two rooms are represented in Figure E-ll 
for which two control surfaces may be employed. Through the smaller the 
momentum flux is due to flow through the first doorway and the resulting 
momentum equation may be written as: 


A (P 
2 V 1 


p ) 
2 


2 

p u A 
h 2 2 2 


(32) 
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Using the larger control surface we find: 



- P 


A + (p' 
2 v 3 


P ) A„ 
4 4 


2 

u 

4 


A 

4 


Method G employs the two equalities below during filling: 


P 

2 


P 

3 


P 

4 


P 


/ 


5 


(33) 


(34) 

(35) 


Equations 
can be written 


governing flow 
as follows: 


into the first of the two connected rooms 


Y_ . 


P 


1 


Y-l p ± 


P 

Y _2 1 

7-i' p 2 + 2 


2 

u 

2 


(36) 


P 


2 




(37) 


To these two equations one of the two possible pressure relations 
represented by Eqs. (32), or (34) must be added to form a complete sys¬ 
tem from which P^, , and for the time interval under consideration 

may be calculated. 


Under Methods F and G, the equations governing flow into the first 
of two rooms are identical to those from which the filling of a single 
room is computed. Calculation of flow into the second or inner room is 
more complicated however. 

Because observed pressure buildup in two connected rooms has not 
generally followed simple calculations based on an assumption of uniform 
pressures over the control surfaces of Figure E-10, provision will be 
made in the analysis at the outset for a correction term A. The quan¬ 
tity A will be chosen to reconcile momentum balances based on inner and 
outer control surfaces such as those shown in Figures E-3 and E-7, that 
is, 
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FIGURE E-10 CONTROL SURFACE USED IN CALCULATING FLOW INTO 
SECOND ROOM 


E-51 



or 




A + A 
2 



P ) 
2 


A 



A 


2 


(38) 


Generalizing, we will take the correction terms, A, equal to the pressure 
differential across the duct times the area of the duct. Thus, the cor¬ 
rected Eq. (33) becomes: 


(p 3 “ V V P 4 “4 A 4 < 39 > 

During filling of the second room, there will be an average wind 
moving through the first room into the second doorway. If u 3 is the speed 
of this wind and A 3 is the cross sectional area of the first room, then, 
in view of the quasi-steady state assumption, conservation of energy im- 
pl ies 


JL- 

Y-l 




2 

u 

3 


P 

_Y_4 1 2 

Y-l P 4 2 4 


(40) 


and conservation of mass implies 


p u A 
3 3 3 


p u A 
4 4 4 


(41) 


and the assumption of isentropy leads to 



(42) 


Eqs. (39), (40), (41), and (42) make it possible to compute by Method F 
the four unknowns: 


P 


4 ’ 


u P and u 
4’ 4 3 
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The resulting solution for can be put in a form similar to that used 
in computing P 

2 


2Y 

Y + 1 


Gt)‘ 


_4 Y - 

P' v + 

3 Y 


/A \ /P \ 

1 I _i \ / _4 \ 

1 " \ A 3 / \ P 3 / 


i j>/ 

1 A 3 \ P 3 / 


l/Y"j2 


If y = — 

* P / 
3 


. v . V 

1 / _4 \ / _4 \ 

1 1 \ A 3 / \ P 3 / 




A P' 
3 \ 3 


and B = 


2Y 

Y + 1 


then Eq. (43) can be written as 


„ !/Y 

By = y + A 


where in general A = A(y) ; however, since P ^ <■ P^ and usually « A^, 
the dependence of A on y is weak; i.e.: 


Use of Eq„ (45) is equivalent to neglecting the average wind speed in the 
first room in Eq. (40) . 


After the duct parameters, i.e., p^, P , u^, p^, u^, and, if desired, 

u , have been found from the foregoing equations, then the new room pres- 
3 

sures P and P and densities p and p are calculated: 

3 5 3 5 
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( -X.!§ + i u 2'L_L 

\ V-l Pj 2 3 j Y-l 


(P 5 - P 5> V 5 


where 


Am 5 = up A . At mass increment in second room 
4 4 4 


P 

— ' — -Am = —- (P - P') V„ + — (P -P') y (47) 

y-l p Y-l 3 3 3 Y-l 5 5 5 


where 


Am = up A At , mass increment in first room 

£ £* Zt 


p = p * + —- 
5 K 5 V n 


Am - Am c 


p = p ' + - 

y 3 h 3 V 


Outflow from both the first and second rooms is treated with similar 
basic equations except Eq. (37) is replaced by 


and Eq. (42) by 




P \ l/Y 


■4k) 


Again, in the computation of outflow parameters in the first doorway 
(1 • e - 1 P 2 > p 2 and u g ) Eqs. (32), (36), and ( 37 ) are used as set forth 
above when only a single room was considered. For outflow through the 
second doorway, two cases should be distinguished: ( 1 ) there is inflow 
through the first doorway and ( 2 ) flow through the first doorway is out¬ 
ward. These conditions affect the form of Eqs. (33) and (40). In case 
(1), u 2 =0, and Eq. (40) simplifies to: 

P' P 

_V_. _3 _ _Y_ . 4 1 2 s 

Y-l P^ Y-l P 4 + 2 U 4 (52) 


but Eq. (33) must contain two correction terms, A and A , to account for 

3 5 

excess force against the righthand and lefthand walls, respectively, of 
the first room: 
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(53) 


( p r P 3> A 2- A 3 + (P 3 " V A 4 +A 5 = P 4 U 4 A 4 
where A g = (P^P') A g and A & = (P' - P'). A 4 . Thus Eq.. (53) becomes: 

P 5" P 4 = P 4 u 4 < 5 

Combining Eqs. (51), (52), and (53), we find: 

1/Y 

By = y + A 


where 


y = F 


_2 y < ^3 

> +1 ’ P 3 <P 5 


P' + A_ 


In case (2), there is net flow through the first room (i.e,, u g / 0) 
Eq. (40) is unchanged but the analog of Eq. (33) is: 

<V P 3> » 2 + < P 3 - V A 4 + K = P 4 U I *4 


where 


A' = (P' - P') A - (P_ - P') A 
5 5 34 1 32 


This value of A' is computed'by equating the force tending to drive the 

air out of the second room, i.e., (P' - P ) A. 

5 4 

to the momentum flux stopped in the first room .and the outside reservoir, 


( P _ p'),a + (P' - P ) A + A' 

X 3 & o 4 Q o 


Solving Eqs. (40), (41), (51), and (58), we find again that 

t_ 1/Y . 

By = y + A 
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where x and B are given by Eqs. (55) and (56), but now 


I P - P A P' + A 

„ 2 X 1 3 2 3 5 2 

(1 - a ) ——;— — + —-7— + a y : 
P' A „ 

5 4 5 


A P' 

4 5 l/v 

a = —. — .v T 
A p' y 
3 3 


Substituting Eq. (59) into Eq. (61), 


Y+l ( , 2 2 1 

A = — | (1 - ar ) + a y j 


Since, ordinarily, A » A , and y < 1, O' « 1, from which we conclude: 

o 4 


An alternative expression for A in case (2) can be found by writing 
the excess force on the wall against which the outflow from the second 
room is directed as the pressure differential across the duct times the 
duct area: 

A "s - <p 5 - P 3> A 4 <60 

Substituting A from Eq. (64) for A' in Eq. (60), we compute for the con- 
5 5 

stant term 

P - P' A "1 I 

. Y+l 1 3 2 2 ( 

A = Z~[ (1 - “ ) -p— ~ + 1 + <2 :y (65) 

I 5 4 l 


for which small a becomes 


P' A 


Y+l 132 
= J -■ -- — + i 

Y-i p; a 


One final case remains to be considered, namely, outflow through the 

first doorway combined with inflow through the second. Conditions in the 

second doorway are computed by solving Eq.,(40) with u = 0, Eq. (42) and 

3 


<P - P') A +(P'-P ) A = p u 2 A 
1 32 3 4 4 4 4 4 
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These equations lead to Eq. (44) in which 



In summary then we see that in every case an equation of the form 
Byl/Y _ y + A must be solved for y. This equation is encountered also 
in the first doorway if Method F is used. In Table E-2 the meanings of 
y, A, and B for each case are listed. Since conditions in the first door¬ 
way under Method F are independent of conditions inside, the definitions 
for the first duct appear separately. 

In Figure E-ll results of a calculation carried out by Method F for 
two rooms are compared with measurement of average pressure in the first 
room. The two rooms consisted of two small models such as those illus¬ 
trated in Figures E-5 and E-6 placed back to back with a connecting door 
exactly like the outside door. The experiment was performed in a shock 
tube in which the wave struck the first doorway head-on. Also in Figure 
E-ll appears the result of a calculation treating the whole volume of the 
two model rooms as if it were in a single room. In Figure E-12 measured 
and calculated results for the second room in the model appear. At least 
in this one example, all three calculated pressure rises, i.e., in the 
first room, in the second room, and in the whole volume of both rooms 
treated as a single room, are quite similar and there appears to be no 
advantage in using the complicated procedures for computing the fill of 
two connected rooms. 

The calculated histories shown in Figures 11 and 12 have not been 
carried beyond the time of equilibrium between inside and outside pres¬ 
sures . 

Figure E-13 shows pressure history calculated by Method D ignoring 
the wall between rooms. The outflow phase (when calculated pressure ex¬ 
ceeds input pressure) especially appears to be without justification. 

For this calculation the discharge coefficient has been set equal to 0.7 
on inflow and 1.0 for outflow. 


* CAVFIL, a FORTRAN program written at IIT Research Institute, was 
used to make the computation 
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Table E-2 




Table E-2 (concluded) 


Outflow Both Doors 


y = 


b = 


A = 


a 


p / p/ 
2y 3 


7 + 1 P 3 P 5 
^ [1 - a 2 (l- y)] 
A p ' 

4 H 5 1 /Y 

- - r y 

a p 

3 3 


or alternatively 


Y + 1 

Y - 1 


' 

"p - P 7 A 

(1 - O' 2 ) 

1 3 2 

-7- -- + 1 

P A 

L 5 4 J 


+ O' 


Inflow First Door/Outflow Second Door 
P 

4 

y = 


B = 


A = 


p' p/ 
2y 3 


Y + 1 P 3 Y 


V - 1 

Y + 1 


Outflow First Door/inflow Second Doorway 

P 

4 

7 = 


B = 


A = 


ix. 


Y + 1 

Y ~ 1 

Y + 1 
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OVERPRESSURE (PSD 



TIME (MILLISECONDS) 


FIGURE E—11 COMPARISON OF OBSERVATIONS IN THE FIRST ROOM OF A TWO-ROOM 
MODEL WITH CALCULATIONS BY METHOD F 
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OVERPRESSURE (PSD 



TIME (MILLISECONDS) 


FIGURE E-12 COMPARISON OF OBSERVATIONS IN THE SECOND ROOM OF A TWO 
ROOM MODEL WITH CALCULATIONS BY METHOD F 
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CALCULATED HISTORIES (G) 

-ALL VOLUME IN ONE ROOM 

-FIRST ROOM OF 2-ROOM MODEL' 

-SECOND ROOM OF 2-ROOM MODEL 

-. CAVFIL (ALL VOL. IN ONE ROOM) 



UNCHOKES 


16 24 

TIME (MILLISECONDS) 


FIGURE E-13 COMPARISON OF OBSERVATIONS IN TWO-ROOM MODEL WITH 
CALCULATIONS OF METHOD G 
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The outside pressure, labelled "input history" in Figures E-ll, E-12, 
and E-13, was measured in a separate experiment. It is in fact the stagna¬ 
tion pressure reported by a pitot gauge, placed in a shock tube having no 
other obstructions, when the incident wave was essentially identical to 

19 it 

the wave which impacted the cube. The source of the curve labelled ex¬ 
ternal history" in Figure E-5 is similar. Thus the short-lived re¬ 
flected wave has been justifiably neglected in these "external histories." 


IV Openings into Different Pressure Fields 


The simple methods considered under the heading General Case in Sec¬ 
tion II above can easily be extended to the case of simultaneous flow into 
a single room through several openings, each of which is exposed to a 
separate outside pressure history. This situation will arise, for example, 
when a nuclear blast wave sweeps over a building, striking one of the four 
walls head-on, two side-on and exposing the fourth to the wake of the wave. 
The roof would be struck partially head-on if the explosion were an air 
burst, or side-on if the burst were on the ground. Head-on impact pro¬ 
duces a strong but brief reflected wave followed by a quasi-steady wind 
and an associated strong drag force superimposed upon the static pressure. 
When impact is not head-on, only the drag and static pressures are ordin¬ 
arily taken into account. References 1 and 3 contain discussions of methods 
of estimating outside pressure in these cases. 


A. Computer Program 

A computer program has been written in time-share FORTRAN to calcu¬ 
late average pressure as a function of time inside an open room exposed 
to a nuclear blast wave. Dynamic pressure in each opening may be obtained 
also. Characteristics of the room, the wave and ambient conditions are 
typed by the operator in response to typed questions. Provision is made 
for a maximum of eight openings and the delays between wave front arrival 
at the leading face of the building and each opening must be supplied. 

Only three different pressure histories, however, are available for as¬ 
signment to each opening: the history associated with head-on impact, 
and two others which contain only the static pressure behind the front 
and a drag force. Three different drag coefficients must be specified, 
one for the history associated with head-on impact and one each for the 
two remaining histories. In order to relate each opening to one of the 
three possible histories a "location code" is specified for each opening. 
This code is the number 1 for the head-on history and the numbers 2 and 
3 for the remaining two histories. 
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Duration of reflected pressure ("clearing time") against the wall 
struck head-on is not calculated in the program but must be entered by 
the operator. Other details of blast history are however calculated ac¬ 
cording to algorithms based on descriptions given in References 1 and 3.. 

The operator may choose certain features of the presentation of the 
program output. He may have outside (or "input") pressure histories at 
each opening listed (pressures as functions of time); he may limit pro¬ 
gram output to peak room pressures (both static and dynamic) or he may 
choose to have only the history in the room during pressure build-up pre¬ 
sented. Finally he must indicate whether or not he wants dynamic pressure 
histories in each opening presented. 

It should be pointed out that the state of knowledge at this time 
does not justify great reliance upon all parts of the numerical solution 
to any given problem. Values of dynamic pressure in particular have not 
been compared with measurements. However, calculated values are conserva¬ 
tive for purposes of damage assessment. Calculated values of room pressure 
during outflow through one or more openings probably do not accurately 
reflect the relatively small actual pressure differential between inside 
and outside. 

The program algorithm is based on Method F. A partial value of Am 
is calculated for each opening corresponding to each time step At and in¬ 
side pressure P 3 . The sum of the partial values for all openings becomes 
Am. The equation 


1/v 

By = y + A 

is solved during outflow in a subroutine. When two solutions for y < 1 
exist the larger is chosen. 

A listing of the computer program appears in Table E-3. 


B. Numerical Example 

As an indication of how the foregoing general procedures may be 
modified to account for several openings a brief numerical example will 
be carried out here. 

Consider a volume contained in two rooms connected to each other by 
an open doorway and consider each room connected to the outside by a 
single doorway, as shown in Figure E-14. The blast wave that sweeps over 
the structure will be characterized by a free field overpressure of 10 psig, 
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BLAST FRONT 



FIGURE E-14 SKETCH ILLUSTRATING NUMERICAL EXAMPLE 


head-on incidence upon the opening to the larger room, and positive phase 

duration of 1 sec. Standard ambient conditions will be assumed, i.e., 

P = 14.7 psia and o = 0.076i lb/ft^. In the front wall the opening 
0 2 ° 

area is 25 ft and there is a like opening in the rear wall. Because the 
door between the rooms has an area comparable to the area of each of the 
outside doors, the presence of the partition will be ignored and the 
volume to be filled taken as 


V = 50 x 20 x 8 = 8000 ft 
3 


As a first step the pressure histories outside the two doors will 
be calculated, according to the procedures recommended by Ref. 1. Since 
the first window is struck head-on, there will be no time delay there 
and the peak pressure will be the reflected pressure which is calculated 
as: 


E-67 



Here 



7 P + 4 P 
o _so 

7 P + P 


o so 


P =10 psi and P =14.7 psi, hence 
so o 

P = 25.3 psi 
R 

This reflected pressure will be felt in declining strength for the dura¬ 
tion of the clearing time, t , which is estimated as 


c c 

where s is a dimension of the wall undergoing pressure clearing and c is 
sound speed. In this case 


= [y P 0 /P 0 ] = [1.4 x 14.7 x 32 x 144/0.076096] 1/2 =1116 ft/sec 


and s = 20 

so that 

t = 53.7 ms 
c 

As an approximation the decline of reflected pressure during the interval 

• 0 £ t £ t 

c 

is treated as linear, that is, at t = t c the pressure P c on the first wall 
is simply the sum of the free field overpressure P and the drag P 
arising from the winds behind the shock front. Since during the time t 
free field pressure has fallen exponentially 1 from P =10 psi ° 
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where t is the free field duration of positive overpressure; in this 
case t Q = 1 sec so that 

P =8.97 psi 
sc 

Peak drag pressure P is computed from the formula 

do 


( P so) ! 


5 \ so/ 5 100 

p = — . - = — • - = 2.21 psi 

do 2 7P+P 2 7 x 14.7 +10 

o so 


so that drag pressure at t = t is: 1 

c 


( c \ (-2 t /t \ 

P =P ( 1 - — eV co/ 

dc do I t / 


Numerically this is 


P =1.78 psi 
dc 


Therefore the pressure outside the first opening at t = t is: 

c 


P = p + p = 8.97 + 1.78 = 10.75 
c sc dc 


and, assuming a linear fall from P so to P c , for 0 St St, pressure P- L 
outside first wall as a function of time becomes: 


P = P + 
1 c 


(*c - «) 


( p - P \ 
V R c/ 


For the remainder of the positive phase duration outside pressure at 
the first opening is simply the sum of the decayed side-on pressure 


P = P 
s so 


i .i\ e - tA o 

*o 


and decayed dynamic pressure 
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so that for t > t > t 
o c 


p 


1 


p + p 
s d 


( 70 ) 


(71) 


For the first wall a drag coefficient equal to + 1.0 has been tacitly 
assumed above but for the second opening this coefficient will be different 
from 1, and according to Ref. 1, a value of -0.4 may be assumed. At the 
rear opening there is no reflected pressure and the delay is equal to the 
time taken by the front to traverse the building (assuming that the open¬ 
ing is already open upon blast arrival or that it is immediately forced 
open by the blast). Blast front speed can be found from the formula 1 


U = c 

o 


1 + 


Y + 


- P 


so 


1/2 


where c = sound speed in ambient air, i.e.. 
o 


c 

o 


y P 32 x 144/p 


o 


1/2 


1116 ft/sec 


Hence 


U = 1116 x 




1 + 


2.4 

2 .8 



1/2 


1148 ft/sec 


so that delay 


at the rear entrance is 


50 

1148 


43.5 ms, 


Beginning at t = 43.5 ms the room starts to fill through the rear opening. 
(Outflow through openings other than the first can ordinarily be neglected 
during the delay period: either the other openings are closed to the blast 
for a certain period or, if not, the blast travel time to them is much 
shorter than the time required to start outflow through them.) Filling 
through the rear opening takes place however from a reservoir at lower 
pressure than that outside the front opening, i.e,, outside pressure P 
at the rear is: 



0.4 P 


(72) 
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The decline in P go and P do which occurs while the blast front travels 
from front to rear opening is negligibly small and may be safely neglected 
for buildings of ordinary size. 

For this sample case, the quantities P-^ (pressure outside the front 
opening) and P-^ r (pressure outside the rear opening) have been calculated 
as functions of time and plotted in Figure E-15. The figure shows the 
discontinuity in the derivative of P^ with respect to time at the point 
(t c , P c ) when the reflected pressure is assumed to disappear and the out¬ 
side pressure takes on its quasi-steady value. 

Also plotted in the figure are inside pressure histories P 3 calculated 
by two methods: the greatly simplified procedure given in section II-A 
of this report and the step-by-step method F explained in Section II-B. 
According to the first of these two procedures the estimated history is 
given by the line segments ODFG obtained as follows. When the blast ar¬ 
rives at the front opening, filling immediately begins along line OA, 
where point A is the intersection of the outside pressure history and the 
abscissa 


V 

3 _ 8000 

2A^ 2 x 25 


160 ms 


Line BC is a similar line representing filling through the rear opening, 
beginning after the delay time of 43.5 ms. Ordinates under the line BC 
have been added to ordinates of OA to produce the line DE. Since areas 
of both openings are equal, the point F is placed halfway between current 
outside pressures outside the front and rear openings and the decline of 
P 3 represented schematically by the line FG. 

The step-by-step calculation results in the curve labelled "Method F" 
in Figure 15. Because of the high reflected pressure during the interval 
0 ^ t ^ t c (which does not influence the results of the simplified method 
in this example) the more careful calculation shows a faster build-up of 
room pressure than the line ODF. To demonstrate the method the first step 
of the stepwise solution will be calculated below. 

Since the least sound transit time across the room is approximately 
20 ms we will choose a value of At = 5 ms. At t = 0 there is only one 
opening, that in the front wall. Outside pressure there at that time is 
P R = 25.3 psig. Inside pressure P 3 = P =14.7 psia and density is 
= 0.076096 lb/ft 3 . 
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1 . 


p 


1 


25.3 + 14.7 = 40.0 psia 


2 * P l = 


~p 

1 

1/y 

40.0 1 

p 

0 

P o = 

14.7 J 


1/1,4 (0.076096) = 0.156 lb/ft 3 


3. P = 0.1912 P = (0.1912) (40.0) = 7.65 psia 

2 1 


4 ‘ P 2 = 


5. u 2 


1/ Y 3 

p = (0.307) • (0.156) = 0.04785 lb/ft 


2 Y 


P 

2 

2.8 

© 

o 

7.65 

Y " 1 

_ P 1 " 

P 2 _ 

= 0.4 

0.156 

0.04785 


32 x 144 = 3.11 x 10 


6 ft 


sec^ 


Cl —- — U 0 — * . x... a -i -\j rt 

m 21 2 2^2 sec 2 


4 lb/ft 

= 7.45 x 10 - 7T (dynamic pressure at front opening) 


7. 


Since P > P 7 : u > 0, i.e., flow is inward. Were P 7 > P , u would 
1 3 2 3 1 2 

be negative 


8 . Am = u p A At = 0.0105 lb 
31 2 2 2 


P 1 ^31 (40.0) (0.0105) 

9 , = --- = -— - • 144 = 1360 ft lb 

31 


(^) ’i 


(0.156) 


10. Since the rear opening is closed at this time: 



Aw 


32 


0 


Were the second opening available, steps 1-9 would be repeated 
using initial outside pressure at rear opening to calculate Am 
and Aw 32 

11. Am = Am + Am = 0.0105 lb 

3 31 32 

12. Aw = Aw + Aw „ = 1360 ft lb 

3 31 32 
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13. P = P' * - X > = 14 . 7 4 wanisgo) 

3 3 Y 3 - 8000 

pressure) 


14.768 psla (room 


14. 


P 3 

in room) 


Am 

/ 3 „ 0.0105 3 

P~ + “ - 0.076096 + = 0.0760973 lb/ft (air density 

*j V oUOO 

o 


15. Set P 7 equal to P 
3 3 

and p 7 equal to p and return to step 1 with value of 

O O 

P at t s it = 5 ms. 


V. Edge Diffraction of an Acoustic Wave 

A weak planar shock striking a semi-infinite wall head-on can be 
treated approximately as a self—similar acoustic wave in the manner demon¬ 
strated by Ludloff , 16 Such treatment is two-dimensional and neglects the 
presence of floor and ceiling. 

In the acoustic approximation all disturbances or effects are propa¬ 
gated with sound speed. Thus after the incoming wave front strikes the 
semi-infinite wall, the influence of the edge will be felt only within a 
cylinder whose axis is the edge and whose radius is 


c t 


where c is sound speed and t is elapsed time after initial impact. 

If a cartesian coordinate system is placed so that the edge becomes the 
z-axis and the negative y-axis lies in the wall the location of the circle 
of disturbance due to the wall is 


2 2 2 2 
x + y = c t 


and the equation satisfied by the overpressure p is 



5* ay 


2 

c 



If a change of variables is made: 


( 73 ) 
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X 


(74) 


T1 = 


x 

ct ’ 


ct 


and tan 0 = 


y = H 

x CT 


then in the new coordinates (7], a, 9) the disturbance is confined to a 
circle of unit radius. If, further, radii are changed in scale according 
to the formula: 


P = 


1 + (1 - 


2 1/2 
r ) 


(75) 


where p is the new radius and r =7] + CT , then the equation satisfied 
by p in the cylindrical polar coordinates (p, 0) can be written 


P 



\ 2 

ap\ ap 

SP J 3S 2 


0 


(76) 


This is Laplace's equation and is satisfied by the imaginary component of 
any analytic function of the complex variable. 


In the case of edge diffraction of a weak 
a solution for 0 ^ r ^ 1 (or 0 £ p ^ 1) . 

- — to —, as illustrated in Figure E-16. 
2 2 


shock the present method yields 
The angle 0 lies in the range 

The origin of coordinates is 


labelled 0. 


The boundary conditions are determined by physical considerations. 
Because the acoustic Eq. (73) or (76) is linear in p, the incident over¬ 
pressure may be taken as unity and the pressure reflected from the wall, 
as 2. The overpressure in the undisturbed air is of course zero. These 
conditions imply the following pressure values on the circumference of the 
circle of disturbance 


TT 


— ^ e ^ o 


0 £ 0 < 


TT 


TT ^ 9 < 


3tt 


j 

p = 2 t 


( 77 ) 


The areas of uniform pressure are marked in Figure E-16. 
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FIGURE E-16 BOUNDARY CONDITIONS ON CIRCLE OF INFLUENCE 


Potential theory guarantees that there is only one pressure distribution 
within the unit circle which satisfies Eq. (76) and meets the foregoing 
boundary conditions. Solution of Eq. (76) will be guaranteed by the fact 
that the distribution function is analytic. 

A further transformation of variable will make the exposition of the 
correct distribution function clear. Let 



In the w-plane the unit circle of disturbance becomes a semicircle of 
radius 1 ; the back side of the wall lies along the line $ = 0 while the 
front^ side falls along $ = tt. The line 0 = 0 in the £-plane rotates to 
$ = -5 and the line 6 = tt is found along § = Thus in the w-plane the 

boundary conditions on the semicircle R = 1 become 
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Any analytic function of w will also be an analytic function of Q. 

TT 
i 4 

The function w - e is represented by the line BA in Figure E-17 
and its argument by the angle a. When A falls on the unit circle, a in- 

TT 3 IT 

creases discontinuously from - - to — as A moves counter clockwise through 

• 4 . 5tt 4 

the point B. The function w - e 4 is continuous within and on the upper 
semicircle DBE. Furthermore the included angle y is 


where Arg 


. TT 
1 4 

w - e I = - 


Y = O' + (3 


a and Arg w - 


= p. But along the arc 


EB, v = —. Hence the function 
1 2 



is zero along EB and 1 along BD. 

By a similar argument based on Figure E-18 it is clear that the 
function 



is zero along the arc EB ' and 1 along B 'd. 
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Thus the sum 



( 80 ) 


meets the required boundary conditions along the semicircle EBD in the 
w-plane. But 


F (w) = 


- < 

TT 


In 


w - 


. 5tt 
1 4 

e w 


7tt“1 


- e 


w - 


TT 

i— 

4 

e w 


, 3tt 


- e — 1 


- ITT 


is analytic in w and z and since f^ + f 2 in Eq. (80) is evidently the 
imaginary part of F(w) , Eq. (80.) provides the sought-for expression giving 

the pressure distribution within and on the circle of disturbance p = R = 1. 


Eq. (80) can be written as 


( R sin $ + . 

J2 \ _ -1 

-x_ . I + tan 


R cos $ - 


■fi 


R Sln § ' J2 

r cos $ - 

J 2 


- tan 


-11 R siJ # - S 

R 00 s ' - J 


- tan 


-1 


R sin $ - 


V2 


R cos $ + ,- 

V 2 , 


- TT 


(81) 


where 


R = p 


1/2 


2 1/2 

1 + (1 - r ) 


1/2 


and 


s Tf 6 

$ = — + — 

4 2 


The denominators of the four arctangents in (81) have real roots in the 
range 0 ^ ^ TT when p 2: 1 /2; and care must be taken to evaluate the in¬ 

verse functions in such a way that p is continuous within the unit circle. 

1 

There are no zeros in the denominators when p < — . 

£ 
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0 and 


If R cos $ + . r 

1 a/2 


R cos t a - ^ 


then — £ $ ^ tt 

4 1 


and 0 ^ $ £ — 

2 4 


If the FORTRAN algorithm is used to evaluate the arctangent terms in (81) 
then for 


the second term, 


§ > and R > 


/ R sin $ + ,- 

-1 / a /2 


R cos $ - 

a/2 


may be increased by TT, and the third term 


-i / R sm * - 7i 


R cos $ - 


must be decreased by tt, or the sum in (81) increased by 2 tt. When 

$ > and R > J 2 WG mUSt add another 2tt to kee P the expression for p 
continuous and insure the existence of its derivatives. 


In order to assure continuity with respect to radius the foregoing 
choices limit the arctangents to certain quadrants where radius R is such 

that no zeros in the denominators exist, i.e., when R < —. Thus, since 
we have chosen to add tt to ^ 


- 1 ,,Mt ./2 

\ R cos | - r 
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for § > $ and R > —, we must, when using the FORTRAN algorithm, add TT 
2 2 ^ 

to the same inverse tangent for all $ when R < — in order to place the 
angle computed in the second quadrant. Similarly, by subtracting tt from 



1 TT 

when $ > § and R > - we are placing the arctangent in the range - — to 
2 2 ^ w 

- tt; hence the same function, when R C -, must be computed in the same.^ 

range by subtracting tt from the angle as computed by FORTAN. When R < - 

the denominator 


R cos $ + .- 


is always positive and when R < — the FORTRAN choices of quadrant for the two 
arctangents containing it are consistent with those for R > —. 

Pressure contours, located as outlined above, are shown in Figure E-19. 
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FIGURE E-19 PRESSURE CONTOURS WITHIN CIRCLE OF DISTURBANCE 


E-82 


NOTATION* 


A area of opening or openings through which room fills 

A constant equal to various functions of y, appearing in 

Equation (5) 

B constant equal to various functions of y and/or of interior 
and exterior conditions, appearing in Equation (5) 

c speed of sound in air 

c v specific heat at constant volume for air 

F numerical factor for consistency of units 

I integral 

K discharge coefficient 

i distance along jet axis measured from opening 
m mass of air 

Am increment of air mass corresponding to time increment At 
P air pressure 

v' air pressure during immediately preceding time step 

P-^ 0 peak overpressure outside front or first wall 
P lr total overpressure outside rear wall 

P c sum of dynamic and side-on air pressure at time t= t c 
P d free-field dynamic pressure 

* Numerical subscripts may be used with variables to refer to a 
particular space or volume; see Preface. 
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P dc dynamic pressure in blast wave at time t= t c 
P do peak dynamic pressure in blast wave 
P Q ambient air pressure 

P R reflected pressure at a wall struck by blast wave 
P s free-field side-on overpressure 

P gc free-field side-on blast overpressure at time t= t c 
P go free-field peak side-on overpressure 

p acoustic (over)pressure 

q dynamic pressure of air 

q (spatial) average dynamic pressure in jet at one cross section 
R gas constant for air (See Ref. 5.) 

R radius of jet of streaming air 
r radial coordinate, i.e., r = T[ 2 + ct 2 

r radial coordinate measured from axis of jet 

S specific entropy of air 
AS increment of entropy 

s wall dimension used to compute clearing time 

T absolute temperature 

AT filling time or interval between first arrival of blast and 
achievement of pressure equilibrium 

t time measured from first arrival of blast at structure 

t c clearing time of structure in blast wave 

t Q duration of positive side-on overpressure 
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At increment of time 


U blast front speed 

u particle speed of air 

u* air particle speed along jet axis 
V volume of room or other space to be filled 
AV increment of air volume 
v speed of object in air stream 
x Cartesian coordinate 

y Cartesian coordinate 

y unknown variable in Equation (5) 

y o a specific solution of Equation (5), i.e., y Q = 0.1912 
w complex variable Re^ 

Aw energy increment 

a angle in complex w-plane 

a a certain function of y 

a acceleration coefficient of object in air stream 

3 angle in complex w-plane 

Y ratio of specific heat at constant pressure to specific heat 
at constant volume 

Y angle in complex w-plane 

A correction term applied to momentum balance associated with 
one control surface 

Q complex variable p e 1 ® 
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T] time independent coordinate, i.e., T| = y/ct 

6 angle between inward normal to surface and the positive direction 
of the x-axis 

8 angular coordinate, i.e., tan 0 = y/x 

p density of air 

p reduced radial coordinate 

p ' air density during immediately preceding time step 
P Q ambient air density 

<J time independent coordinate, i.e., c = x/ct 

T time required to transmit a sound signal over the longest room 
dimension 

/ 

cp angle in complex w-plane 
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F^tr va/an/ S/s?y/a- Z/eyraa -e~/~ Freac/om *fyf4rm 

Mass /a c/or^ — O. i> ** 

/oaS /sc^ 72/_ - 0.44 *" 

Farir/ama -Ac/or 4/ e ~ k/ L - 0.44 * 


e e £ 4 
r-4/6) 


^ - *2 K m • OZ./S3)f.s) = 4.09 / S *- src 2 ////,». »;&A 
kz = k ^ = ( 4 . 4f/Sx / 0 +)a*h 2.874£4>x/4>*-*/4t 


ct/$c? 7~4 L&/c /si 0 
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T- — 27T 


„ 5 

277 Yz.£7S&k/0*- 


T? - O- £>9/ sec 

£ = '• ° 

= _^L£ - /o. 99 
7^ -*>9/ 

fbr ^U. = 3.0; 7^'" e = O. £7 

= (/XStJ Kl 


f* = 


( 7 S.S 4 ) X u 

O.S 7 


J&L. 


/S. £■4 
O.S7 


= /<£. 2 ps c 


£$.9-/ 


T/ie^e-/■£>'’€, a/er/^s? 'S ac/ee t&Ye £/?ce 

t/ie ees/j-/#/*ce /egc//ee <2 -/or Syne-/n Sc. /e>er/ 2s /ess 
//a*, / Ae e> es/s/orrtre /~e.^c<c r-ec/ /or- aSes/y/y ('see 

step 4J. Sec /' £>*■/ PV/ // /-ye/ 4>e neSa-^/y ne at /or 
p/ise/soseS <?/ cost <es// sm? 7/ ay. 

8. 22e Aoc/eyS 


’ 792 y sooJ/J c £ 

_J2 2.9 2 -< / 2f 

(42, Soo) (2 o)('S. s) 


- £>.08/ sec 


Sc/./jp 8-// 


v _ / o 

T 2 . 08 ! 


= /z. 3 
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■X 


gcr = 3.0 

(f>‘ 2S O. 25 ^ * 

<f>' = (. 2 S)(?.o) - O.SO 

Aj~ Coos){/2)(/S.s) 

A; « o. 93 c\/ft 

Vsc #S(P /o"; Aj = O.9S0"/s? 

*?■ 5/<zA 5u'/nma.e'tj 
■fc - 3ooo fsc. 
d «= /£.S" 
t = /7.0 " 

Z = 22-// {<2/eras’ Sf4*t) 

Aj — ^// @ 3" M//9A evxrp o/Aof s*e Aar- At-st / ts 

■for- e/iiXJ gpsu}/ /esiSror) . 

Aj — ^ S @ /O' J /xfs?~t 0/7 (9 a si c. Aa> oee/ /n Ae>/A 
JOfjPe’n 9 «s*//s. 

As — £>.0020 9 A /Zcf. 3 

= €>. A/ c"//9 

Uaz * 4 - €> //% " <s ’cZc/t -/ace 

Wt. re A ass = (s.S/sJ^j + (2 A 7 °)(-f 3 )-/-(• 

= /7.3£ ps/ 


f-a/- s/snp/e scs/0/?t>s'9s t 


_r_ = 

'Si " & 
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JT. EXT£ZJOe tV/UL 

At. l—oaX/sij —/VaAA A^as/ysi A 

Assosve £A/>A AAo ey.Aes/es cAystasi/e JtC/A />/”<?rfc/se <te£s e>S 
a, /Scr'yi A) </<s A^/ A^fecA AcacA cs> AA<r iascAA^ iaz/A/cuA 
a.7 f fe#t'&X'e>~ «,/'£A (AeyAA. *5 /sice. £ 4 e+s /s s>s> /sAss/es 

ss tv- c Ao/AaA~e/y /A £As £>4ss^is/* £, /As A os/syostAo / 
c/ys)fl,/rj/c jD/'tsscsse s>/7 £Ae> k/oA/ /j £4Ac*? as 

K 0 f v te£/,f.s--S4 


J / )i is assume/ £Aa A /£ 0 « £>. Sos- a,* c/s7S4.-/cr*Ss<A 
ct>A<?srA<? so// 0 / +?reo/,‘u#>? £0 sA/// ct>*ss/fA<=*c*/j £‘t£i j /’- A~£> 7 

ana' £A a/ -fee, A Ass, £ A<z yeoA Aos/^es Ag/ se>; A 

yts'essc/rs d/t/t r £0 AA e AAasA //>a oA/*7^ /s 

^ = 6*X' S ) = 7 ' S />**■ 


7 ~Ae s Ac A/e sc/S yeessese a£ /As Ay 0/ AAs h/o/A^ 
ih/AA} AAs £y> &£ A A/ a- A /As A /Aces’ 3 At} A 4. A ysa/r^ /J 

/ n)£ I Ac£ £ 

y?. = l 3 °AUL =r O. 2 / X>S( Sec. 23 'C 

r* /AA ' 


TAe sAaA/c sc// yross^e o. £ AAc AcAAo/n 0/ AAc? 
L/joA/j Acs 0. c/y /A £0 AAer £ oy c A AAe £oc>A //7 y 
»/ /£'-//"; JS 

= = A/A ps£ 

/AS¬ 



IA / £ /’s assa/*>ao/ £Aa£ £ As s/a he so// yssssc/ss ' s 

Sda/ t/a/seA Ao O' t**i//>ssrt Ay e/ijAsi AcAecA (Ays? d/n/c Sfi/A 


yressveo 


f 


/%- f- is* ,,i. 


/go./:/) />.£-£ 
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Men AA<r />e«A f reside ^ c/es*f» *>/ /A<? 

iA/a// / s 

^ = 7^/ 4 - = g.9 />si. 

ff. /-^/^e /im/ /i JA/& A/ 3/esi^s? ^ / 

A// e.t'fericr' i^a//s hs'AA jc// 3a c A-fs// a*-e a/ 
s^i/ae cc*fsA^vc/;e>*. 2>e:/f» *! r**-*/*#, r^f'/ 

scy>4>£>*~'A&(/ **" AA ~ •$• 

/ geeeesr/ JFe^/iAasc^y J/ee/ A stv***e e*./s3~f wo//(A/aA^es* ,s 
adcgvo/e 4r>J co/co/a/e refvieec/ s/ee/Oreo. ■ See f/efS.As/o* 

3, stA,™ /C Ma/s’s'/*/ /ei 

See f- / • 

3. A/e^csi'-ea/ Pef/sAa^ice 



TAe e/ear y?o*? o/ / Ae wo//, Se/weew £Ar A/ret 
//cos s/4* *»</ AAe /*/> *S AAe AeeS'A? , vArSes 

/2-4" to /3-//"/*<=* /* P-/C). TAeee/oe* assi/Aie a. s/aaei 
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/,<SV7 / ^ ^ /-3- 9 fc//y>£>j<?j. 7~Ae /~ e ^c7^ e /7 

ft-rteof tem/Ae j/ ce / ; s 

^ =• ^°- 7 ' ) //3.7S x/ 2 1 2 

(■ O 7ZJ( 4- 2, ooo) ( /O.s J 

0 C = O.gJ 

A $ = f>A<A = (oe> 87 j(/ 2 )(/o.sJ ~ /. 70 #'//? 

*8 $>8'/z"; 4 - /,/2 er'/Mj & = o.S9 


4. SA «■ Of- CA<?cA 


c/_ _ /£>. S 

4 (/ 3 . 75 ) 0 *) 


- o. 0&3A4- 


gor- A^ — 3£>e& yOSt 

S8f>st ))> &. 9 

T/>rrr£>*’<” sect/os> a.</c?g c/at*? Aoe sAe 

8. A7/a.jcei &/ Te^s/c CAc*cA 

8e>r~ A At? <£>/.) sY/nj us#// *v'7~A ^ A (£) /<3 
= sp'AcA = o./3 

/vAeetr <A - 72 ~ (2 7 - 9 , 78 ‘ 

TAe/-<s-/o/-<; / 


oe. 


Air A /) T-i'C}. 


/> 


77.73 


A?.7S)0V 


o. 00 ,4 & /' = o. // 


9-3 
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^ = //. 9 /ssC > '#■ 7 />■** 


TAcre&re, Scc/t^ *cfej**& </•'*?***/ 


&. Z>y^a, m ic Ana/^sis oS P/*c/t*i''»*ry //3>CJtf/i */ 


T= - 


(/3.7S7.J z) 


(42, 50C)('*89)' /,l ( /£> - 


— Q. O(od sc c 


A S'* 


A// /A As - /• & 3CC -for" c/y/raoi'C ft" A/ 0+7 £>/ /aa cJ/sif; 


a. A*t«/yr>s As **■£ i^acca^/oc/ »4v c*ft 


e s// yOCJCJ ■ *3cc />• ^ 4, 


7 Aousi <7 

= / &_ _ /S. 4 

7 - .0&5 

A~0c ~ 43. O 

p ~ 0.2 5-^ c &f.t,Fy.B- / 0 

£ r tis~£/'ij re//iAc?fCe.rnCn4 /*-! tv&// />/-&?'{Acs 

d> = 4 = 

r .85 

TAerefae, s3^cc sot/ J*c/A// veo/d/?'~o V ,Wc c^stdc^A/e 
ce/sctstd rcs/s/otice, 7A& <ex.is//st? s/<?c / /s 

asset**co/ Ac/sgu#/e. 


F-ll 





^TX-terr/fix Ida// /)<x/x>sxi A Xij 

SJ = 3c&0 p>sl 

d = /-O.S* 
t = /2 " 

= @ <£ A S/sjte ^ — &. & 9(p. / r -<g)^ pxav/(A<s s&jst/■&?£//& 

IxXsA x^/xx/ox^O'stupx/ <Z/ A><^x?<A‘x7^ up CXC xy //) ix>c/ /oxip X’<?A>£fX, 

A'; = *4® /£" 

/Zos'/y ox/Aa/ Aj — Z 4- (a) /2 sac A -Agee ('cx/at/o^J 

JV{. xeAaxs^ (2.670)^ 2J+ 6£)(l2_J p (. £ (gjfe) 

- d.SS' fj-f aZ K/a// 

TAe Ae/jA ox c/?4x yJox/^ e>A A As* <^xAxx/Ox 
w<?// i/t?x/es <Avo -6o cA/AZeoex/cos //i s/<?A A As'cAxxsj 
4 1 Ao//onxJ ' 

(a, ) Airs ex xo Z/cox>/ 

I =(/* '-// ")-(/ L s ") = /3-£* 

(6) SAop J/oxa^c AFoo/n M- 

L =(/4'-/S"J -(/-O'j = /3~/;" 

(c) SAe>p /Foos/i ■* 

L = {U'-//"J-(/-3# "J = /3'-7% " 

(Aj Aesosts (A A A Soxx/dox -j^- 

z = c , ft,, *j - (/ '- 2 y z j" =• , 3 '_ s / 2 " 

* ^/c.u/ a A; 0 xi Zox f* e s/e .p s av ^ AAesx xo^ /S „ e -f 
txc/isc/ct </ /'« AA>s <zxiw?’p>/o, 
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AAecA Pi/a// 2/cssys? 4/ /as Acs/? ce> / / cog/ 

A/a^c /ssi esss? Kir s//'ca/ /ce?cZ access; /■n w?// 

suffas //siy S/a A Acs sys // ssj /Pesos se A/coai ctst s/ 

Is csfcs#/ /a 

/= G.39 <? / &• // /- 2c/ A 

l = (3 2'- 2‘ l )-(7-7“) - 24 ~ 7" f. A-2 

^ - /f.o fsc 

<2 rr sgl = (/f.ej(2< /z)=st os ~^y//s7. S/ss/P/j 

f = a/p,? / -f/S,Sj(24.S83K/z) =466/ */„. */c¥4 

J/ ;/ is assusned' /A#/ A Ac fee?/ oeac/soa access; 
fsPcs- A a asry a/ysi/Z/cOs A c/rcesy ss/ AAc A/a> A ts/ase 
/Ac#, 

/ = /. 39J(S6OS)+(.///466/J = 26 99 *//a. iwfA/A 

AAasu a A Acs- Ksa// asrO /yscj y e P Cast Ae caoc/ac/cG/ 

PAa/ sc ses-A/Ct? / /aac/ a-/ AA/s ssutya/Za</e essau/eZ 

fsc A Si A/y / s?Csf4J C. A Ac ess a// sCS ss/a*rCC a.AevC 

/A# A Aos M Zsss?/ /oac/ os Ay • TAoso/osc c / / j 
asse/sssec/ Acs /Ass c.Ka.ss?f/c PA a/ AAc PA#// 

Pcssys? Af/ / s e? c/c^oe? / Ac s ess's/ AAc. estrs/sca/ 

Ace?/ csr/y, «*td A Ac cc*s AsAc Z A?/Ass?/ as / iscs/sce? / 
/ce?c/ s?S A A CtsZ SK a c/s /sea/s as - 
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/O. /-£>oA/s?q /Ac S AS 


(&■) A/Zovvi} A/c So// ffear/n^ Press use 

Assume sc m^cA/om sZ/AAc/ay ee J#'?cAy e/ay 


S/o, A/ 'd 


yZy^Zo&o /-2oc/-/ psf (max-&OOOfsZ) Pef. S 
' Sec. 2 80 S 

U/Aere /-/ - a/ey>/A Ae/oru m/A/m um c/ey>/A 

2/ /s-/ 'y / 4 ' 

= 2o£>0 -t- (2£>o)(/4) =■ 4 ££>0 f>sA‘ 


Syr> Aurj i c 

A//o VSO A/^ </yS 1 &*-T//C. SO// Aepr/Uiy / s A4 A**"*? 

(IS /ov/ce AAe jAoA/c uo/ue y>/us /Ae yep As 

/roe- A/e/c/ sse// / &f'<rssvs'<P* 

-A> - llKt^sL + /S ^ £2 />jL 
/<* A r 


AeA. / y y>. * 9-/4 


(h) 7/r/a/ A~o>ec<r esu, M/#// Aoo A/e>y 


S/zA 


A - /7 y l = 2 4-7‘ 


UA. = A/7j(/){'*) A/4S) (2AS£3 1 2 /0 &//*. uvo/A (v/a 

/723 12 / 


/-■/. — O (Assume A removec/ Ay A/<?sAJ 


PVa.// 


A - /2"; Z = /3- 6 " 

D.l. = ( /l )(Q(' z ) A/As)(/.?.s) = / A£ ^//ei, us / A A A 

/72S 
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Pysiatir/ic £^cj/ra/<Srt~/ Re>n 22. 2. RcR. / /? 9 4 

p, - ( 2/0 * H3)(£.) = 3// */a>. "//RAJ 

7~£>/g/ iRrf/i'ce/ Lcod on W#// /-ee>-//s?g 

/y- — R ■/- 1/ (n3j/tc "//sty Roo R/ny 

4 - 3 // 42099 = 30/0 ^//n. k/c/RA. 

(c) PV#// £bc'S/**y /1//</^4 

0*je e. C*-//nocui Iaj 4// sysftfcA -/’eoR/nJ . 

Re^o/n ed w/c/AA/ 

<■£• W-*’- 


2sc 1 - 3# 

yemscsre yA>, = -LRR'.S- - ~R 9 R J( - 


(d) tVf// PeoR/ny 7~A'cAs?e S 5 

Sfcif) Rnom waAA -Race 

_ 3 8-/2- _ /^ M 

2 

p .C.C,6/ e 


&ff j/ p.8-*7 




<A = 6.£" t -A2 


$? ~ 


f.O/e)(42, coo) 


—■) = 0.244 

£. 3 / 
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A s - ~(.0024-t){/z)(£. sj - 0.25 o'/U/ 

Use *4 @ 9 A ", 4 <0. 25 a/jj 

(^) U/'a jcsia/ 7<?s/r/£>+i C /sec /b 

ft)** 2S\/£~J‘ (t) &£/,/>. 

£-ee ^ — 2e>oo ^jc 

/y — (2 5 ){.2Hk 2oo a )' / i('-&f) 2 
/»JC J> z 7 9 ^ 5 C 


TAee-e/cs-e^ SCC-//0+, sj aUe^uaS -Ap? U/ajssia/ 
6<£s? s /&*?. 


({-) IsUa // /-£>£>// Stssrisr? as*y 


-, {l' = Zooo psi 
d = S.5' // 


t 1 = /2 " 

L - S3" 

As* *4@ 9/z" 


7 s”a*r verse A s = 0 . 002.0 t; & 

= O. 9/ a" 

Use: 3 —# 5" A s - O. 93& " 

Etewe/s J 3-Q /osKp T?C> stiq./c/i us)// s-e/sstrs 

Use ^ S <& 3fa <ZssU ' ?5> V (3) /S " 

Wt reSw ^^sY0)(ff) d {/.o 43)(l)+(** 7oJ(^(s) 

d(.US)(^j (3) 

= t£.45 pUf 

= U.sz />sd 


'-45 
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£*7 


M. INTERIOR S+/ELTEZ MALI 

A. ~ W«// Pes^y, *2 

Assume /7a7 Ate A/asA /oat/Mf ZAe s«/«r,*r sde/drr 

usa//, Aue Ac A A? Ads;d /sressore do,A<d- <y> /« 7/e ror- 
s/or/ed A*-,* *>/ AAe Aas******, 7s cj**d d A Ac 

,A c / Ac* / S7ojd overpressure . dAc c,d so A wo/A 
/c ac J,~f 7S dAe A*-/™/ A/aj/ Aesc/Sr? OLc/to? Ur 

AAe W/ «/*»e w'AAouA coro cdvrfr? ycr/Ccod A****S *» 

AAe */«// dor do od*A reacd,e„s. 

3. Eve/f~, ,r oef do Arr,er SAr/drr W*// &esif» 

Assume ddr </*<>/« pDoryOSJfJ /Aod a// Ar dory or M4//J 
«»c/ 0 tS»f AAe sdc/dee a re a *r-o fJsoj/ar. /><ss/f» as 

or?e- way "*4// */»**/»$ dr/yveer d > 7rsf ant/ 4*te*»«r* 
-/7c Or s/vASj vs 7 A ^AC -S.£>. 

I 3rr,.<r»/ 3e/,/orcrr,J S/co7 J Assume <?c =2.0 (o?A.«.) 3r£//>.*-* 

2. Oyoccrrc Ma/<rr«?7 3eo/> crA/rs 

£ — 42ooo />’<■ j/eisc/tsra 7 jrac/c 

°y r 

/ / - 2 y 7s/o 3 S ^ / dps core red vs>AA / c '~ 3 e CO />-*<■ 

Cf'c 


J?_ 3c^ Of red 3e s/s /an c r 

^ = 7%a fsi. 


1* 


d/ 2 


A (<?.*72g£ 


p>. A ~ 7 

Ae/./jf.S-S 


L =(/4'-0") ~ (/'-2'dd') = 72-9 /z " 

(Madsror, c/ear Sf>a» dor corr/dor */*//)' 
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d - (/2.79Z k /2)f - __) ^ 

I (o. 07Z)(2.o)(420COJ 

d = A 6 7/ 

Ose. </ = <£ 7 ; ft c ~ 7. d 3 

A 3 - d c/ = (e/?3j(,z)(s.7sj= 2. 03 o "/ft7 

17s£- ^ / / (a) 9 " y 4 S - 2. o£ a "//3 

t = S. 7S +- o. 7S■+ = 70.2 " 

2 . 

C7 Je > 70 " 


4. 37>e a*- Cdec/L 

— - -^LZ_£__ C. 0£70 

7- 02.7?z * /2.J 

^ ^°/0-’c !> 7S.S fsL 

77>e^c7'c>s-e J s^c 7/Or? &dsg isa i4 dds jAeao 

~ 7 . / 3 mj) 0 s?!)/ Tens/on C/?ec£ 

3-Or ft' = 0,2 5 ft c ft v - 0.2 S (o>,s7.J 


7?c7 /j 2ia. 9-3 


£?/./j />. £-& 


sgy = 70 0^j~-~jf/+(d2Illf^ft43?£^(/£3x3ooo) // * ftoS7o) 1 


= /J.3 


/9.0 /?SL 


/7evseveOj m/nisz'U/n 


s>.£-d 


sgj = (3.s)(3 e><PG) //z - ft,0S7£>) = 'C.d < /P.O^j7 

77tes-efths-f' j,Tcd/o+i /os odd/c/r^r3 3-or- d/ti^oo#/ /ons/an. 
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T 


/ c 0>s' ^y - ^ 9. ^ / ,Sc 

/9.0 = ^ + 2 ^ 42 - o£> ° \([. 93 k 3ocJ)Y-0 s 7 °) 


= £>. £ 7 

S/sice? — / 93 , / p/~ < #v/<2 / <?= regcr/re c7 sst?7 s'f/'/sAorc ess?es>f 

Jty A^s/cZ/ss^ cs/> cstrsy &//<?*■ /0s>^/'^e/<//>7a/ /"Cdas'. 


£>ys7a/ru'c Asia/ys'f <cf As< e //'rrssiasy Wa // Pejs^si 


__ (P7.7?2x/2j Z 

- £>. 04 6 ' 3 ‘ crc 

/>. >=■-/ 

/V ^ soo)(A 93j Vz 78. 7SJ 

A// A A ^/ = / <P j'<rc 

c/y S 6 "i /'c ^>d>sA/ 0 sj &A 

/<5><S (A/sjJ J 

a o/ysiAsri/c. Asui/ys'S / s 

S7K>/ M/0SS'a» 7e <7 AOs' 

C £3S A 

<S s 7/s*/ 4 A'Si^ A ^"'A g,Jes ' 

Strcr /?. 



~7. TAe&oosio 1 

<6 = AA = 2/.T 

7" • <c ’^ 

7~Os~ ~ -3. O 

^>'c2 0,25^ c O.AS Pef/jP'j.E 

A: = ^'£4 = f.e>o4S)(/z){£.75) - -7/^/Z 

tocr ^ (a> /O^" } As - O.SO <*"/&■ 


■/o 


F-19 



£S. do Aere/oe H/a// 3i/sr>/77ary 

/c - 3 COOpsi 

d - 3 3 A" 


A s = *//<@ ?" zo/'AA eve-s-y o/Afe oeAor A oof csp? 

too d/o peno / te/is/o*). 

A‘ s = ** <£ " 

f-/oe/ep Oe fa/ /f^ — 0.0020 f b 

- < 0.2 A - a "/ft 

(Ase & 4 @> /&'* j e*tcA -/ace 
Wt. rebars = {£. 3/3)(/J:J y- {/■ So2)^dL_j ) 

- 9. <b A p s d 

AAo Ae/aA't oA AAe Jeif<ee</c>e t/i/a// {doe cos/ esf/er/oA/AjJ 

roe/ es doe Ao AAo d/ A/erenc e /A s/oA ( AbcA/iess Q s 
do Afows: 


(<l) £eieroe A? i? oer/ 

L - /3-C" 

(b) AT Aevador dorr/' dor — ddocr/A 

Z = (/*'-//*) ~ {/- O "J =r /3 -//" 

(c) Ar/eva/cr doer/doe — AAor/A 

l = (//'-//'')- (7*/*"J = /d'-3/A 

(d) Ado/A does'/’ 0 l> cr 

<L — ' 3 '-&•/* " 


2- d-AO 


p. d-/0 
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9. 3/> t fc.A lA/a.// Design #2 Son !/er/, ct>//oa.2 

A/la*/snusn sen//eg/ /e ><9 2 occc/rs on /si t e -r /on utsa // 
suppor //rig 2/o/> Des/g si &/ so/ Feses^re Foom / an 2 
/ s eg'uo/ /o 

V = 2 £99 #//n. «///HA p- 2 -// 

This / 04 c/ /s net cr/Z/'ca/j As?2 T/r is-/ 4// /j o sstsss/ec/ 
ac/eeua/O /o res/si i/e iserZ/oa/ Da 2 w/Z/pot n>, o 2 /Fca//or>. 


10 . Foe /isig Design 

(cu) A//o*>aJ>/e So// Bearing Pressure 
P = F 2 psi 


p.F-/2 


(6) To//)/ Force om H/a// Feo/sesg 


5’/at h 


D.l. - 2/0 *//s?. 


/Va// 


p. F-/2 


/ — /o" J Z = /3'-0" 

£>./. = ('°X ' 2 L' 2 L ///£)(/ 3 . 0 ) = /36 */"?■ 

/ 72 S 

P r = 2C99-/ (2/0 +/3£)(f j * 298 7 *//»■ p. F-/3 

(c) For Trig I/// 9 /A 
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Appendix G 


TYPICAL DESIGNS - PRELIMINARY DESIGN TECHNIQUES, 
COMPUTER PROGRAMS AND OTHER DATA 


By H. L. Murphy and J. R. Rempel 


The purpose of this appendix is to present additional information 
relative to the Typical Designs section of Chapter 6 - information that 
includes backup and other data to allow the guide user to extrapolate/ 
interpolate for additional values, delve into other situations, and gain 
further understanding of the approach and techniques used in the Typical 
Designs section, without further enlarging that section. 
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Appendix G 


TYPICAL DESIGNS - PRELIMINARY DESIGN TECHNIQUES f 
COMPUTER PROGRAMS AND OTHER DATA 


I. One-Way Slabs - Simply Supported 

Preliminary Design Procedure 

This procedure is briefly described in the Slab Design section of 
Chapter 6 (and illustrated in Appendix F) , and was used for preliminary 
design of several one-way slabs, based on the following common parameters: 

f ; = 3 ksi; f' c = 1.25 f'; f dy = 42 ksi; P SQ = 15 psi; 

t. = 1 sec; (X = 3. < 

The design steps were as follows: Using a step loading pulse, re¬ 
quired unit resistance q (also termed an equivalent static load) equals 
p /(I - 1/2p) . Flexural requirements were met by assuming d and using 
Equation 8-4* to obtain the slab clear span L. T was calculated from 
Equation 6-4, then Figure B-10* was entered with p and t d /T (t d = t ± =1 sec) 
to determine the required elastic rebound steel ratio p'. Diagonal ten¬ 
sion requirements then dictated the required web steel ratio p v by use of 
Equations 6-9a and 6-9b (first two equations of Figure 8-2b*; Equations 
8-6 and 8-7* are also applicable but in a less convenient form). Pure 
shear capacity proved to be much larger than required, in all cases. 

One-way slabs were similarly designed with p and p' adjusted so that 
the required web steel ratio p v equaled zero. 

Arithmetic plots of L versus d, L versus A g , and L versus A g all 
yielded straight lines for each of the two families of slabs. One small 
and one large value of d for each slab family were, therefore, selected 
for further work (steel was detailed sufficiently to determine slab thick¬ 
ness); descriptive values for the four one-way slabs are: 


* Reference 2 of Volume 1. 


G-5 



Slab 1 


Slab depth d# 


7" 

Slab span L 


130" 

Tens, steel p 

At 

mid¬ 

0.02 

Compr. steel p ' 

span 

.0052 

Web steel p v (required) 

.0091 

Slab thickness t# 

8-5/8" 

Ratio t/d 


1.233 


Slab 2 

Slab 3 

Slab 4 

to 

o 

10" 

27.5” 

360" 

135" 

360" 

0.02 

0.014 

0.014 

.0052 

.0029 

.0029 

.0091 

— 

— 

23-5/8" 

11-9/6" 

29-1/2" 

1.148 

1.156 

1.072 


Based on the steel detailing of these few slabs, approximate values 
of slab thickness versus depth were developed from slabs 1 and 2 into an 
upper limit equation for rough cost estimating and from slabs 3 and 4 into 
a lower limit equation that can be considered in estimating mass that is 
used later in equations of motion.+ The resulting equations may be useful 
until a more complete study is made; they are as follows: 

t = 1.27d - 0.006d^, but t not less than l.ld (for rough cost estimating) (6-7) 
2 

t = 1.20d - 0.005d , but t not less than l.ld (for least mass estimate) (g-1) 


Modified Preliminary Design Procedure 

The preliminary design procedure described just above may well provide 
designs adequate for final design use - a few design comparisons are included 
below - but there were three changes in the procedure that were recognized 
and may be used to improve the resulting designs (that is, as a part of 
Step 7 of Reference 2, art. 9.2), to make them give results closer to those 
from a final design technique that includes all refinements reasonably 
available at this time (Chapter 6). 

The blast load pulse assumed to apply to one-way slabs over base¬ 
ments has a zero rise-time and an exponential decay from its peak value 


* Such odd dimensions would certainly not be used in practice; however, it 
was the t/d ratios that were sought for rough cost and mass estimating 
purposes. 

t Equation 6-7 was used for both purposes in the typical designs in 
Chapter 6. 
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(Eq. 3.51.1 of Ref. l). The convenient chart solutions for single-degree- 
of-freedom dynamic systems, Figure 6-1 (p. 11-3), are exact only for 
triangular load pulses; however, techniques are available for using the 
chart solutions for decay pulses consisting of two or three straight 
lines (Figures 3-8 and -9*), thereby more closely approximating the 
exponential decay pulse. Some comparisons were made between the results 
from triangular and exponential decay pulses, using the positive phase 
duration times from Figures 3-7 and 3-5, respectively, of Reference 2. 
Solutions were obtained using dynamic load factors, Table 6.3, and the 


Newmark (3 Method 21 

with a one msec time 

interval 

and g = 

1/6. Using a 

single problem example about in the middle of a 

range of 

one-way 

slab 

designs and a peak 

blast overpressure p 

so = 15 psi, the maximum 

resis- 

tance q and ductility p values varied as 

follows: 




Umax(psi) 

23.6 17.0 

16.0 

15.8 

14.5 

13.8 

^triangular 

1.24 2.69 

3.43 

3.71 

6.27 

8.83 

^exponential 

1.25 2.75 

3.57 

3.88 

6.86 

10.03 


.995 .977 

.961 

.956 

.914 

.880 


These tabulated values vary over a range of p values, but are related to 
one midrange beam design situation. The following tabulated values com¬ 
pare |jLt with p e for a range of preliminary one-way slab designs, but with 
p varied only over a narrow range (again, p go = 15 psi): 


d ;L 

7" ;138" 

20" ;398" 

7" ;151” 

20” ;440 

f c ;f dy (kSi) 

3 ;44 

3 ;44 

5 ;52 

5 ;52 

Umax(P si ) 

16.2 

16.0 

16.1 

15.5 

M-triangular 

4.31 

3.42 

4.41 

3.71 

^exponential 

4.40 

3.57 

4.53 

3.92 


.979 

.960 

.974 

.945 


Using the 10 pairs of values tabulated, a curve-fitting equation for 
"adjusting" pt (calculated using a triangular decay pulse approximation, 
tgo °i Fig ur e 3-7**) to an estimated p g for an exponential decay pulse was 
developed ; it is reiterated that the equation is for estimating only, and 
then only in the particular realm of applications from which it was derived 


* Reference 2 of Volume 1. 

t Using more precise values of p-j. and p g than those tabulated. 
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P e = p t /(-.012747 p t + 1.00822) 

where significant figures shown are simply those resulting from a library 
computer program run, not an implication of accuracy. The equation is 
most useful with computer programs; a simple figure is generally more 
useful (Figure G-l). This change in p is the first of the three changes 
mentioned above. 

The second change that may be readily applied as a part of "Step 7" 
is to use a "weighted" dynamic load factor. Table 6.3, tailored to the 
p, desired. In review, it should be recalled that only one (combined) 
dynamic load factor K M /K L need be used, .applied to the mass as follows: 

If the fundamental equation of motion 

p(t) - kx = ma 

includes dynamic load factors, it becomes 

k L P(t) - K R kx = K m ma , 

Recognizing that K R = K L * and dividing through by K L leads to 

p(t) - kx = (K m /K l ) ma 

The combined factor K m /K l is given directly in Table 6.3. The unanswered 
question is how K M /K L elastic and plastic values are to be "weighted" (a 
procedure suggested without explanation by at least two sources 22 ' 33 ). 

Two methods of weighting the factors were tested by use, first in a 
chart solution, discussed earlier in this section, then with the results 
compared to those from a full j3-method solution (triangular decay load 
pulses used in both solutions). The two "weighting" methods were: 

1. Apply elastic and plastic factors in accordance with relative 

areas under the elasto-plastic resistance plot (center figure 
page 11-4) -for example, using K e (=K M /K L elastic) and K p (= K M /K L 
plastic) and p = 7, the weighted factor K = (0.5/6.5)K + 

(6/6.5) K p . e 

2. Apply elastic and plastic factors in accordance with the relative 
displacements for elastic and plastic phases _ for example, using 
K e and K p with p = 7, the weighted factor K = (l/7) K e + 6/7) K p . 

The latter is recommended for use, because it is easier to apply and 
no discernible difference was found in the test applications of the two 
weighting methods. 


* Reference 34 (Equation 7.17, p. 148) of Volume 1. 
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The third change that may be readily applied as a part of "Step 7" 
is to use an ultimate moment capacity formula that takes into account 
rebound steel acting as compressive reinforcement. The change may be 
seen by comparing Equations 6-1 and 6-6; the effect is to allow a slighter 
smaller section depth d for a given clear span L. 

The first two changes described for consideration involve changing 
a preliminary design procedure that already has certain elements incorpo¬ 
rated in it that at least partially correct for the discrepancies dis¬ 
cussed ; for example, refer to the discussion of the calculation of the 
natural period of vibration T that follows Equation 6-5. On the other 
handj incorporating a final "refinement" such as the third change dis¬ 
cussed does not necessarily upset the assumptions made in the preliminary 
design procedure. For these reasons, the usefulness of the third change 
was tested first in developing the modified preliminary design procedure 
that follows. 


A preliminary design (first cut) might then be 
improved ("Step 7") by using instead the following modi- 
fied procedure: 

1. Given values for the following parameters: 

f c ( ksi ) ) f dc (ksi) J fcj y (ksi) ; p so (psi) ; 
t Q (sec); L (in.); p; b (in.) 

2. Obtain t Q0 (sec) from Figure 3-7.* 

3. Assume values for p, and T (sec); calculate 
too/T. 

4. Using Figure 6-1, page 11-5, read P m /q ; cal¬ 
culate q, using p m = p so . 

5. Calculate L/d, using the following modifica¬ 
tion of Equation 6.1: 

(L/d) = [8 pbf^/aq] 

[1- 0.554 pf dy /(0.85 f d ' c )] (G-2) 

obtaining L/d, then calculating d (in.). 

6. Using Equation 6-4, calculate T (sec), then 


Numerical 
Example: 

3 ; 3.75 ; 44; 15 
1.55; 255 ; .02; 

0.71 

6; .08; 

8.9 

.985 ; 

15.23 


19.788; 12.891 

.0839 ; 

8.46 


* Reference 2 of Volume 1. 

t Carrying more decimal places than justified, for demonstration purposes 
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Numerical 
Example: 


7. Using Figure B-10^ with p and t 0 o/T, read 
r/q (equals p7p); calculate p 7 (=p(r/q), 
but not less than .0025), then p (= 0.1/ 
(p-p 7 )); iterate this sequence until com¬ 
patible values are found for p ' and p. 


6; 8.46 


.0046; 6.5 


8. Using Equations 6-7 and -8, calculate 

"corrected" p go to include effect of slab 
weight. 


15+1.23 = 16.231 


9. 


Using corrected values for p and tQQ/T, repeat 
item 4 for new P m /q and, with new p m , calcu¬ 
late new q (psi) ; repeat item 5 for a new 
d(in.). Repeating items 6-8 should be 
unnecessary. 


6.5; 8.46 
1.00; 16.23; 
16.23 
13.31 


10. Web steel to meet diagonal tension require¬ 
ments should be calculated, and a check for 
pure shear strength should be made, using 
Equations 6-9 and -10, respectively. 


Another example with the same parameters except for changing f<( (5 ksi), 
f d ' c (6.25 ksi), f dy (78 ksi), and L (22’2”) led to a design d = 10.21". t 


A comparison with results from calculations using the final design 
procedure (Chapter 6) was then made at this point: 





dt 

L 

/ 

P 

P 

q 

Trial 1 

, modified 

design ' 

13.31" 

21’3" 

. 0046 

6. 5 

16.1 

Final 

design 


12.82" 

21'3" 

.0048 

6.6 

15.9 

Trial 2 

, mod ified 

design 

10.21" 

22’2" 

.0054 

6.9 

15.4 

Final 

design 


9.86" 

22'2" 

.0054 

6.9 

15.2 


# Reference 2 of Volume 1. 

t Carrying more decimal places than justified, for demonstration purposes. 
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Of the changes discussed earlier, the third one can be an inserted 
item in the modified procedure: 

9.1 Adjust d slightly to take into account the rebound steel 

acting as compressive reinforcement, by using the following 
form of Equation 6-6 (with d appearing on both sides of the 
equation, a cut-and-try solution is indicated; however, for 
most designs it should be sufficiently accurate to simply use 
the d from item 9, in the right-end term below): 

(L/d) 2 = [8 (p-p') bf dy /aq] [l - 0.554 (p-p / )f dy /(0.85 f^)] 

+ [8 p ; bf dy /aq] [l - d'/d]* (G-3) 

A new comparison of results was made, adding those where step 9.1 
had been applied: 



dt 

L 

/ 

P 

M- 

q 

Trial 1, modified design 

13.31" 

21’ 3" 

/ 

.0046 

6.5 

16.1 

Modified design with 
step 9.1 

12.92" 

21’ 3" 

.0046 

6.5 

16.1 

Final design procedure 

12.82" 

21’3" 

.0048 

6.6 

15.9 

Trial 2, modified design 

10.21" 

22'2" 

.0054 

6.9 

15.4 

Modified design with 
step 9.1 

9.91" 

22 * 2" 

.0054 

6.9 

15.4 

Final design procedure 

9.86" 

22’2" 

.0054 

6.9 

15.2 


The close agreement of results from the modified preliminary design 
procedure with step 9.1 inserted and the final design procedure results 
made it appear that further work on the modified form of the preliminary 
design procedure was unwarranted. If further design refinement is thought 
to be warranted, the final design procedure (Chapter 6) may be used. 


* As in final design procedure (Chapter 6), d' assumed as one in. 
t Carrying more decimal places than justified, for demonstration purposes. 
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Tabular Data on Typical Designs 


Table G.l contains tabular data on the typical designs, made using 
the final design procedure of Chapter 6, for preparation of Figures 6-3 
and 6-6. 

Computer Program Used for Typical Designs 

Table G-2 shows a listing of the computer program used in calculating 
the typical designs for which data are shown in Table G.l. The program¬ 
ming language is Dartmouth BASIC as used by one of the largest commercial 
time-sharing organizations. For ease of handling data and plotting 
results, d was assumed and L calculated (rounded-off to a whole number, 
in inches). In retrospect, the beta method* time interval of one msec 
was smaller than necessaryj one might calculate a time interval based on 
the estimated T in the program, using 0.05 T or 0.1 T (accuracy has been 
estimated to be within about 1% or 3.6%, respectively, of results obtained 
from a very small time interval, 0.001 T).t Also, the time interval could 
be increased considerably once the plastic phase is reached. 


* Reference 21 of Volume 1. 
t Reference 20 (p. 33) of Volume 1. 
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Table G.l 


DATA FOR TYPICAL DESIGNS - ONE-WAY SLABS, SIMPLY SUPPORTED 

p = .02 


f'/f 

C dy 

, d 

L 

p' 

P v 

M 

max 

q 

X 

e 

X 

m 


t 

m 

P 

m 

(ksi) 

(in.) 

(in.) 

(%) 

(%) 

(ft-kip) 

(psi) 

(in.) 

(in.) 


(msec) 

(psi) 

3/52 

6 

126 

.38 

.97 

31.8 

16.0 

.896 

5.6 

6.25 

89 

15.59 


10 

214 

.45 

1.04 

90.4 

15.8 

1.58 

10.3 

6.54 

148 

15.96 


14 

304' 

.52 

1.10 

179.9 

15.6 

2.31 

15.6 

6.76 

204 

16.32 


18 

397 

.61 

1.18 

301.8 

15.3 

3.125 

22.5 

7.20 

265 

16.66 


21 

471 

.71 

1.25 

416.4 

15.0 

3.82 

29.2 

7.65 

315 

16.91 

4/52 

6 

128 

.38 

.77 

32.9 

16.0 

.869 

5.5 

6.28 

88 

15.59 


10 

217 

.45 

.83 

93.2 

15.8 

1.52 

10.1 

6.60 

147 

15.96 


14 

308 

.53 

.88 

185.1 

15.6 

2.22 

15.2 

6.86 

202 

16.32 


18 

402 

.62 

.94 

309.4 

15.3 

2.97 

21.8 

7.35 

262 

16.66 


21 

476 

.73 

1.01 

425.6 

15.0 

3.615 

28.4 

7.85 

313 

16.91 

5/52 

6 

130 

.38 

.62 

33.6 

16.0 

.845 

5.3 

6.33 

88 

15.59 


10 

218 

.45 

.67 

94.9 

15.9 

1.465 

9. 3 

6.36 

140 

15.96 


14 

311 

.53 

.72 

188.0 

15.6 

2.137 

14.9 

6.95 

201 

16.32 


18 

405 

.63 

,77 

313.9 

15.3 

2.858 

21. 3 

7.47 

261 

16.66 


21 

479 

.74 

.82 

430.8 

15.0 

3.457 

27.7 

8.02 

311 

16.91 

3/72 

6 

146 

.41 

1.15 

41.9 

15.7 

1.58 

10.3 

6.48 

118 

15.59 


9 

223 

.48 

1.23 

96.7 

15.5 

2.53 

16.3 

6.45 

173 

15.87 


12 

305 

.56 

1.30 

175.6 

15.1 

3.62 

25.2 

6.95 

237 

16.14 


15 

395 

.71 

1.43 

282.6 

14.5 

4.995 

39.3 

7.86 

319 

16.41 


18 

501 

.99 

1.65 

424.0 

13.5 

6.99 

68.0 

9.74 

440 

16.66 

4/72 

6 

150 

.42 

.99 

44.0 

15.7 

1.58 

10.2 

6.48 

118 

15.59 


9 

229 

.49 

1.05 

101.0 

15.4 

2.52 

17.0 

6.77 

178 

15.87 


12 

311 

.58 

1.12 

182.7 

15.1 

3.55 

24.8 

7.00 

237 

16. 14 


15 

401 

.73 

1.21 

291.6 

14.5 

4.82 

38.5 

8.00 

317 

16.41 


18 

499 

.95 

1.35 

430.1 

13.8 

6.37 

58.4 

9.16 

410 

16.66 

5/72 

6 

152 

.42 

.85 

45.2 

15.7 

1.56 

10.1 

6.50 

118 

15.59 


9 

231 

.49 

.90 

103.5 

15.5 

2.45 

16.0 

6.51 

172 

15.87 


12 

314 

.58 

.96 

186.7 

15.1 

3.455 

24.4 

7.07 

235 

16.14 


15 

405 

.74 

1.05 

296.7 

14.5 

4.655 

37.8 

8.13 

315 

16.41 


18 

502 

.97 

1.17 

435.4 

13.8 

6.095 

57.2 

9.39 

408 

16.66 
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Table G.l (continued) 


p = .015 


f'/f 

d 

L 

P / 

P 

M 

q 

X 

X 

u 

t 

P 

c dy 




V 

max 


e 

m 


m 

m 

(ksi) 

(in.) 

(in.) 

(%) 

(%) 

(ft-kip) 

(psi) 

(in.) 

(in.) 


(msec) 

(psi) 

3/52 

6 

112 

.27 

.76 

24.7 

15.7 

.664 

5.3 

8.05 

93 

15.59 


11 

210 

.34 

.83 

84.8 

15.4 

1.296 

11.4 

8.79 

170 

16.05 


16 

311 

.42 

.90 

182.2 

15.1 

1.98 

18.6 

9.37 

244 

16.495 


20 

395 

.53 

.99 

288.8 

14.8 

2.60 

25.9 

9.95 

306 

16.83 


24 

480 

.63 

1.07 

420.9 

14.6 

3.245 

32.6 

10.05 

358 

17.15 

4/52 

6 

114 

.27 

.57 

25.3 

15.7 

.64 

5.2 

8.12 

92 

15.59 


11 

212 

.34 

.63 

86.7 

15.4 

1.24 

11.1 

8.93 

169 

16.05 


16 

314 

.43 

.69 

185.8 

15.1 

1.89 

18.1 

9.58 

242 

16.495 


20 

397 

.53 

.76 

293.2 

14.9 

2.44 

24.0 

9.82 

296 

16.83 


24 

482 

.63 

.82 

426.1 

14.7 

3.025 

30. 3 

10.02 

347 

17.15 

5/52 

6 

114 

.27 

.50 

25.7 

15.7 

.62 

5.1 

8.19 

92 

15.59 


11 

213 

.34 

. 50 

87.8 

15.5 

1.19 

10.2 

8.55 

162 

16.05 


16 

314 

.43 

.53 

187.8 

15.2 

1.81 

16.8 

9.29 

234 

16.495 


20 

399 

.54 

. 59 

296.0 

14.9 

2.34 

23.5 

10.03 

294 

16.83 


24 

482 

.63 

.64 

429.1 

14.8 

2.87 

28.3 

9.88 

337 

17. 15 

3/72 

6 

131 

.30 

.98 

32.9 

15.4 

1.20 

10.0 

8.33 

124 

15.59 


10 

224 

.38 

1.06 

94.0 

15.0 

2.176 

19.2 

8.82 

205 

15.96 


13 

298 

.47 

1.15 

162.0 

14.6 

3.023 

28.7 

9.51 

272 

16.23 


16 

376 

.59 

1.25 

250.6 

14.2 

3.99 

40.1 

10.04 

340 

16.495 


20 

482 

.74 

1.38 

400.5 

13.8 

5.37 

54.5 

10.14 

418 

16.83 

4/72 

6 

133 

.30 

.80 

34.1 

15.4 

1.18 

9.9 

8.37 

124 

15.59 


10 

227 

.38 

.87 

96.8 

15.0 

2.15 

18. 9 

8.92 

204 

15.96 


13 

302 

.49 

.95 

166.3 

14.6 

2.92 

28.2 

9.67 

271 

16.23 


16 

378 

.59 

1.02 

255.4 

14.3 

3.77 

37.3 

9.90 

329 

16.495 


20 

483 

.74 

1.12 

405.7 

13.9 

5.013 

50.9 

10.15 

406 

16.83 

5/72 

6 

134 

.30 

.67 

34.8 

15.4 

1.159 

9.8 

8.42 

123 

15.59 


10 

229 

.39 

.74 

98.5 

15.0 

2.063 

18.6 

9.01 

203 

15.96 


13 

304 

.49 

.80 

168.8 

14.6 

2.825 

27.7 

9.82 

270 

16.23 


16 

380 

.60 

.86 

258.4 

14.3 

3.632 

36.7 

10.10 

328 

16.495 


20 

483 

.74 

.95 

408.9 

14.0 

4.756 

47.8 

10.05 

395 

16.83 




Table G.1 (continued) 


p = .01 


f'/f 
c dy 

d 

L 

P' 

P 

V 

M 

max 

q 

X 

e 

X 

m 


t 

m 

P 

m 

(ksi) 

(in.) 

(in.) 

(%) 

(%) 

(it-kip) 

(psi) 

(in.) 

(in.) 


(msec) 

(psi) 

3/52 

6 

93 

.25 

.53 

17.0 

15.6 

.421 

4.10 

9.79 

87 

15.59 


' 12 

189 

.25 

.57 

69.5 

15.5 

.887 

8.70 

9.77 

160 

16.14 


18 

286 

.27 

.60 

157.8 

15.4 

1.363 

13.5 

9.88 

224 

16.66 


24 

385 

.34 

.66 

283.5 

15.3 

1.87 

18.5 

9.92 

282 

17.15 


30 

485 

.42 

.73 

447.2 

15.2 

2.396 

23.9 

9.96 

337 

17.626 

4/52 

6 

94 

.25 

.50 

17.2 

15.6 

.402 

4.0 

9.91 

86 

15.59 


12 

191 

.25 

.50 

70.5 

15.5 

.847 

8.4 

9.94 

158 

16.14 


18 

288 

.28 

.50 

159.5 

15.4 

1.301 

13.1 

10.08 

222 

16.664 


24 

387 

.35 

.50 

286.6 

15.3 

1.774 

18.0 

10.17 

280 

17.15 


30 

485 

.42 

.52 

450.8 

15.3 

2.247 

22.2 

9.89 

327 

17.626 

5/52 

6 

94 

.25 

.26 

17.4 

15.6 

.389 

3.9 

10.00 

85 

15.59 


12 

191 

.25 

.26 

71.0 

15.5 

.818 

8.2 

10.06 

157 

16.144 


18 

289 

.28 

.50 

161.1 

15.4 

1.255 

12.9 

10.24 

221 

16.664 


24 

387 

.34 

.50 

288.2 

15.4 

1.694 

16.8 

9.94 

271 

17.15 


30 

487 

.42 

.50 

452.9 

15.3 

2.157 

21.8 

10.10 

326 

17.626 

3/72 

6 

110 

.25 

.77 

23.0 

15.3 

.788 

8.0 

10.21 

117 

15.59 


11 

205 

.25 

.80 

79.0 

15.1 

1.527 

15.0 

9.79 

195 

16.054 


15 

283 

.30 

.85 

148.8 

14.9 

2.165 

21.2 

9.81 

254 

16.408 


20 

384 

.40 

.95 

269.3 

14.6 

3.051 

30.6 

10.01 

327 

16.830 


25 

488 

.49 

1.03 

426.5 

14.3 

4.00 

40.7 

10.16 

398 

17.228 

. 4/72 

6 

111 

.25 

.60 

23.5 

15.3 

.762 

7.9 

10.34 

117 

15.59 


11 

207 

.25 

.63 

80.5 

15.1 

1.474 

14.6 

9.93 

194 

16.054 


15 

285 

.31 

.68 

151.4 

14.9 

2.083 

20.8 

9.99 

252 

16.408 


20 

385 

.40 

.75 

272.5 

14.7 

2.885 

28.5 

9.87 

317 

16.83 


25 

489 

.49 

.82 

430.2 

14.4 

3.758 

38.1 

10.13 

387 

17.228 

5/72 

6 

111 

.25 

.50 

23.7 

15.4 

.736 

7.1 

9.71 

110 

15.59 


11 

208 

.25 

.50 

81.5 

15.1 

1.434 

14.4 

10.04 

193 

16.054 


15 

286 

.31 

.54 

152.9 

14.9 

2.018 

20.5 

10.14 

251 

16.408 


20 

386 

.40 

.60 

274.5 

14.7 

2.782 

28.0 

10.06 

316 

16.83 


25 

488 

.49 

.66 

432.4 

14.5 

3.584 

35.8 

10.0 

377 

17.228 


26 

508 

.50 

.67 

468.2 

14.5 

3.736 

36.8 

9.85 

385 

17.305 

5/52 

12.2 

195 

.25 

.26 

73.4 

15.5 

.832 

8.4 

10.0 

159 

16.16 


12.4 

198 

.25 

.26 

75.9 

15. 5 

.847 

8.5 

10.0 

161 

16.18 


12.6 

201 

.25 

.26 

78.4 

15.5 

.861 

8.6 

10.0 

163 

16.20 


12.8 

204 

.25 

.50 

80.9 

15.5 

.875 

8.8 

10.01 

165 

16.215 


12.9 

206 

.25 

.50 

82.2 

15.5 

.882 

8.8 

10.0 

166 

16.22 


13.0 

208 

.25 

.50 

83.5 

15.5 

.889 

8.9 

10.0 

167 

16.23 


14.0 

224 

.25 

.50 

97.0 

15.5 

.961 

9.5 

9.93 

177 

16.32 
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Table G.l (concluded) 


p = .005 


£'/f 

d 

L 

P' 

P 

M 

q 

X 

X 

M- 

t 

P 

c dy 




V 

max 


e 

m 


m 

ni 

(ksi) 

(in.) 

(in.) 

(%) 

(%) 

(f t-kip) 

(psi) 

(in.) 

(in.) 


(msec) 

(psi) 

3/52 

6 

66 

.25 

.26 

8.5 

15.8 

.176 

1.8 

10.29 

60 

15.59 


■ 16 

179 

.25 

.28 

63.7 

16.0 

.518 

5.1 

9.82 

142 

16.495 


25 

280 

.25 

.29 

157.4 

16.1 

.824 

8.3 

10.02 

204 

17.228 


36 

402 

.25 

.30 

328.5 

16.3 

1.187 

12.0 

10.08 

270 

18.15 


43 

477 

.25 

.30 

469.7 

16.5 

1.407 

14.0 

9.94 

306 

18.764 


47 

520 

.25 

.31 

561.7 

16.6 

1.532 

15.2 

9.93 

326 

19.114 

4/52 

6 

66 

.25 

.26 

8.5 

15.8 

. 168 

1.7 

10.38 

59 

15.59 


16 

179 

.25 

.26 

63.9 

16.0 

.494 

4.9 

9.97 

140 

16.495 


26 

291 

.25 

.26 

170.9 

16.2 

.813 

8.0 

9.80 

203 

17.305 


36 

401 

.25 

.26 

329.4 

16.4 

1.125 

11.1 

9.89 

261 

18.151 


43 

478 

.25 

.26 

471.1 

16.5 

1.342 

13.7 

10.18 

304 

18.764 


47 

519 

.25 

.26 

563.4 

16.7 

1.452 

14.3 

9.83 

317 

19.114 

5/52 

6 

66 

.25 

.26 

8.5 

15.8 

.162 

1.7 

10.43 

58 

15.59 


16 

179 

.25 

.26 

64.0 

16.0 

.478 

4.8 

10.09 

139 

16.495 


26 

291 

.25 

.26 

171.2 

16.2 

.787 

7.8 

9.93 

202 

17.305 


36 

401 

.25 

.26 

330.0 

16.4 

1.088 

10.9 

10.05 

260 

18.151 


43 

477 

.25 

.26 

471.9 

16.6 

1.290 

12.9 

9.97 

296 

18.764 


47 

520 

.25 

.26 

564.3 

16.7 

1.404 

14.0 

10.00 

316 

19.114 

3/72 

6 

77 

.25 

.50 

11.7 

15.7 

. 336 

3.2 

9.62 

77 

15.59 


14 

185 

.25 

.54 

66.9 

15.6 

.876 

8.5 

9.74 

168 

16.321 


22 

294 

.25 

.55 

167.6 

15.5 

1.423 

14.2 

9.97 

245 

16.992 


31 

417 

.25 

.56 

335.0 

15.4 

2.046 

20.8 

10.16 

323 

17.713 


36 

485 

.25 

.56 

452.8 

15.4 

2.387 

24.2 

10.15 

362 

18.15 


40 

539 

.26 

.58 

560.2 

15.4 

2.664 

27.0 

10.15 

393 

18.501 

4/72 

6 

77 

.25 

.32 

11.7 

15.7 

. 322 

3.1 

9.73 

76 

15.59 


14 

186 

.25 

. 38 

67.2 

15.6 

.837 

8.3 

9.92 

166 

16.321 


23 

308 

.25 

.39 

184.1 

15.5 

1.425 

14.5 

10.15 

251 

17.072 


32 

430 

.25 

.40 

358.6 

15.5 

2.008 

20.1 

10.01 

321 

17.801 


36 

484 

.25 

.40 

454.6 

15.5 

2.267 

22.8 

10.04 

353 

18.15 


40 

539 

.26 

.41 

562.3 

15.5 

2.528 

25.5 

10.07 

383 

18.501 

5/72 

6 

77 

.25 

.26 

11.8 

15.7 

. 312 

3.1 

9.80 

75 

15.59 


14 

186 

.25 

.27 

67.4 

15.6 

.811 

8.1 

10.04 

165 

16.321 


22 

294 

.25 

.28 

168.6 

15.6 

1.309 

12.9 

9.85 

235 

16.992 


31 

417 

.25 

.29 

337.1 

15. 5 

1.882 

19.1 

10.17 

312 

17.713 


36 

483 

.25 

.29 

455.7 

15.6 

2.181 

21.5 

9.86 

344 

18.151 


40 

538 

.26 

. 30 

563. 5 

15.6 

2.4 32 

24.1 

9.92 

374 

18.501 
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NOTATION 


See Notation section at end of Chapter 6. 
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